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EXPERIMENTAL  AND  SYSTEMS  DESIGN,  APPARATUS  AND  EXPERIMENTAL 

PROCEDURE  CHAPTER 


THEORETICAL  BACKGROUND 

Evaluation  of  candidate  human  protective  devices  against  impact  acceleration 
'requires  human  simulation  by  such  analogs  as  cadavers,  dummies,  mathematical 
models  and  various  primates.  The  human  analog  is  required  because  man  cannot 
tolerate  the  extreme  impact  forces  against  which  such  systems  are  designed  to  afford 

protection. 

**  •  *  \ 

All  human  analog  testing  and  data  derived  therefrom  suffers  from  one  principal 
defect.  Its  correspondence  to  living  man  is  not  quantitatively  demonstrable,  largely 
because  of  the  almost  complete  absence  of  quantitative  human  dynamic  response 
data,  and  injury  criteria  pertaining  thereto.  As  a  consequence,  dummies  are  designed 
to  reproduce  static  rather  than  dynamic  anthropometric  measurements. 

The  use  of  living  volunteers  to  test  protective  systems  is  dangerous,  expensive 
and  difficult  to  interpret.  If  a  subject  is  not  used  as  his  own  control,  variance 
between  subjects  makes  comparisons  difficult.  Adjustments  of  certain  protective 
restraint  systems  are  a  very  large  factor  in  their  ability  to  protect,  but  are  difficult 
to  reproduce  validly  for  crash  tests  even  for  the  same  subject. 

Development  of  a  valid  human  analog  would  result  in  the  capability  to  reproduce 
man's  inertiai  dynamic  response  in  a  standard  lest  object.  Such  an  analog  does  not 
exist  today.  Its  development  requires  determination  of  the  following: 

1 .  Mass,  center  of  mass  location  and  principal  moments  of  inertia  of  the 
important  body  segments.  (Some  of  these  data  are  available) 

2.  Functional  relationships  between  those  body  segments  in  the  dynamic 
environment.  The  goal  is  to  be  able  to  predict  dynamic  response  to  a  given  segment, 
given  the  input  accelerations  to  that  segment. 

The  following  logic  was  used  to  attempt  to  measure  response  data  which  would 
be  useable  for  analog  design  and  validation: 

1.  Man's  dynamic  response  to  small  experimentally  produced  increments 
of  sled  acceleration  (1)  peak,  (2)  rate  of  onset,  (3)  duration  at  peak  and  (4)  vector 
directions, is  safely  measureabSe  up  to  the  point  of  voluntary  tolerance. 

2.  An  analog  which  accurately  reproduces  man's  dynamic  responses  to 
simple,  experimentally  produced  individual  input  accelerations  will  yield  a  valid 
total  dynamic  response  to  complex  input  accelerations  such  as  those  of  a  crash. 

It  was  believed  that  such  responses  could  only  be  measured  by  instrumenting 
the  living  human  being. 
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Successful  measurement,  instrumentally,  of  human  head  and  neck  response 
requires  the  ability  to  mount  transducers  on  the  man  to  achieve  (a)  reproducibility 
of  placement  between  runs;  (b)  lack  of  spurious  movement  of  the  transducer  relative 
to  the  man;  (c)  results  which  could  be  related  to  the  human  population  rather  than 
to  the  individual  alone;  (d)  location  of  the  transducers  in  the  plane  of  expected 
motion;  and  (o)  repeated  runs  on  a  single  instrumented  individual  to  determine  the 
repeatability  of  any  data  collected,  as  well  as  to  eliminate  inter-sub{ect  variation. 

The  majority  of  pre-1967  instrumentation  effort  has  usually  been  directed 
toward  the  vehicle.  Instrumentation  of  Stapp's  experiments  was  limited  by  the  state- 
of-the-art  of  inertial  measurement  transducers,  especially  accelerometers  and  rate 
gyroscopes.  In  those  experiments,  Stapp  noted  that  head-mounted  transducers  as 
mountea  on  a  helmet  or  a  "bite-plate"  could  not  give  useable  results.  (1) 

If  instrumentation  can  be  firmly  mounted  on  the  man,  changes  in  his  dynamic 
response  due  to  changes  in  any  of  tne  input  as  measured  on  the  subject  can  be 
precisely  determined.  These  changes  with  increasing  input  accelerations  can  be 
exactly  measured  by  using  each  subject  as  hie  own  control .  Differences  between 
instrumented  subjects  can  be  determined  by  comparisons  of  dynamic  response 
obtained  from  identical  input  acceleration  pulses  applied  to  them. 

A  joint  Army-Navy  project  was  initiated  in  1966  at  the  Naval  Aerospace  Medical 
Institute,  Pensacola,  Florida  to  make  these  measurements.  The  last  run  was  completed 
on  27  April  1969. 

Dynamic  response  of  the  head  and  neck  to  experimentally  applied  impact 
acceleration  was  the  first  subject  of  study.  Previous  studies  (1,  2)  have  demonstrated 
the  difficulties  inherent  in  measuring  accelerations  and  displacements  of  the  head 
and  neck,  as  referenced  to  the  seat.  One  difficulty  has  been  in  the  means  of 
generating  the  input  acceleration.  Many  studies  have  been  performed  on  acceleration 
facilities  in  whicn  an  initial  acceleration  was  given,  followed  by  a  "constant"  velocity 
phase  (actually  a  deceleration  of  less  than  0.5  G),  followed  by  a  terminal  shaped 
deceleration  pulse. 

One  disadvantage  of  such  an  accelerator  is  the  necessity  for  establishing  a 
condition  of  zero  dynamic  response  prior  to  entry  into  the  deceleration  phase.  If  the 
head  and  neck  were  responding  to  the  initial  acceleration  or  subsequent  velocity 
deccy  at  the  time  of  initiation  of  the  input  acceleration,  an  artifact  would  be 
introduced  into  the  measured  response  to  the  deceleration  pulse. 

Another  disadvantage  is  that  physiological  data  collected  during  an  experiment 
cannot  be  related  to  a  single  pulse  but  instead  are  related  to  a  series  of  pulses.  It 
would  appear  easier  to  determine  the  effect  of  a  single  pulse  if  the  experimental 
pulse  were  induced  with  the  subject  at  both  inertial  and  physiological  "rest". 

Another  difficulty  lies  in  the  complexity  of  the  forces  acting  simultaneously  on 
the  head.  A  man  seated  on  a  sled  and  restrained  by  pelvic  and  shoulder  harness 
when  subjected  to  a  -G^  acceleration  transmitted  to  his  body  by  an  interaction  of 

his  body  with  the  seat  back  and  restraint  harness .  (2)  The  sled  moves  away  from  the 
pelvis  in  the  direction  of  acceleration,  but  the  pelvic  restraint  which  couples  the 
man  to  the  seat  transmits  the  acceleration  to  the  pelvis  and  thus  to  the  torso.  The 
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torso  then  attempts  to  rotate  around  the  restrained  pelvis  in  the  mid-sagittal  plane  but 
is  prevented  from  doing  so  by  the  shoulder  restraint  harness.  Both  the  pelvic  and 
shoulder  harnesses  stretch  in  response  to  to  the  dynamic  load  of  the  subject.  A 
reversal  of  torso  trajectory  then  occurs,  resulting  in  the  torso  being  forced  down  and 
back.  Concurrently,  the  torso  is  being  moved  through  space  with  the  acceleration 
vehicle. 

The  head  and  neck  also  respond  by  rotating  through  space  about  a  center  of 
rotation  in  the  neck  or  neck-torso  juction  area.  Schulman  et  al.  have  shown  that 
th  ’  neck  is  capable  of  considerable  stretch.  (2)  Therefore,  displacements  of  the  Head 
as  measured  from  a  seat  reference  point  are  difficult  to  interpret  due  to  the  complexity 
of  motion  of  the  entire  sled-torso  -  head  system. 

The  purposes  of  the  present  study  were  fourfold:  (1)  to  measure  precisely  the 
complete  in  >ut  acceleration  to  the  head  and  neck  measured  at  the  first  thoracic 
vertebra;  (2)  to  measure  precisely  the  dynamic  response  of  the  head  and  neck  to  the 
input  acceleration;  (3)  to  develop  a  method  of  obtaining  the  data  in  such  a  format 
that  automatic  data  processing  may  be  used  extensively;  (4)  to  develop  and  validate 
a  general  method  for  the  experimental  measurement  of  the  bioengineering  characteris¬ 
tics  of  the  human  body  with  such  precision,  accuracy  and  repeatability  that  a  mathe¬ 
matical  model  or  other  analog  of  the  human  dynamic  response  to  impact  acceleration 
can  be  constructed.  Such  a  model  has  many  important  applications,  but  previous 
modeling  attempts  have  been  handicapped  by  a  paucity  or  data  of  necessary  quality, 
resulting  in  models  of  limited  utility. 

EXPERIMENTAL  DESIGN 

The  following  procedures  were  adopted  in  order  to  achieve  the  stated  goals: 

1 .  Displacement  was  measured  by  double  integration  with  respect  to  time 
of  acceleration  measured  on  the  head  and  Tj  spine  by  accelerometers  and  rate  gyros 
rigidly  mounted  on  the  man  in  completely  reproducible  locations  whose  relationship 
to  bony  anatomy  was  measured  by  x-ray  anthropometry, 

2.  A  second  and  independent  means  of  measuring  displacement  consisted 
of  sled-mounted,  pin-registered,  high-speed,  high-precision,  16  mm  cine  cameras. 
Data  derived  from  photographic  data  reauction  were  used  to  validate  transducer- 
derived  data. 

3.  All  data  were  collected  in  a  format  compatible  with  automatic  data 
processing  by  digital  computer. 


the  track. 


4.  All  data  of  primary  interest  were  collected  on  the  acceleration  end  of 


In  the  design  of  these  experiments,  certain  simplifying  assumptions  were  made, 
the  first  was  that  all  accelerations  acting  on  the  head  and  neck  of  the  seated  pelvic- 
and  torso-restrained  individual  subjected  to  a  -Gx  acceleration  would  be  detectable 
at  the  first  thoracic  vertebra  (Ti)  or  base  of  the  neck  since  all  accelerations  acting 
on  the  head  and  neck  must  be  transmitted  through  the  vertebral  column.  Only  those 
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additional  accelerations  transmitted  by  the  soft  tissues  of  the  neck  could  be  a  source 
of  error.  Therefore,  the  acceleration  measured  at  T |  was  presumed  to  serve  as  input 
information  te  the  head  and  neck.  The  second  assumption  was  that  the  head  is  a 
rigid  body,  and  the  third,  that  all  significant  head  and  neck  motion  would  be  coplanar 
with  the  plane  of  induced  acceleration,  i.e.,  the  mid-sagittal  plane. 

In  order  to  describe  clearly  the  complex  motion  of  the  head  and  neck,  all  motion 
was  considered  as  occurring  in  four  two-dimensional  coordinate  systems,  as  described 
in  the  Theoretical  Kinematics  Chapter: 

1 .  The  laboratory  coordinate  system  was  selected  as  the  position  occupied 
by  the  sled  coordinate  system  immediately  prior  to  the  onset  of  acceleration  to  the 
sled. 

2.  The  sled  coordinate  system  was  selected  as  an  arbitrary  photographic 
reference  system  fixed  to  the  chair  back  with  a  known  orientation  to  gravity. 

3.  The  T]  coordinate  system  was  established  as  a  two-point  photographic 
reference  system  on  the  Tj  mount,  as  shown  in  Figure  DP-9  with  axes  directions 
further  defined  electronically  as  described  under  the  section  on  Transducer  System. 

4.  The  head  coordinate  system  was  established  by  photographic  targets 
on  the  head  (bregma)  mount  and  the  mouth  mount  as  shown  in  Figure  DP-9  with 
axes  directions  further  defined  electronically  as  described  under  the  section  on 
Transducer  System. 

Two  redundant  systems  tracked  the  coordinate  systems.  The  first  system  used  an 
inertial  reference  platform  mounted  on  each  coordinate  system  of  the  head  and  Tj 
with  miniature  transducers  to  measure  acceleration  and  angular  velocity.  In  the 
second,  high-precision  cinematography  was  used  to  measure  displacements  of  the 
coordinate  systems.  Velocities  and  accelerations  also  were  determined  from  the 
measured  displacements  by  differentiation  with  respect  to  time. 

Precise  anatomical  displacement  information  can  be  obtained  from  the  measured 
coordinate  system  displacements  if  several  requirements  are  met:  (1)  The  transducers 
and  the  photographic  targets  must  be  uniquely  fixed  to  each  individual;  this  was 
satisfied  by  the  construction  of  mouth,  head  and  spine  mounts  that  can  be  placed  only 
in  one  position  on  each  man.  (2)  The  position  of  each  photographic  target  and  each 
transducer  must  be  recorded  for  each  individual;  this  was  accomplished  by  x-ray  and 
photographic  anthropometry  of  each  subject's  head  and  neck  with  the  mounts  in  place. 
(3)  The  position  of  each  transducer  and  target  in  space  must  be  determined  fust  prior 
to  the  onset  of  acceleration;  this  requirement  was  met  by  data  from  the  lateral  and 
posterior  cameras  and  by  the  x-ray  anthropometry.  Once  all  anatomical  displacements 
were  determined,  a  comparison  was  made  to  validate  the  results  obtained  by  the  use 
of  either  the  photographic  or  transducer  systems. 

Each  system  used  to  accomplish  the  entire  experiment  was  carefully  integrated 
with  each  other.  The  equipment  and  techniques  will  be  described  system  by  system. 


4 


SYSTEMS  DESIGN  AND  APPARATUS 
Transducer  System 

Transducer  systems  forming  Inertia!  reference  platforms  were  established  at  three 
locations:  The  spine  (T])  system  was  positioned  over  the  posterior  spinous  process 
ofT1  and  consisted  of  two  accelerometers  mounted  with  tneir  sensitive  axes  at  right 
angles  in  the  mid-sagittal  plane,  and  a  rate  gyroscope  with  the  sensitive  axis 
perpendicular  to  the  mid-sagittal  plane  of  the  subject.  The  force-balance  accelero¬ 
meters  weighed  14.2  grams  each,  while  the  rate  gyroscopes  weighed  approximately 
94  grams.  The  head  system  was  positioned  over  the  back  of  the  shaven  nead  of  the 
subject  and  consisted  of  two  accelerometers  mounted  with  their  sensitive  axes  at  right 
angles  to  one  another  in  the  mid-sagittal  plane.  The  mouth  system  was  positioned  by 
a  bite-block  and  consisted  of  a  duplicate  set  of  the  transducers  mounted  at  Tj .  The 
transducer  placement  is  illustrated  in  Figure  1. 


ACCELEROMETERS  — ► 

RATE  GYROS  O 

FIGURE  1 

The  head  and  mouth  system  outputs  were  combined  to  form  a  single  head  output 
as  described  under  OUTVAR  in  the  Data  Processing  Chapter.  The  outputs  of  the  head 
system  and  those  of  the  Tj  system  were  used  for  tracking,  each  with  respect  to  the 
laboratory  coordinate  system.  The  displacements  of  Tj  were  then  subtracted  from  those 
of  the  head  to  produce  differential  head  displacement  relative  to  Ti . 
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The  selection  of  transducers  and  the  remainder  of  the  electronic  data  system  was 
made  to  limit  the  error  of  the  overall  system  to  0.5%,  This  criterion  dictated  the 
choice  of  linear  accelerometers  that  had  a  maximum  nonlinearity  error  of  0. 1%  full 
scale  and  ±  5000  deg/sec  rate  gyros  with  a  maximum  error  of  0,5%  full  scale  within 
the  range  of  expected  angular  rates.  The  accelerometers  responded  as  a  second  order 
system  with  natural  frequency  ranging  from  112  Hz  to  165  Hz  and  damping  ratios  from 
0,5  to  0,8,  The  rate  gyroscopes  had  a  natural  frequency  of  64  Hz  ±  4  Hz  with 
damping  of  0,6  ±0,1, 

Transducer  outputs  were  hard  wired  to  operational  amplifiers  at  the  analogue  tape 
recorder.  Time-locking  of  the  transducer  recording  with  the  photographic  data  system 
was  accomplished  by  an  electrical  pulse  (T  pulse),  recorded  on  one  of  the  tape 
channels,  that  reset  the  conti  nuously-counfir  g  camera  timers  to  zero. 

The  time  axis  of  the  experiment  was  established  by  the  100  kHz  crystal  oscillator 
in  the  tape  recorder  phase-lock  system.  All  transducer  outputs  recorded  on  the 
analogue  tape  were  phased  to  this  time  reference. 


A  mount  consisting  of  three  modules  was  required  at  each  transducer  system 
location.  One  module  was  required  to  obtain  conformity  to  the  subject's  anatomy 
in  a  completely  repeatable  way;  the  second  was  required  to  maintain  the  transducers 
in  a  known  standard  configuration;  and  the  third  was  required  to  mate  the  individual 
anatomical  modules  to  the  common  transducer  modules. 

The  mouth  mount  and  head  mount  were  designed  to  compress  the  skull  between 
them,  thus  preventing  motion  relative  to  the  head  during  acceleration.  The  connection 
between  the  head  and  mouth  mounts  was  established  and  maintained  under  considerable 
constant  pressure  by  a  harness. 

A  pressure  mold  of  the  shaven  bregma  of  the  subject  was  used  to  construct  the 
head  anatomical  module.  The  mouth  anatomical  module  was  constructed  by  making 
a  metal  casting  of  the  teeth,  hard  palate,  ami  upper  jaw. 

The  Tj  anatomical  module  was  constructed  by  making  a  pressure  mold  of  the 
spinal  column  of  the  subject  from  Cg  through  T3  and  cutting  it  to  fit  anatomically. 

It  was  secured  by  a  harness  that  mamtained  the  mount  on  the  selected  site  under 
considerable  pressure.  The  acceleration  pulse  drove  the  mount  against  its  spinal 
anchor  point,  thus  improving  stability  during  the  dynamic  event. 

The  three  mounts  were  thus  installed  in  a  precisely  repeatable  location  for  each 
run.  After  the  anatomical  mounts  were  constructed  for  each  subject,  the  transducer 
modules  were  installed  on  the  center  line  so  that  each  was  coplanar  with  the  other 
two.  Once  these  relationships  were  established,  changing  the  mounts  from  one 
subject  to  another  was  a  brief  task. 

Immediately  prior  to  each  run,  the  subject  was  positioned  so  that  all  transducer 
systems  were  coplanar  in  the  mid-sagittal  plane  which  was  the  plane  of  expected 
induced  movement. 
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Motion  of  the  mounts  relative  to  the  anatomical  mounting  site  would  produce  \ 
suprious  displacements  and  accelerations  as  measured  both  instrumental!/  and 
photographically.  It  was  thus  necessary  to  determine  the  presence  or  absence  of  such 
relative  motion.  Motions  of  the  Tj  mount  relative  to  the  spine  were  determined 
photographically  by  the  posterior  and  lateral  cameras.  While  some  motion  was 
detected,  particularly  on  early  high-G  runs,  later  mount  improvements  reduced  these 
movements  to  relative  insignificance.  Motions  of  the  head  mount  relative  to  the 
mouth  mount  were  determined  photographically  on  the  frame  of  film  taken  at  the 
instant  of  highest  head  acceleration. 

The  extent  of  stability  of  the  head  and  mouth  mount  is  illustrated  by  review  of 
the  headspan  plois  (HSP),  plots  612  through  632.  The  change  in  dimension  throughout 
the  runs  is  minimal  • 

Photographic  Data  System 

Pin-registered,  sled-mounted  16  mm  movie  cameras  were  selected  to  avoid  the 
parallax  problem,  as  suggested  by  Shulman's  work.  (2)  To  prevent  loss  of  high 
frequency  response  data,  500  frames/sec  were  exposed.  One  camera  was  mounted 
on  the  sled  at  shoulder  level  and  to  the  right  of  tne  subject,  with  the  film  plane 
paralleling  the  plane  of  expected  motion.  A  second  camera  was  located  directly 
to  the  rear  of  the  subject  at  shoulder  level .  Each  camera  timed  the  center  of  shutter 
opening  of  each  frame  of  film  to  within  -±0.1  msec.  The  timers  produced  positive 
time-locking  of  the  cameras  and  the  analogue  data  on  the  tape  recorder. 

Photographic  lighting  during  acceleration  was  provided  by  sled-mounted  flood 
lights. 

Figure  DP-9  is  an  example  of  the  photographic  information  taken  at  the  point  of 
maximum  differential  head  displacement  during  a  sled  run. 

Data  Train  System 

The  electronic  data  train  used  to  record  and  compute  the  outputs  of  the  transducers 
is  diagrammed  in  Figure  2,  The  train  was  designed  to  minimize  the  manhours  required 
to  obtain  the  desired  measurements  and  to  improve  their  accuracy.  Although  high¬ 
speed  photography  is  simple  and  useful  for  data  collection,  analysis  of  the  films 
requires  an  inordinate  amount  of  time  and  does  not  lend  itself  to  the  precise  deter¬ 
mination  of  instantaneous  acceleration. 

The  redundant  transducer  system  was  chosen  to  permit  the  instantaneous  assessment 
of  the  differential  relationships  between  the  head  and  T]  of  the  subject  and  permitted 
tracking  the  linea  and  angular  components  of  the  head  and  Tj  of  tne  subject  in 
inertial  space. 
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Schematic  Representation  of  Data  Acquisition  and  Processing  System  for 

Transducer  Data 

FIGURE  2 

Signals  leaving  the  transducers  were  amplified  and  scaled  by  integrated  circuit 
operational  amplifiers.  These  units  were  designed  to  permit  the  experimenter  to 
expand  the  range  of  the  transducer  outputs  so  that  the  expected  range  of  experimental 
output  corresponded  roughly  to  full  range  on  the  tape  recorder.  The  setting  was  done 
witn  switchable  resistor  networks  and  the  stability  of  the  amplifiers  was  such  that 
monthly  calibration  was  sufficient  to  guarantee  a  maximum  error  of  0.1%. 

The  greatest  source  of  error  in  the  system  was  contributed  by  the  analog  magnetic 
tape  recorder.  An  Ampex  FR  1800L  tape  recorder  having  an  overall  signal -to-noise 
ratio  of  53  dB,  or  about  0.5%  root  mean  square  noise,  was  selected.  The  remainder 
of  the  data  train  was  chosen  to  reduce  all  other  noise  and  inaccuracy  to  an  insigni¬ 
ficant  level.  The  recorder  was  calibrated  at  the  time  of  each  experimental  run  in 
order  to  achieve  the  full  benefit  of  the  recorder  specifications. 

A  phase-locked  tape  speed  servo  was  used  in  the  analogue  tape  recorder  to  pro¬ 
vide  short-term  constancy  and  linearity,  but  the  electronic  stability  ,vas  such  that 
frequent  recalibration  was  necessary.  In  addition,  the  Ampex  FR  1800L  used  to 
reproduce  the  analogue  signals  for  plotting  and  A/D  conversion  was  also  used  by 
another  experimenter  in  the  laboratories,  which  caused  frequent  calibration  changes. 
In  order  to  minimize  the  time  and  effort  required  to  retain  system  accuracy,  a 
calibration  sequencer  was  built  into  the  system  immediately  prior  to  each  experimental 
run.  This  sequencer  recorded  a  series  of  constant  voltage  levels  on  each  channel  of 
the  tape  recorder.  During  subsequent  data  reduction  (see  Data  Processing  Chapter) 
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the  output  voltage  levels  of  the  calibration  sequence  were  used  to  generate  a  first- 
order  transformation  to  calibrate  the  incoming  data. 

The  data  transmitted  to  a  remote  site,  played  back  through  an  identical 
recorder,  were  sampled  at  2000  samples  per  second  with  a  10— bit  accuracy,  and 
transmitted  to  the  core  of  the  computer  for  processing.  The  sampling  rate  was 
selected  to  be  at  least  10  times  the  expected  maximum  frequency  of  the  response. 

The  mounts  provided  eight  data  channels  which,  along  with  a  chair  accelerometer 
channel,  were  digitized.  Once  the  data  were  resident  in  the  computer,  a  series  of 
processing  programs  was  applied  to  them  as  described  in  the  Data  FVocessing  Chapter. 

Computer  Systems 

The  UNIVAC  1108-3G  multiprocessor  system  installation  at  the  NASA  Slidell 
Computer  Complex  was  utilized  in  the  processing  of  data  and  calculation  of  derived 
variables.  This  system  is  a  general  purpose,  high  performance  processor  consisting 
of  six  tapes  of  subsystems: 

1 .  Central  Processor  Units  (CPU) 

2.  Input/Output  Controller  (IOC) 

3.  Main  Storage  Modules 

4.  Auxiliary  Storage  Subsystem 

5.  System  Interconnection  Components 

6.  Peripheral  Subsystems 

The  UNIVAC  1108  multiprocessor  system  includes  three  central  processors, 
each  of  which  can  perform  all  functions  required  for  the  execution  of  instructions, 
arithmetic,  logic,  input/output,  and  executive  control.  Each  central  processor 
unit  contains  control  registers,  arithmetic  section,  control  section,  input/output 
section,  indexing  section,  and  storage  class  control  section. 

This  system  has  two  input/output  controllers.  The  IOC  is  an  independent  device 
that  controls  the  operation  of  up  to  sixteen  peripheral  subsystems  under  the  direction 
of  the  three  central  processing  units.  Once  an  I/O  request  has  been  issued  to  an 
IOC,  the  IOC  controls  the  operation  and  transfer  of  data  between  main  storage  and 
eight  high-speed  data  channels.  Each  channel  can  transfer  data  at  speeds  up  to 
250,000  words  per  second. 

The  main  storage  module  is  a  magnetic  cor;?  consisting  of  32,768  thirty-six  bit 
words.  This  storage  area  is  expandable  in  pairs  to  a  maximum  of  slightly  more  than 
262  K  words.  The  effective  instruction  cycle  time  for  main  storage  is  750  nanoseconds. 

The  auxiliary  storage  subsystems  consist  of  three  types  of  magnetic  drums:  three 
high-speed  FH-432  drums  with  a  total  of  786,432  words  capacity,four  medium-speed 
FH-1782  drums  with  a  total  of  8,388,608  words  capacity.  Each  word  of  storage 
consists  of  36  bits.  In  the  worst  case  (Fastrand  II  arums)  the  average  access  time  is 
92  milliseconds  with  a  word  transfer  rate  of  25,590  words  per  second. 
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Peripheral  subsystems  include  twenty-two  Uniservo  VIII  -  C  seven  track  tape 
units  having  a  maximum  transfer  rate  of  96,000  six-bit  characters  per  second, 
two  card  punches,  two  card  readers,  and  five  high  speed  printers  with  a  printing 
rate  of  1200  lines  per  minute.  Also  included  are  five  UNIVAC  1004  remote  batch 
terminals,  and  a  number  of  UNIVAC  Uniscope  300  CRT  terminals,  UNIVAC  DCT 
500  data  terminals,  and  IBM  1050  data  terminals. 

A  UNIVAC  500  DCT  demand  terminal  is  located  at  the  NAMRL  Detachment  of 
the  NASA  Michoud  Assembly  Facility,  The  DCT  500  efficiently  transfers  data  between 
the  terminal  and  the  centrally  located  UNIVAC  1108  processor.  The  system  consists 
of  a  132  position  per  line  printer  and  an  American  Standard  Code  Character  set.  It 
can  operate  at  speeds  of  up  to  300  bits  per  second  and  selectable  printer  speeds  of  ten, 
fifteen  or  thirty  characters  per  second. 

The  plotting  of  digital  information  as  presented  in  this  report  is  accomplished  by 
the  use  of  a  Stromberg-Carlson  Model  S-C  4020  microfilm  printer/plotter.  Input 
to  the  S-C  4020  is  in  the  form  of  1/2  inch  magnetic  tape  end  is  an  off-line  process 
utilizing  tape  created  by  the  UNIVAC  1 108  processor,  A  set  of  special  software 
routines  developed  and  maintained  by  the  computer  installation  is  used  in  the  genera¬ 
tion  of  the  tapes  for  plotting. 

The  S-C  4020  utilized  IBM  seven  track  tapes  in  densities  of  200,  556,  or  800 
bits  per  inch.  It  will  accept  data  from  magnetic  tape  at  input  rates  up  to  62,500 
six-bit  characters  per  second  and  print  these  data  at  speeds  in  excess  of  17,000 
alpha-numeric  or  symbolic  characters  per  second.  The  plotting  of  scientific  curves 
can  be  accomplished  at  a  speed  of  12,500  plotting  points  per  second,  or  if  being 
used  exclusively  as  a  printer,  will  print  7,000  lines  per  minute.  Graphs  or  frames 
combining  characters,  vectors,  and  curves  vary  with  complexity  but  an  average 
annotated  graph  can  be  recorded  in  milliseconds. 

The  S-C  4020  printer/plotter  is  capable  of  translating  high-speed  digital  data 
into  alpha-numeric  characters,  symbols  ana  lines.  This  is  accomplished  by  means  of 

the  CHARACTRON^  Shaped  Beam  Tube,  This  tube  contains  a  built-in  character¬ 
forming  section  consisting  of  an  electron  gun,  character-selection  plates  and  a 
character  Matrix.  The  electron  gun  directs  an  electron  beam  towards  the  Matrix 
through  the  character-selection  plates.  The  beam  is  deflected  by  these  plates  so  that 
it  passes  through  the  desired  Matrix  aperature.  As  the  beam  leaves  the  Matrix,  it 
has  been  shaped  in  the  form  of  the  character  aperature  through  which  it  passed.  The 
beam  is  then  deflected  to  the  proper  position  on  the  tube  face.  Images  on  the  tube 
face  are  recorded  by  either  a  35  mm  or  a  16  mm  camera  positioned  above  the  tube. 

The  S-C  4020  displays  characters  on  an  area  of  1024  horizontal  and  1024  vertical 
spaces  ("rasters")  resulting  in  over  one  million  addressable  positions.  After  the  desired 
margins  have  been  reserved,  the  maximum  range  of  the  plot  is  scaled  to  the  remaining 
rasters  and  the  grid  is  drawn.  Retaining  the  same  scale  factors,  all  data  points  are 
converted  to  its  equivalent  rasters  and  projected  on  the  screen  and  film.  The  center 
of  any  available  character  may  be  projected  on  any  of  the  addressable  positions.  The 
minimum  area  occupied  by  a  character  in  the  grid  is  three  horizontal  by  three  vertical 
spaces . 
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The  following  equations  are  used  to  compute  the  scale  factors.  Variables  A 
and  B  represent  the  X  axis  scale  factors  arJ  C  and  D  represent  the  Y  axis  scale 
factor-. 

(1023,  -  FMR)  -  RML 
A  =  XR-XL 

B  =  FM.  -  A  x  XL 

(1023.  -  FMT)  -  FMB 
C  YT  -  YB 

D  =  FMB  -C  +  YB 

where: 

FMR  =  spaces  reserved  for  right  margin 

FML  =  spaces  reserved  for  left  margin 

XR  =  right  limit  of  the  scaled  plotting  area 

XL  =  left  limit  of  the  scaled  plotting  area 

FMT  =  spaces  reserved  for  top  margin 

FMB  =  spaces  reserved  for  bottom  margin 

YT  =  top  limit  of  the  scaled  plotting  area 

YB  -  bottom  limit  of  the  scaled  plotting  area 

These  equations  are  used  to  convert  data  coordinates  X  and  Y  into  rosier 
coordinates  for  plotting. 

IX  =  AxX  +  B 

IY  =  C  x  Y  +  D 

A  Xerox  Copyflo  11-3  printer  is  available  to  produce  single  copies  continuously 
from  either  16  mm  or  35  mm  roll  microfilm.  These  hard  copies  are  produced  at  the 
rate  of  28  frames  (copies)  per  minute. 

Accelerator  System 

The  WHAM  II,  an  acceleration  device  designed,  constructed  and  operated  by 
Wayne  State  University  Biomechanics  Research  Center  was  modified  for  use  in  this 
investigation  by  adding  the  numerous  safety  redundancies  required  for  man-rating. 
The  device  permitted  use  of  the  acceleration  end  for  producing  the  experimental 
acceleration  pulse  and  hence  the  subject  was  at  rest  during  his  initial  condition  as 
required  by  the  instrumentation  system. 
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A  sled  was  designed  and  constructed  by  WSU  especially  for  this  project.  It  has 
a  steel  seat  which  permits  forward  and  backward  tilting  in  the  mid-sagittal  plane  as 
well  as  rotation  in  the  horizontal  plane.  A  modified  A-l  aircraft  lap  belt  and 
shoulder  harness  with  an  inverted  V  and  a  chest  safety  strap  formed  tne  restraint  system. 
This  was  secured  as  tightly  as  the  subject  would  allow  for  each  run,  with  the  sole 
exception  that  the  chest  safety  strap  was  simply  made  snug,  with  very  little  tension. 

Seat  back  and  shoulder  harness  heights  were  made  adjustable  to  the  individual, 
thus  achieving  maximum  subject  safety  and  restraint  during  acceleration  and 
deceleration.  This  also  permitted  unhampered  photography  of  the  Tj  mount  and  the 
rear  of  the  head  and  neck. 

Experimental  Controls  and  Determinations 

The  following  methods  were  used  to  determine  the  effects  of  different  values  of 
a  single  variable  and  to  control  the  effects  of  other  variables  believed  to  influence 
the  dynamic  response: 

1 .  To  control  the  effect  of  rate  of  onset  of  acceleration  it  was  kept  con¬ 
stant  at  a  nominal  250  G/sec  for  the  initial  series  of  runs  on  each  subject. 

2.  To  determine  the  effect  of  rate  of  onset  of  acceleration  the  first  step  was 
to  repeat  the  entire  initial  profile  of  increments  of  peak  sled  acceleration  at  a  rate 
of  onset  double  that  of  the  initial  series;  that  is,  at  a  nominal  500  G/sec. 

3.  To  control  the  effects  of  individual  variations  in  muscle  mass,  strength 
and  other  anthropometric  and  anthropomorphic  sources  of  variance,  each  subject 
was  used  as  his  own  control. 

4.  To  determine  the  effects  of  the  subject  size,  subjects  were  selected  on 
the  basis  of  their  sitting  height. 

5.  To  control  the  effect  of  variations  in  peak  sled  acceleration  the 
accelerator  control  settings  were  reproduced  as  nearly  as  possible  for  each  run  at  a 
given  peak  acceleration  level . 

6.  To  determine  the  effect  of  variations  in  peak  sled  acceleration  each 
subject  was  exposed  to  individual  runs  at  1-G  increments  from  3  G  through  10  G. 

7.  To  control  the  effect  of  sled  acceleration  vector  direction  it  was  kept 
constant  for  all  runs. 

8.  To  control  the  effect  of  variations  in  distance  between  accelerometers 
in  the  mount  and  that  of  spatial  relationships  between  accelerometers  and  rate 
gyroscope,  the  accelerometers  and  rate  gyroscope  were  locked  into  place  on  the 
transducer  module,  and  the  entire  transducer  module  was  shifted  to  the  corresponding 
mating  and  anatomical  module  of  subsequent  subjects.  Thus  the  spatial  relationships 
of  the  accelerometers  and  rate  gyroscopes  relative  to  each  other  for  a  given  mount 
were  constant  for  the  entire  series  of  runs. 
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9.  To  control  the  effect  of  differing  distances  between  subjects  from 
anatomical  landmarks  on  any  individual  to  the  transducer  module  ana  photographic 
targets,  precise  photographic  and  x-ray  anthropometry  was  performed  on  all  subjects 
with  this  anatomical  and  transducer  modules  in  place.  Distances  from  external  and 
internal  landmarks  to  transducers  arxl  to  the  photographic  targets  were  precisely 
known  and  did  not  vary  from  run  to  run  on  a  given  subject.  Therefore,  transformation 
of  the  measured  accelerations  and  velocity  from  the  point  of  measurement  to  the 
center  of  mass  of  the  head  of  a  given  subject  or  to  any  specified  location  on  the  head 
of  that  subject  was  possible  within  a  sagittal  plane. 

10.  To  determine  the  effect  of  instrumentation  weight  and  its  distribution 
each  subject  was  exposed  to  several  runs  at  both  rates  of  onset  with  all  mounts  and 
photogruphic  targets  but  with  no  instrumentation  mounted  except  in  the  Tj  mount. 

These  runs  were  termed  unencumbered  runs,  and  photographic  data  derived  therefrom 
will  be  compared  to  those  derived  from  fully  instrumented  runs.  These  data  have  not 
yet  been  analyzed. 

11.  To  determine  run-to-run  repeatability  for  a  given  set  of  run  parameters 
on  individual  subjects,  multiple  runs  were  performed  on  selected  subjects  at  selected 
sled  peak  acceleration  and  rate  of  onset  of  acceleration  values.  These  data  are 
presented  in  plots  517  through  611 . 

12.  To  determine  the  charge  in  distance  between  the  Tj  and  head  anatomical 
coordinate  systems  due  to  the  dynamic  event,  the  change  which  occurred  during  a 
normal  head  nod  of  the  fully  instrumented  and  restrained  volunteer  were  measured 
photographically.  The  reduced  data  from  the  special  version  of  TROUTP  (known  as 
TROUTH)  were  compared  with  those  from  the  dynamic  event  and  the  true  cnange  due 
to  the  acceleration  determined. 

The  effect  of  the  restraint  variable  is  assumed  to  be  independent  of  the  funda¬ 
mental  relationship  between  the  Ti  acceleration  and  the  head  acceleration.  However, 
as  a  practical  matter,  each  individual  restraint  system  was  adjusted  as  tightly  as 
possible  for  each  run,  both  as  a  safety  measure  and  to  keep  the  photographic  targets 
within  the  field  of  view  of  the  sled-mounted  cameras.  Parenthetically  it  should  be 
noted  that  the  restraint  adjustment  procedure  was  determined  experimentally  to  result 
in  apparent  excellent  control  of  the  restraint  variable,  as  will  be  shown  later.  An 
analysis  of  the  precise  effects  of  variations  in  the  restraint  of  an  individual  subject 
on  the  dynamic  response  of  the  head  and  neck  must  be  performed,  but  has  not  yet 
been  accompl  ished , 

Subjects  and  Procedure 

Volunteer  subjects  from  the  U.  S.  Army  were  initially  screened  carefully  by 
medical  history,  examination  of  personnel  and  health  records,  and  independent 
interviews  with  two  flight  surgeons. 

Those  prospective  subjects  resulting  from  the  screening  were  measured  anthropo- 
metrically  and  preliminary  selection  was  based  on  the  single  criterion  of  sitting 
height. 
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Candidates  surviving  the  second  selection  were  then  examined  thoroughly  by 
specialists  in  aerospace  medicine,  dentistry,  orthodontics,  radiology,  otorhinolaryn- 
gology,ophthalmology,  neurology  and  psychiatry  and  vestibular  physiology.  All 
car  didates  who  qualified  in  these  examinations  became  subjects  in  the  study;  only 
abut  5%  of  the  volunteers  successfully  completed  the  entire  selection  program. 

Two  flight  surgeons  of  the  research  team,  who  were  physically  qualified  by  the 
same  standards,  volunteered  for  duty  and  made  several  experimental  runs  as  subjects. 

A  careful  interim  medical  history  was  obtained  from  each  subject  prior  to  a  run, 
and  a  physical  examination  and  urinalysis  were  performed. 

Following  each  run,  the  subject  was  carefully  examined  and  studies  made  of  his 
vestibular  function,  cardiac  and  neurological  status,  and  urine.  A  24-hour  post-run 
history  was  taken  prior  to  his  next  run.  One  run  per  subject  per  day  was  the 
maximum  exposure  rate. 

The  seated  individuals  were  subjected  to  a  -Gx  impact  acceleration  at  the 
accelerction  end.  The  end  of  the  experimental  conditions  of  interest  occurred  when 
dynamic  response  was  substantially  complete,  which  occurred  in  every  instance  prior 
to  the  end  of  the  initial  acceleration  pulse.  The  sled  achieved  a  peak  velocity 
which  remained  relatively  constant  until  the  sled  brakes  were  activated.  Sled  breaking 
produced  a  smooth  2-G  deceleration  until  zero  sled  velocity  was  attained. 

Acceleration  pulse  shape  for  each  run  was  triangular  with  a  long  decay. 

A  complete  experimental  run  profile  of  nominal  values  for  each  subject  included: 

1.  3G  through  10G  at  250  G/sec  in  1-G  increments. 

2.  3G  through  8G  at  500  G/sec  in  1-G  increments. 

3.  6G,  250  G/sec;  10G,  250  G/sec;  and  6  G,  500  G/sec  unencumbered 

runs. 

This  involved  a  design  goal  of  17  runs  for  each  subject. 

Of  the  236  human  exposures  to  -Gx  impact  acceleration  in  this  program,  30  runs 
were  incompletely  instrumented  and  7  were  conducted  on  subjects  wearing  helmets; 
therefore,  199  human  runs  fulfilled  the  requirements  for  data  analysis. 

Six  subjects  completed  the  stated  17-run  profile.  An  additional  two  subjects 
completed  runs  through  6G,  500  G/sec,  omitting  only  the  7G  and  8G,  500  G/sec 
runs,  for  a  total  of  15  runs  per  subject.  Of  these  eight  subjects,  two  also  completed 
9G,  500  G/sec  runs.  Additional  runs  were  made,  ranging  in  number  from  1  through 
1 1  by  various  other  subjects. 

Due  to  the  extensive  instrumentation,  data  retrieval,  and  data  processing 
arrengements  made  for  each  experimental  exposure,  large  masses  of  data  are  available. 
It  was  necessary  to  select  for  tne  present  report  data  that  were  available  in  nine 
categories  established  by  peak  sled  acceleration  and  rate  of  onset  from  12  subjects. 

The  selection  was  not  a  random,  but  was  made  on  the  basis  of  data  available  to  meet 
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the  requirements  of  this  report.  A  listing  of  the  runs  selected  is  presented  in 
the  Results  Chapter. 

The  effort  has  resulted  in  the  formulation  of  an  extensive  body  of  experimental 
data.  This  report  is  a  limited  selection  from  the  data  base.  Data  presented  consists 
of  kinematic,  and  anthropometric  variables  requested  by  the  Department  of 
Transportation  along  with  a  detailed  discussion  of  the  methods  of  acquiring,  processing, 
and  analyzing  the  variables.  Validation  procedures  and  variables  specifically 
selected  to  demonstrate  validation  are  also  presented. 
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THEORETICAL  KINEMATICS  OF  THE  EXPERIMENTS  CHAPTER 

Daniel  J.  Thomas 


GENERAL 

Within  the  context  of  the  experimental  program  the  effort  to  measure  the  exact 
kinematics  of  the  human  head  and  neck  is  of  paramount  importance.  The  results 
of  the  effort  constitutes  the  major  portion  of  this  report.  A  detailed  description  of 
the  theoretical  kinematics  to  understand  and  describe  the  measurements,  derived 
variables,  and  validation  procedures  is  required. 

The  study  of  the  mechanics  of  human  subjects  can  be  organized  according  to 
classical  mechanics.  The  size  and  shape  of  tne  human  beings  constitute  geometry. 
The  mass  distribution  parameters  associated  with  the  size  and  shape  constitute  the 
statics  of  human  subjects.  The  motions  of  anatomical  segments  constitute  the  kine¬ 
matics.  The  forces  and  torques  underlying  these  motions  constitute  the  dynamics. 

The  classification  of  the  mechanics  of  human  subjects  into  geometry,  statics, 
kinematics,  and  dynamics  is  completely  sufficient  provided  the  analysis  is  restricted 
to  the  consideration  of  rigid  body  mechanics.  For  these  experiments  the  rigid  bodies 
are: 


1.  The  head  with  attached  mouth  and  head  mounts,  and  their  attached 
targets  and  transducers. 

2.  The  first  thoracic  vertebra  with  attached  spine  mount,  transducers 
and  targets.  The  kinematic  measurements  of  the  experiment  documents  the  motions 
between  these  two  rigid  bodies  and  is  valid  to  the  degree  that  each  body  remains 
rigid.  This  is  what  is  meant  by  the  measurement  of  the  response  of  the  head  and 
neck.  In  order  to  systematically  provide  for  these  measurements  and  to  make  the 
coordinate  system  transformations  necessary  to  analyze  the  data  a  standard  procedure 
for  defining  coordinate  systems  and  kinematic  variables  and  a  minimal  set  of  defined 
coordinate  systems  is  required, 

DEFINITIONS  OF  COORDINATE  SYSTEMS 

Although  the  experiments  to  date  were  confined  to  two  dimensions,  the  approach 
to  the  description  of  the  coordinate  systems  end  variables  is  applicable  in  three 
dimensions.  This  was  achieved  by  first  defining  a  set  of  basic  orthogonal  three 
dimensional  coordinate  systems.  The  variables  are  all  expressed  in  terms  of  these 
coordinate  systems.  Although  a  standard  procedure  for  defining  coordinate  systems 
has  been  adopted  by  the  laboratory,  this  system  was  not  in  effect  at  the  time  of 
collection  and  initial  processing  of  the  two  dimensional  data  from  the  human  runs 
conducted  at  Wayne  State  University  (1).  Therefore,  the  current  program  and 
variables  which  were  measured  or  used  in  intermediate  calculations  do  not  conform 
to  the  laboratory  standard.  However,  all  variables  which  are  plotted  for  this  report 
have  been  expressed  with'n  the  laboratory  standard  as  a  final  step  prior  to  plotting. 


16 


The  minimum  rules  for  coordinate  system  definition  which  were  identified  and 
selected  are:  . 

1.  Select  an  origin. 

2.  Select  the  first  and  second  axes  of  the  system,  which  must  be  at  right 
angles  to  each  other. 

3.  Define  the  third  axis  and  all  operations,  parameters  and  variables 
within  the  coordinate  system  by  use  of  a  right-hand  rule. 

4.  +0  degrees  to  +90  degrees  (first  quadrant)  will  always  be  subtended  by 
the  positive  axes  in  the  plane  in  which  the  angle  Is  described. 

For  this  report  the  notation  X,  Y,  Z  is  used  to  specify  the  1st,  2nd,  and  3rd 
axis  of  any  defined  coordinate  system. 

The  units  of  measure  for  all  kinematic  variables  will  be  centimeters,  radians, 
and  milliseconds,  with  the  exception  that  categorization  and  identification  of  runs 
uses  G  for  acceleration  and  G  per  second  for  rate  of  onset.  One  G  is  equivalent 
to  .000930616  centimeters  per  millisecond  squared.  Also,  values  extracted  from 
the  plots  for  discussion  in  thi  text  have  been  converted  into  more  readily  compre¬ 
hended  units,  in  which  case  the  units  are  given. 

The  coordinate  systems  required  for  these  experiments  are  grouped  into  four 
categories: 

1 .  Reference  -  is  any  coordinate  system  in  which  measuremehts  are 
directly  taken, 

2.  Instrumentation  -  is  any  coordinate  system  established  by  instrumen¬ 
tation  or  photographic  target  location. 

3.  Anatomical  -  is  any  coordinate  system  established  by  the  anatomical 
landmarks  of  a  test  subject. 

4.  Principal  -  is  any  coordinate  system  established  by  determination  of 
the  center  of  gravity  and  principal  axes  of  the  moment  of  inertia  of  a  defined  rigid 
segment. 

The  defined  coordinate  systems  remain  rigid  throughout  the  experiment. 

The  individual  coordinate  systems  required  for  these  experiments  are  further 
defined  below. 

The  reference  coordinate  systems  (figure  1)  are: 

1.  The  inertial  reference  coordinate  system  is  one  in  which  the  linear 
and  angular  ccceleration  is  zero.  The  origin  at  data  analysis  time  zero  (DATZ)  is 
at  zero  velocity  at  the  origin  of  the  laboratory  reference  system  defined  below. 

The  +Z  axis  is  opposite  in  direction  to  gravity  and  the  +X  axis  is  opposite  in  direction 
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FIGURE  1 


to  the  acceleration  stroke,  This  coordinate  system  translates  relative  to  the  laboratory 
reference  system  at  a  velocity  dependent  on  the  acceleration  of  gravity. 

2.  The  sled  reference  coordinate  system  has  its  origin  fixed  to  the  sled 
chair  target  center  with  axes  parallel  to  the  inertial  reference  system. 

3.  The  laboratory  reference  coordinate  system  has  its  origin  fixed  at  the 
point  occupied  by  the  sled  chair  target  center  prior  to  the  onset  of  acceleration 
with  axes  parallel  to  the  inertial  sled  reference  system. 

4.  The  Telereadex  reference  system  is  established  by  the  film  reader  device 
at  the  sled  chair  origin  with  axes  slightly  rotated  from  the  slea  reference  axes.  This 
coordinate  system  is  not  illustrated  because  of  its  minor  differences  from  the  sled 
reference  coordinate  system. 

The  instrumentation  coordinate  systems  are: 

1.  The  photo  head  mount  instrumentation  coordinate  system  (Figure  2)  has 
its  origin  at  the  midpoint  between  the  head  mount  target  and  the  mouth  mount  target. 
The  +X  axis  is  in  the  direction  of  the  line  between  the  targets  from  the  head  mount 
target  to  the  mouth  mount  and  the  +Z  axis  is  perpendicular  to  the  line  between  the 
targets  in  the  superior  direction  of  the  man. 

2.  The  photo  spine  mount  instrumentation  coordinate  system  (Figure  3) 
has  its  origin  at  the  midpoint  between  the  spine  mount  targets.  The  +X  axis  is  in 
the  direction  of  the  line  connecting  the  targets  in  the  anterior  direction  of  the  man 
and  the  +Z  axis  is  perpendicular  to  this  line  in  the  superior  direction  of  the  man. 

3.  The  transducer  mouth  mount  instrumentation  coordinate  system  (Figure  2) 
has  its  origin  at  the  midpoint  between  the  mouth  accelerometers.  The  +X  axis  is 
parallel  to  the  superior  surface  of  the  mouth  mount  in  the  anterior  direction  and  the 
+Z  axis  is  perpendicular  to  the  +X  axis  in  the  superior  direction. 

4.  The  transducer  head  mount  instrumentation  coordinate  system  (Figure  2) 
has  its  origin  at  the  midpoint  between  the  top  of  the  head  accelerometers.  The  +X 
axis  is  parallel  to  the  superior  surface  of  the  mount  in  the  anterior  direction  and  the 
+Z  axis,  is  perpendicular  to  the  +X  axis  in  the  superior  direction. 

5.  The  transducer  center  of  head  instrumentation  coordinate  system  is 
identical  to  the  photo  head  mount  instrumentation  coordinate  system  by  virtue  of  the 
fact  that  the  mouth  mount  target  was  placed  at  the  midpoint  between  the  mouth 
mount  accelerometers  and  the  head  mount  target  was  placed  at  the  midpoint  between 
the  head  mount  accelerometers. 

6.  The  transducer  spine  mount  coordinate  system  (Figure  3)  has  its  origin 
at  the  midpoint  between  the  spine  mount  accelerometers.  The  +X  axis  is  parallel  to 
the  superior  surface  of  the  mount  in  the  anterior  direction  and  the  +Z  axis  is  perpen¬ 
dicular  to  the  +X  axis  in  the  superior  direction.  This  coordinate  system  has  its  origin 
identical  to  the  photo  spine  mount  instrumentation  coordinate  system  with  a  fixed 
angle  (0.2  radians)  between  respective  +X  axes. 
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INSTRUMENTATION  AND  ANATOMICAL 
COORDINATE  SYSTEMS  (HEAD) 


INSTRUMENTATION  COORDINATE  SYSTEMS 
P= PHOTO  HEAD  MOUNT 
A=  TRANSDUCER  CENTER  OF  HEAD 
M“ TRANSDUCER  MOUTH  MOUNT 
D  =  TRANSDUCER  HEAD  MOUNT 


ANATOMICAL  COORDINATE  SYSTEM 

H  =  HEAD  ANATOMICAL 


FIGURE  2 
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INSTRUMENTATION  AND  ANATOMICAL 
COORDINATE  SYSTEMS  (SPINE) 


V  =  FIRST  THORACIC  VERTEBRA 

INSTRUMENTATION  COORDINATE  SYSTEMS 
K  =  PHOTO  SPINE  MOUNT 
T  =  TRANSDUCER  SPINE  MOUNT 

ANATOMICAL  COORDINATE  SYSTEM 

T,  *  SPINE  (T.) 


FIGURE  3 


The  anatomical  coordinate  systems  are: 

1 ,  The  head  anatomical  coordinate  system  (Figure  2)  is  derived  from  an 
anatomical  plane  which  is  specified  by  the  superior  edge  of  each  auditory  meatus 
and  by  the  infraorbital  notcnes.  The  origin  is  at  the  midpoint  of  a  line  connecting 
the  superior  edges  of  the  right  and  left  auditory  meati.  The  +X  axis  is  from  the  origin 
through  the  midpoint  of  a  line  connecting  the  infraorbital  notches  in  the  anatomical 
plane.  The  +Z  axis  is  from  the  origin  in  the  superior  direction  perpendicular  to  the 
anatomical  plane.  For  the  purpose  of  this  work,  the  XZ  plane  thus  described  .is 
considered  a  precise  definition  of  the  mid-sagittal  plane. 

2.  The  spine  (Tl)  anatomical  coordinate  system  (Figures  3  and  4)  has  its 
origin  at  the  anterior  superior  corner  of  the  first  thoracic  vertebral  body  (T]).  The 
+X  axis  is  defined  by  a  direction  connecting  the  midpoint  of  the  superior  corner  and 
the  inferior  corner  of  the  posterior  spinous  process  of  Tl  to  the  anterior  superior 
corner  of  Tl .  The  +Z  axis  is  perpendicular  to  the  X  axis  in  the  plane  of  the  above 
three  points  and  in  the  superior  direction. 

All  principal  coordinate  systems  must  be  experimentally  determined  and  therefore 
are  expressed  in  terms  of  the  cnosen  anatomical  coordinate  system.  The  principal 
coordinate  system  is: 

1.  The  head  principal  coordinate  system  is  expressed  in  terms  of  the  head 
anatomical  coordinate  system  and  has  its  own  origin  at  the  center  of  gravity  of  the 
head  as  defined  in  the  Tulane  and  NAMRL  studies,  both  of  which  are  still  unpub¬ 
lished  (  ).  There  is  insufficient  information  to  specify  the  principal  axes  of  the 
moment  of  inertia  so  the  axes  remains  undefined.  The  origin  (C.G.)  is  illustrated 
in  Figure  5. 


ANATOMICAL  COORDINATE 
SYSTEM  (SPINE) 


T,  =  SPINE  (Ti)  ANATOMICAL  COORDINATE  SYSTEMS 
V  =  FIRST  THORACIC  VERTEBRA 


DEFINITIONS  OF  PHYSICAL  TERMS  USED  TO  DEFINE  VARIABLES 

The  kinematics  of  human  subjects  constitute  the  major  portion  of  the  data 
presented  in  this  report.  The  anatomical  segments  of  interest  are  treated  as  rigid 
bodies.  Therefore,  the  terms  used  to  define  the  variables  are  taken  from  rigid 
body  kinematics.  The  definition  of  physical  terms  is  as  follows: 

1 .  Displacement  is  a  vector  quantity  describing  the  position  of  one  point 
relative  to  another  point. 

2.  Angular  displacement  is  a  transformation  describing  the  orientation  of 
one  coordinate  system  relative  to  another  coordinate  system.  For  two  dimensional 
motion  where  the  angular  displacement  is  restricted  to  a  single  plane,  the  angular 
displacement  is  identical  to  an  angle  defined  by  the  angular  position  of  one  vector 
direction  relative  to  another  vector  direction. 

3.  Velocity  is  a  vector  quantity  describing  the  time  rate  of  change  of 
displacement. 

4.  Angular  velocity  is  a  vector  quantity  describing  the  instantaneous 
angular  time  rate  of  change  of  the  angular  displacement  of  one  coordinate  system 
relative  to  another  coordinate  system.  For  two  dimensional  motion  in  which  angular 
displacement  is  an  angle,  angular  velocity  is  the  time  rate  of  change  of  the  angle. 

5.  Acceleration  is  a  vector  quantity  describing  the  time  rate  of  change  of 

velocity. 

6.  Angular  acceleration  is  a  vector  quantity  describing  the  time  rate  of 
change  of  angular  velocity. 

For  this  report  all  points  selected  for  displacement  measurements  are  origins  of 
defined  coordinate  systems.  Since  all  motions  discussed  in  this  report  are  two 
dimensional  and  constrained  to  the  mid-sagittal  plane  of  the  subject,  description  of 
vector  directions  are  sufficient  to  define  the  angle  identical  to  angular  displacement. 

All  vector  directions  selected  for  angular  displacement  measurements  are  the 
directions  of  the  axes  of  defined  coordinate  systems.  In  general,  all  terms  are  vector 
quantities  except  angular  displacement.  In  particular  for  this  report  since  the  motions 
are  two  dimensional,  the  linear  quantities  are  vectors  described  by  two  components 
and  the  angular  quantities  can  be  treated  as  scalers  requiring  only  a  magnitude.  In 
most  cases  the  magnitude  of  the  linear  quantities  is  also  given. 

DEFINITION  OF  VARIABLES 

The  kinematic  variables  presented  in  graphical  form  are  defined.  Selected 
variables  that  are  needed  to  clarify  the  definitions  of  variables  which  are  graphi¬ 
cally  presented  are  also  defined.  The  definitions  consist  of  vector  and  scaler 
magnitudes,  vector  Cartesian  components  and  one  polar  angle  and  are  arranged  by 
computer  program  names  or  stage  of  processing. 
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A  variable  of  the  same  name  is  fundamentally  the  same  physical  quantity  with 
necessary  processing  changes  introduced  by  each  program,  the  definitions  will 
reflect  this.  The  groupings  of  the  variables  in  order  are: 

PLOT 

INTERMEDIATE 
TROUTP 
PHD  AT  A 

REDUCED  PHOTOGRAPHIC  DATA 

TROUT 

OUTVAR 

D  SCALE 

BETH02 

BETHBH 

ANALOGUE 

ICONDS 

ANXRAY 


The  definition  of  the  variables  includes  a  statement  of  its  physical  nature  along 
with  a  statement  of  source  where  considered  necessary  for  clarity.  A  variable  once 
defined  under  a  grouping  will  only  have  its  changes  noted  within  succeeding  groupings. 

The  PLOT  grouping  consists  only  of  variables  appearing  on  the  plots.  The 
INTERMEDIATE  grouping  consists  of  two  simply  computed  variables  that  ere  required 
to  compute  the  plotted  variable.  The  other  groupings  refer  to  the  analogue  data, 
programs  or  photographic  data  reduction  which  is  explained  in  the  appropriate  sections. 
Selected  important  angular  and  linear  displacement  variables  are  illustrated  by 
drawings.  The  source  of  all  variables  which  are  presented  is  illustrated  in  flow  chart 
format.  The  charts  showing  the  source  of  the  plotted  variables  do  not  show  P|.OT 
blocks.  These  blocks  are  implied. 

Two  important  constants  fixed  by  the  transducer  mount  design  are  given  after 
the  variable  definitions. 


PLOT  Variables 
DATZ 


HSP 

MCT 


Definitions 

is  Data  Analysis  Time  Zero  which  is  36  milliseconds 
prior  to  the  statistically  determined  sled  peak 
acceleration. 

is  the  magnitude  of  distance  of  the  head  target 
center  from  the  mouth  target  center,  Figures  6 
and  8. 

is  the  magnitude  of  displacement  of  the  head 
principal  coordinate  system  origin  (C.G.)  relative 
to  the  spine  (Tl)  anatomical  coordinate  system 
origin,  as  computed  from  PLOT  XCN  and  PLOT 
ZCN,  Figures  7  and  8. 
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TARGET  (3) 


MOUTH 
TARGET  14) 


ILLUSTRATION  OF  HSP  AND  SUBJECT 
MOUNTED  PHOTOGRAPHIC  TARGETS 


FIGURE  6 


SOURCE  OF  VARIABLES  XCN,  ZCN,  MCI,  MCX,  MCZ,  HSR 


ILLUSTRATION  OF  VARIABLE  MCT,  COMPONENTS 
MCX.MCZ  IN  THE  T,  ANATOMICAL  COORDINATE 
SYSTEM  DIRECTIONS 


FIGURE  9 


FIGURE  12 


\  Ti  ANATOMICAL 


FIGURE  17 


MCX  is  the  component  of  displacement  of  the  hoad  (C.G  J 

relative  to  the  (T \)  anatomical  coordinate  system 
origin  along  the  +X  axis  of  the  T ]  anatomical 
coordinate  system,  Figures  8  and  9. 

MCZ  is  the  component  of  displacement  of  the  head  (C.G  j 

relative  to  the  Tj  anatomical  coordinate  system 
origin  along  the  +Z  axis  of  the  Tj  anatomical 
coordinate  system,  Figures  8  and  9. 

NNV  is  the  magnitude  of  the  velocity  of  the  head(C.G.) 

relative  to  the Tj  anatomical  coordinate  system 
origin.  Figure  10. 

is  the  magnitude  of  the  angular  displacement  of  the 
+X  axis  of  the  T  i  anatomical  coordinate  system 
relative  to  the  +X  axis  of  the  laboratory  coordinate 
system.  Figures  11  and  12. 

is  the  magnitude  of  the  angular  displacement  of  the 
+X  axis  of  the  head  anatomical  coordinate  system 
relative  to  the  +X  axis  of  the  laboratory  coordinate 
system  as  computed  from  photographic  data. 

Figures  11  and  13. 

is  the  magnitude  of  the  angular  displacement  of  the 
+X  axis  of  the  head  anatomical  coordinate  system 
relative  to  the  +X  axis  of  the  Tj  anatomical 
coordinate  system  as  computed  from  photographic 
data.  It  is  PLOT  PDH  minus  PLOT  NRN, 

Figures  12  and  14. 

RAH  is  the  magnitude  of  the  angular  acceleration  of  the 

angle  from  the  +X  axis  of  the  head  anatomical 
coordinate  system  relative  to  the  +X  axis  of  the 
laboratory  reference  coordinate  system  as  computed 
from  the  rate  gyroscope.  It  is  the  time  rate  of 
change  of  PLOT  RVH,  Figure  15. 

RDH  is  the  same  as  PLOT  PDH  except  that  it  is  computed 

from  the  rate  gyroscope  data.  Figure  16. 

RHN  is  the  magnitude  ot  the  angular  velocity  of  the 

angle  from  the  +X  axis  of  the  head  anatomical 
coordinate  system  relative  to  the  +X  axis  of  the  T  j 
anatomical  coordinate  system  as  computed  from 
the  rate  gyroscopes,  Figure  17. 


NRN 


PDH 


PHN 
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is  the  magnitude  of  tho  angular  velocity  of  the 
angle  from  the  +X  axis  of  the  head  anatomical 
coordinate  system  relative  to  the  +X  axis  of  the 
laboratory  reference  coordinate  system  as  computed 
from  the  rate  gyroscope ,  Figure  1 5 .  “  •  . 

is  the  component  of  displacement  of  the  head  C.G. 
relative  to  the  T]  anatomical  coordinate  system 
origin  along  an  axis  parallel  to  the  +X  axis  of  the 
sled  reference  coordinate  system.  Figures  7  and 
8.  It  is  TROUTP  XCL  minus  TROUTP  XNA. 

is  the  component  of  the  velocity  of  the  head  C.G. 
relative  to  the  sled  reference  coordinate  system 
origin  along  an  axis  parallel  to  the  +X  axis  of  the 
sled  reference  coordinate  system.  It  is  the  time 
rate  of  change  of  TROUTP  XCL,  Figure  10. 

is  the  displacement  of  the  sled  coordinate  system 
origin  relative  to  the  laboratory  coordinate  system 
origin.  This  displacement  is  constrained  to  the 
laboratory  -X  axis.  Figure  18. 

is  the  component  of  the  velocity  of  the  head  C.G. 
relative  to  the  T]  anatomical  coordinate  system 
origin  along  an  axis  parallel  to  the  +X  axis  of  the 
sled  reference  coordinate  system.  It  is  XCV  minus 
INTERMEDIATE  XV A,  Figure  10. 

Is  the  velocity  of  the  sled  coordinate  system  origin 
relative  to  the  laboratory  coordinate  system.  This 
velocity  is  constrained  to  the  laboratory  -X  axis. 
Figure  18. 

is  the  magnitude  of  the  velocity  of  the  head  C.G. 
relative  to  the  sled  reference  coordinate  system 
origin  as  computed  from  PLOT  XCV  and  PLOT  ZCV, 
Figure  10. 

is  the  component  of  the  displacement  of  the  head 
C.G.  relative  to  theTj  anatomical  coordinate 
system  origin  along  an  axis  parallel  to  the  +Z  axis 
of  the  sled  reference  coordinate  system  data.  It 
is  TROUTP  ZCL  minus  INTERMEDIATE  ZNA, 
Figures  7  and  8. 

is  the  component  of  velocity  of  the  head  C.G. 
relative  to  the  sled  reference  coordinate  system 
origin  along  an  axis  parallel  to  the  +X  axis  of 
the  sled  reference  coordinate  system.  It  is  the 
time  rate  of  change  of  TROUTP  ZCL,  Figure  10. 


ZNV  is  the  component  of  velocity  of  the  head  C.G. 

relative  to  the  Tj  principal  coordinate  system 
origin  along  an  axis  parallel  to  the  +Z  axis  of  the 
sled  reference  coordinate  system.  It  is  PLOT  ZCV 
minus  INTERMEDIATE  ZVA#  Figure  10. 

4FD  is  the  same  as  PLOT  4FM  except  it  is  computed  from 

the  head  accelerometers.  Figure  19. 

4FM  is  the  component  of  acceleration  of  the  head  C.G. 

relative  to  the  inertial  reference  coordinate 
system  origin  along  the  head  anatomical  coordinate 
+X  axis  system  as  computed  from  mouth  accelero¬ 
meters,  Figure  15. 

4SD  is  the  same  as  PLOT  4SM  except  it  is  computed 

from  the  head  accelerometers,  Figure  19. 

4SM  is  the  component  of  acceleration  of  the  head  C.G. 

relative  to  the  inertial  reference  coordinate  system 
origin  along  the  head  anatomical  +Z  axis  as 
computed  from  the  mouth  accelerometers.  Figure  15. 

4XM  is  the  component  of  acceleration  of  the  head  C.G. 

relative  to  the  inertial  reference  coordinate  system 
origin  along  the  inertial  reference  coordinate 
system  +X  axis  as  computed  from  the  mouth  accelero¬ 
meters,  Figure  15. 

4XZ  Is  magnitude  of  the  acceleration  of  the  head  C.G. 

relative  to  the  inertial  reference  coordinate  system 
origin  as  computed  from  4XM  and  4ZM,  Figure  15. 

4ZM  is  the  component  of  acceleration  of  the  head  C.G. 

relative  to  the  inertial  reference  coordinate  system 
origin  along  the  inertial  reference  coordinate  system 
+X  axis  as  computed  from  the  mouth  accelerometer. 
Figure  15, 


INTERMEDIATE  Variables  Definitions 


XVA  is  the  component  of  velocity  of  T i  anatomical 

coordinate  system  origin  relative  to  the  sled 
coordinate  system  origin  along  the  sled  coordinate 
+X  axis.  It  is  the  time  rate  of  change  of  TROUTP 
XNA,  Figure  10. 
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is  the  component  of  velocity  of  the  Tj  anatomical 
coordinate  system  origin  relative  to  the  sled 
reference  coordinate  system  origin  along  the  sled 
reference  coordinate  system  +Z  axis.  It  is  the  time 
rate  of  change  of  TROUTP  ZNA,  Figure  10. 

Definitions 


is  the  same  as  PLOT  MSP,  Figures  6  and  8, 

is  the  component  of  displacement  of  the  posterior 
(left)  spine  target  center  relative  to  the  sled 
reference  coordinate  system  origin  along  the  sled 
reference  coordinate  system  X  axis#  Figures  6  and  8. 

is  the  Z  component  associated  with  TROUTP, 

Figures  6  ana  8, 

is  the  component  of  displacement  of  the  anterior 
(right)  spine  target  center  relative  to  the  sled 
reference  coordinate  system  origin  along  the  sled 
reference  coordinate  system  X  axis.  Figures  6  and  8. 

is  the  Z  component  associated  with  TROUTP  X2, 
Figures  6  and  8. 

is  the  component  of  displacement  of  the  head  target 
center  relative  to  the  sled  reference  coordinate 
system  X  axis.  Figures  6  and  8, 

is  the  Z  component  associated  with  TROUTP  X3, 
Figures  6  and  8. 

is  the  component  of  displacement  of  the  mouth 
target  center  relative  to  the  sled  reference  coordin¬ 
ate  system  origin  along  the  sled  reference  coordinate 
system  X  axis.  Figures  6  and  8. 

is  the  Z  component  associated  with  TROUTP  X4, 
Figures  6  and  8. 

is  the  component  of  displacement  of  the  photo¬ 
graphic  head  mount  instrumentation  coordinate 
system  origin  relative  to  the  sled  reference  coor¬ 
dinate  system  origin  along  the  sled  reference 
coordinate  system  +X  axis.  Figures  8  and  DP-11 . 

is  the  Z  component  associated  with  TROUTP  X34, 
Figures  8  and  DP-1 1. 


XN 

is  the  component  of  displacement  of  the  photo¬ 
graphic  spine  mount  instrumentation  coordinate 
system  origin  relative  to  the  sled  reference 
coordinate  system  origin  along  the  sled  reference 
coordinate  system  X  axis.  Figures  8  and  DP-1 1, 

ZN 

is  the  Z  component  associated  with  TROUTP  XN, 
Figures  8  ana  DP-1 1 .  I 

XCL 

is  the  component  of  displacement  of  the  head  C.G, 
relative  to  the  sled  reference  coordinate  system 

X  axis,  Figure  8. 

ZCL 

is  the  Z  component  associated  with  TROUTP  XCL, 
Figure  8. 

XNA 

is  the  component  of  displacement  of  the  spine  T] 
anatomical  coordinate  system  origin  relative  to 
the  sled  reference  coordinate  system  origin  along 
the  sled  reference  coordinate  system  X  axis , 

Figure  8, 

ZNA 

is  the  Z  component  associated  with  XNA,  Figure  8. 

ACT 

MCT 


is  the  polar  angle  of  TROUTP  MCT  displacement- 
vector  relative  to  the  laboratory  reference  coordin¬ 
ate  system  +X  axis.  Figure  8, 

is  the  same  as  PLOT  MCT,  Figures  6,  7  and  8, 


NRN 

RN 

ADH 

PDH 

PHD  AT A  Variables 

X1,Z?,X2,Z2 
X3,  Z3,  X4,  Z4 


is  the  same  as  PLOT  NRN  except  it  is  negative. 
Figures  11  and  12. 

is  the  magnitude  of  the  angular  displacement  of  the 
photographic  spine  mount  instrumentation  coordinate 
system  +X  axis  relative  to  the  laboratory  +X  axis, 
Figures  12  and  DP-1 1. 

is  the  magnitude  of  the  angular  displacement  of  the 
photographic  head  mount  instrumentation  system 
-X  axis  relative  to  the  laboratory  reference 
coordinate  system  +X  axis.  Figures  13  and  DP-11. 

is  the  same  as  PLOT  PDH  except  it  is  negative, 
Figures  11  and  13, 


Definitions 


are  the  same  as  TROUTP  variables. of  the  same  name 
except  they  are  expressed  in  the  Tslereadex  ref¬ 
erence  coordinate  system.  Figures  12  and  13. 
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REDUCED  PHOTOGRAPHIC  Definitions 

DATA  •  ■  . . rr . 

XI,  Zl,  X2,  Z2  ore  the  same  as  PH  DATA  variables  of  the  same 

X3,  Z3,  X4,  Z4  name  except  they  are  unsealed  and  are  stored 

on  data  cards  ana  card  Image  tapes.  Figures  12 
.  and  13. 

TROUT  Variables  . . ■; "  '  Definitions  : J 

1FD  is  the  component  of  acceleration  of  the  head 

anatomical  coordinate  system  origin  relative 
to  the  inertial  reference  coordinate  systom  origin 
along  the  head  anatomical  coordinate  system  +X 
axis  as  computed  from  the  head  mount  accelero¬ 
meters,  Figure  19. 

1SD  is  the  Z  component  associated  with  TROUT  1FD, 

Figure  19. 

1FM  is  the  same  as  TROUT  1FD  except  it  is  computed 

from  the  mouth  mount  accelerometers.  Figure  15. 

ISM  is  the  same  as  TROUT  1SD  except  it  is  computed 

from  the  mouth  mount  accelerometers.  Figure  15. 

4FD,4SD,  4FM  are  the  same  as  PLOT  variables  of  the  same  name, 

4SM,4XM,4ZM  Figures  15  and  19.  ' 

QUTVAR  Variables  Definitions 

AFM  is  the  component  of  acceleration  of  the  transducer 

mouth  mount  instrumentation  coordinate  system 
origin  relative  to  the  inertial  reference  coordinate 
system  origin  along  the  transducer  mouth  mount 
instrumentation  coordinate  system  +X  axis. 

Figure  15. 

ASM  is  the  Z  component  of  acceleration  associated  with 

OUTVAR  AFM,  Figure  15. 

AFD  is  the  component  of  acceleration  of  the  transducer 

head  mount  instrumentation  coordinate  system 
origin  relative  to  the  inertial  reference  coordinate 
system  origin  along  the  transducer  head  mount 
instrumentation  coordinate  system  +X  axis.  Figure 
19. 

ASD  is  the  Z  component  associated  with  OUTVAR  AFD, 

Figure  19. 
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RPH _ _ _ _ _ Is  the  magnitude  of  the  angular  displacement  of  the 

photographic  head  instrumentation  coordinate 
system  -X  axis  relative  to  the  laboratory  reference 
coordinate  system  +X  axis.  Figure  15. 

RAH  is  the  negative  of  PLOT  RAH,  Figure  19. 


XVC,  XDC 

PSCALE  Variables 
AFM 

ASM 

— AFD 


ASD 


RVH 


AXC 


RVN 


BETHQ2  Variables 

AFM,ASM,A FD, 
ASD, RVH, RVN, 
AXC 


are  the  same  as  PLOT  variables  of  the  same  name 
except  they  are  not  in  standard  coordinate  systems, 
Figure  18.  \ 

Definitions  _ 

is  the  scaled  digital  value  of  the  t.  outh  mount 
anterior  accelerometer  output.  Figure  15. 

is  the  scaled  digital  value  of  the  mouth  mount 
superior  accelerometer  output.  Figure  15. 

is  the  scaled  digital  value  of  the  head  mount 
anterior  accelerometer  output.  Figure  19. 

is  the  scaled  digital  value  of  the  head  mount 
superior  accelerometer  output.  Figure  19. 

is  the  scaled  digital  value  of  the  mouth  mount 
rate  gyroscope  output.  It  is  the  negative  of  the 
plotted  RVH  which  has  been  converted  to  standard 
format.  Figures  15,  17  and  19. 

is  the  scaled  digital  value  of  the  sled  accelerometer, 
Figure  18. 

is  the  negative  of  the  magnitude  of  the  angular 
velocity  of  Hie  angle  of  the  transducer  Tj  mount 
coordinate  system  +X  axis  relative  to  the  laboratory 
reference  coordinate  system  +X  axis.  It  is  the 
scaled  digital  value  of  the  Tj  spine  mount  rate  gyro- 
scope  output.  Figure  17. 

■  i 

Definitions 


are  the  same  as  DSCALE  variables  of  the  same 
name  except  they  are  unsealed.  Figures  15,  17 
and  18. 
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BETHBH  Variables  ••  Definitions 

are  variables  identical  to  BETH02  variables  of  the 
vT»v.' KvH , RVN  same  name  except  they  are  in  unsorted  digital  tape 

AXC  ,  format.  Figures  15,  1£  18  and  19. 

ANALOGUE  Variables  —  Definitions  . —  - 

AFM,ASM,  AFD  are  variables  identical  to  the  DSCALE  variables 

AS D ,  RV H ,  RV  N  of  the  same  name  except  they  are  in  analogue 

AXC  voltage  format.  These  variables  are  not  illustrated 

— . . .  separately.  ""'7  . - . 

ICONDS  Variables  Definitions 

ADH  is  the  initial  value  of  TROUTP  ADH,  Figure  15. 

RN  is  the  initial  value  of  TROUTP  RN,  Figure  12. 

ALPHA  is  the  magnitude  of  the  angular  displacement  of  the 

+X  axis  of  the  Telereadex  reference  coordinate 
system  relative  to  the  +X  axis  of  the  laboratory 
reference  coordinate  system.  Figures  12  and  13. 

XI, Zl,  X2,  are  initial  values  of  PHDATA  variables  of  the 

Z2,  X3,  Z3,  same  name.  Figures  12  and  13. 

X4/  Z4 

ANXRAY  Variables  Definitions 

XAFM  is  the  component  of  displacement  of  the  center  of 

the  anterior  mouth  accelerometer  relative  to  the 
head  anatomical  coordinate  system  origin  along 
the  head  anatomical  coordinate  system  +X  axis. 
Figures  8  and  15. 

XASM  is  the  component  of  displacement  of  the  center  of 

the  superior  mouth  accelerometer  relative  to  the 
head  anatomical  coordinate  system  origin  along  the 
head  anatomical  coordinate  system  +X  axis.  Figures 
8  and  15. 

ZAFM  is  the  Z  component  associated  with  XAFM,  Figures 

8  and  15. 

ZASM  is  the  Z  component  associated  with  XASM,  Figures 

8  and  15. 
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OAFM 

XAFD 

XASD 

ZAFD 

ZASD 

OAFD 

AHS 

XAFN 

XASN 

ZAFN 

ZASN 


is  the  magnitude  of  the  angular  displacement  of  +X 
axis  of  the  transducer  moutn  mount  instrumentation 
coordinate  system  relative  to  the  +X  axis  of  the  head 
anatomical  coordinate  system  as  computed  from 
x-ray.  Figure  15. 

is  the  component  of  displacement  of  the  center  of 
the  anterior  head  accelerometer  relative  to  the  - 
head  anatomical  coordinate  system  origin  along 
the  head  anatomical  coordinate  system  +X  axis. 
Figures  8  and  19.  • —  - -  . 

is  the  component  of  displacement  of  the  center  of 
the  superior  head  accelerometer  relative  to  the 
head  anatomical  coordinate  system  origin  along 
the  head  anatomical  coordinate  system  +X  axis. 
Figures  8  a, id  19. 

is  the  Z  component  associated  with  XAFD,  Figures 
8  and  19. 

is  the  Z  component  associated  with  XASD,  Figures 
8  and  19. 

is  the  magnitude  of  the  angular  displacement  of  +X 
axis  of  the  transducer  head  mount  instrumentation 
coordinate  system  relative  to  the  +X  axis  of  the 
head  anatomical  coordinate  system  as  computed 
from  x-ray  data.  Figure  19. 

is  the  magnitude  of  the  angular  displacement  of  the 
+Z  axis  of  the  head  anatomical  coordinate  system 
relative  to  the  +X  axis  of  the  photographic  head 
instrumentation  coordinate  system  as  computed  from 
x-rays.  Figures  15  and  16. 

is  the  component  of  displacement  of  the  center  of 
the  anterior  Tj  accelerometer  relative  to  the  Ti 
anatomical  coordinate  system  origin  along  the  T] 
anatomical  coordinate  system  +X  axis.  Figure  8. 

is  the  component  of  displacement  of  the  center  of 
the  superior  Tj  accelerometer  relative  to  the  Ti 
anatomical  coordinate  system  origin  along  the  Tj 
anatomical  coordinate  system  +X  axis,  Figure  8. 

is  the  Z  component  associated  with  XAFN,  Figure  8 

is  the  Z  component  associated  with  XASN,  Figure  8 
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OAFN  is  the  magnitude  of  the  angular  displacement  of  +X 

axis  of  the  transducer  spine  mount  instrumentation 
coordinate  system  relative  to  t;he  +X  axis  of  the  Tj 
anatomical  coordinate  system  as  computed  from 
x-ray  data.  Figure  12. 

CONSTANT  Definitions 

.02  radians  fixed  angle  between  the  photo  spine  mount 

. _ “ _ _  instrumentation  coordinate  system  +X  axis  and  the 

transducer  spine  mount  instrumentation  coordinate 
system  +X  axis.  Figure  3. ,, . 

.90199  centimeters  one  half  the  distance  between  the  centers  of 

adjacent  accelerometers  within  a  given  mount. 

This  constant  was  identical  for  all  transducer 
mounts.  Figure  DP-13 , 
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ANTHROPOMETRY  CHAPTER 

Charming  L.  Ewing 
CAPT  MC  USN 


One  of  the  questions  which  continually  arises  in  human  biodynamics  research 
and  development  is  the  effect  of  differences  in  subject  size  or  of  sizes  of 
components  or  organs  (such  as  Hie  head  or  the  arm)  or  of  their  mass  distribution. 

The  experimental  design  included  consideration  of  this  question  in  several 

ways:  •« - 

1 .  Determination  of  the  effects  of  sled  peak  acceleration  and  of  rate 
of  onset  of  acceleration ,  bv  subjecting  each  volunteer  to  increments  of  peak 
acceleration  of  1G ,  while  keeping  the  rate  of  onset  constant.  Then  the  sled 
peak  acceleration  was  kept  constant  and  the  rate  of  onset  of  acceleration  was 
changed.  Thus  the  effect  of  the  acceleration  pulse  could  be  determined  com¬ 
pletely  independently  of  physical  size  for  each  volunteer,  by  each  acting  as  his 
own  control. 

2.  The  subjects  were  especially  selected  originally  because  of  their 
physical  size  in  the  parameters  of  sitting  height  alone.  (1,2,3)  Using  the  U.S. 
Navy  Anthropometric  Survey  of  1964  as  a  guide,  subjects  were  selected  roughly 
within  the  fifth,  fiftieth,  and  ninety-fifth  percentiles  of  sitting  height.  (4)  Compari¬ 
sons  between  individuals  and  between  groups  could  then  be  made  to  determine  the 
effects  of  sitting  height  on  run  results. 

3.  Comparisons  between  individuals  and  between  groups  could  be  made  to 
determine  the  effects  of  other  external  anthropometric  variables. 

4.  Determination  of  location  of  each  transducer  and  photographic  target 
relative  to  an  anatomical  coordinate  reference  system.  (See  Data  Processing  chapter) 

5.  Determination  of  the  location  of  head  C.G.  on  each  volunteer,  by 
extrapolation  from  cadaveric  measurement.  (See  X-ray  Anthropometry  data  reduction 
in  Data  Processing  chapter) 

These  headings  will  be  discussed  in  the  same  ord  jr. 

1 .  Subject  acting  as  his  own  control -observation  of  the  presented  data  will 
reveal  the  utility  or  absence  thereof  of  this  approach, 

2.  Grouping  of  volunteer  data  by  subject  size  will  reveal  the  utility  of 
this  approach  also.  This  is  facilitated  by  the  tabular  material  in  the  Results  Chapter. 

3.  Externd  Anthropometry  -  Clauser's  group  at  the  USAF  Aerospace  Medical 
Research  Laboratory  was  requested  to  perform  external  anthropometry  on  all  subjects 

in  the  study  on  completion  of  the  run  series  in  April  1969.  (5) 


50 


...  ......  ..... 


The  data  obtained  are  presented  in  the  accompanying  tables.  Table  1  lists 
all  Hie  measurements  and  the  estimated  centers  of  mass  and  the  estimated  weights. 
Table  2  lists  fust  those  measurements  applicable  to  this  report.  All  linear  and  circum¬ 
ferential  measurements  are  given  in  centimeters.  All  weights  are  given  in  grams, 
save  where  noted.  Those  values  preceded  by  *  were  calculated  according  to  the 
methods  in  reference  (6), 

The  anthropometric  techniques  used  are  those  presented  in  his  monograph.  (6) 
Those  procedures  were  described  for  use  on  cadaveric  material  and  were  modified 
by  Clauser  as  stated  below  to  measure  the  volunteers. 

A  brief  description  of  the  particular  measurements  given  in  Table  2  is  in  order. 
Head  breadth  is  obtained  by  using  spreading  calipers  to  measure  the  maximum  hori¬ 
zontal  breadth  of  the  head  (Figure  1).  Head  circumference  is  obtained  by  using  a 
flexible  tape  measure,  passed  above  the  brow  ridges  and  parallel  to  the  Frankfort 
plane,  to  measure  the  maximum  circumference  of  the  head.  (Figure  2).  The 
Frankfort  plane  is  a  horizontal  plane  passing  posteriorly  from  the  inferior  most  surface 
of  the  orbit  to  the  superior  most  surface  of  the  auditory  meatus.  (Figure  6)  Head 
length  is  obtained  by  using  spreading  calipers  to  measure  the  maximum  length  of 
the  head  between  the  glabella  and  the  occiput.  (Figure  3)  The  height  of  the  head 
is  a  derived  dimension  calculated  by  subtracting  the  chin-neck  intersect  height  from 
stature.  (Figures  4,  5)  The  chin-neck  intersect  is  defined  as  the  intersection  of  the 
following  two  planes:  (1)  from  the  menton  posteriorly,  parallel  to  the  Frankfort 
plane  and  (2)  vertically  along  the  anterior  neck  surface  upward.  (Figure  6)  The 
chin-neck  intersect  height  is  obtained  by  measuring  the  distance  from  the  floor  to  the 
chin-neck  intersect  with  the  subject  standing  erect.  (Figure  5) 

Sitting  height  is  obtained  by  measuring  the  height  from  the  sitting  surface  to  the 
top  of  the  head,  with  the  subject  sitting  erect.  (Figure  7)  Stature  is  obtained  by 
measuring  the  height  from  the  floor  to  the  top  of  the  head,  with  the  subject  standing 
erect,  (Figure  5) 

The  total  body  weight  was  obtained  by  weighing  the  nude  subject  to  the  nearest 
gram.  The  weight  in  pounds  is  listed  for  convenience  and  was  obtained  by  multi¬ 
plying  the  kilogram  weight  by  2.2. 

The  weight  of  the  head  (from  Clauser)  and  the  weight  of  the  head  and  neck 
(from  NAMRL  load  cell)  are  estimates  and  should  not  be  taken  as  exact  measured 
values.  The  value  given  by  Clauser  has  been  corrected  by  Clauser  by  multiplying 
the  value  obtained  from  the  predictive  equation  (7)  by  a  correction  factor  in  order 
to  allow  for  the  observation  that  the  sum  of  the  parts,  as  calculated  from  values 
obtained  from  the  predictive  equations,  is  slightly  greater  than  the  observed  total 
body  weight.  Clauser  attributes  the  discrepancy  to  the  loss  of  body  fluids  when  the 
cadaver  was  dismembered  prior  to  obtaining  individual  segment  weights.  The  same 
correction  technique  is  applied  to  the  live  subject  data.  It  should  be  emphasized 
that  all  the  calculated  data  presented  in  both  tables  have  been  corrected  such  that 
the  sum  of  the  parts  will  be  equivalent  to  the  measured  total  body  weight.  If  the 
regression  equations  given  in  reference  (6)  are  used  in  an  attempt  to  verify  the 
tabular  values,  an  error  will  at  first  be  apparent  unless  the  foregoing  explanation  is 
recalled. 


i  ■ 


•i  ■ 


,i  ; 

I  "i 


I  \ 

-1 


1 


51 


4.  Location  of  transducers  relative  to  anatomy  -  Inertial  sensors  mounted 
on  anatomical  locations  of  biological  specimens,  including  living  human  beings, 
measure  the  inertial  events  occurring  at  that  location  only.  To  be  able  to  generalize 
to  a  papulation  larger  than  the  instrumented  individual,  it  would  appear  desirable 
and  inaeed  necessary  to  be  able  to  specify  the  inertial  events  occurring  at  a  site 
standard  for  all  subjects.  It  is  hereby  proposed  that  such  a  procedure  might  reduce 
variance  in  the  data,  caused  by  anatomical  variation  in  the  population. 

The  standard  site  should  be  one  with  a  known  quantitative  relationship  to  the 
anatomical  location  of  the  center  of  gravity  of  Hie  instrumented  body  segment  of 
interest.  For  purposes  of  the  present  experiments,  it  is  also  important  to  be  able 
to  establish  a  Known  quantitative  relationship  to  a  coordinate  system  fixed  to  the 
body  segment  of  interest. 

investigations  of  all  cadaveric  anthropometric  data  available  in  the  literature 
indicate  that  no  head  center  of  mass  data  referred  to  a  completely  defined  anatomical 
coordinate  system  and  satisfactory  for  use  in  this  series  of  experiments  was  available. 
Clauser's  work  is  the  best  available.  (6) 

A  study  was  undertaken  at  the  request  of  the  principal  investigator  by  Walker  and 
Harris  at  Tulane  University.  (8) 

As  a  means  of  reporting  the  data  measured  by  Hie  Tulane  study,  a  series  of 
anatomical  reference  coordinate  systems  was  developed  by  the  NAMRL  research 
team.  In  order  to  make  and  represent  quantitative  measurements  in  the  mid-sagittal 
plane,  it  was  necessary  to  develop  a  system  of  x-ray  anthropometry  which  could  be 
used  to  relate  two  dimensional  mid-sagittal  plane  head  C.  G.  location  data  derived 
from  cadaveric  studies  to  similar  head  C.G.  locations  on  living  volunteers. 

The  volunteers  who  made  the  acceleration  runs  reported  herein  were  subjected  to 
x-ray  anthropometry  while  wearing  mouth  and  head  mounts,  and  the  spine  mount, 

A  lead  bar  of  known  length  was  placed  in  the  mid-sagittal  plane  coplanar  with  the 
transducer  mounts.  An  x-ray  was  then  made  of  the  head  facing  to  the  right  and  of 
the  first  thoracic  vertebra.  These  x-rays  permitted  good  visualization  of  the  external 
auditory  meati;  the  infraorbital  ridges;  the  angular  and  linear  location  of  the 
anatomical  and  transducer  mounts  relative  to  alt  observable  anatomical  features;  and 
a  good  measure  of  skin  thickness  between  anatomical  mounts  and  the  nearest  body 
point.  These  features  could  be  visualized  at  leisure  and  without  injury  or  discomfort 
to  the  subject.  Furthermore,  for  these  experiments,  a  three  dimensional  problem  was 
reduced  to  a  two  dimensional  one. 

The  anatomical  coordinate  systems  were  established  by  selecting  anatomical 
features  which  would  be  visible  by  x-ray  and  which  would  be  as  rigid  as  possible, 
dynamically  or  statically.  The  anatomical  features  also  had  to  be  easily  locatable 
by  the  same  means  in  cadavers  and  should  not  change  due  to  death  or  cadaveric 
preservation  techniques. 
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For  the  head,  the  Frankfort  plane  was  selected  as  the  basis  of  such  a  coordinate 
system.  For  TI  the  anterior  superior  corner  of  the  vertebral  body  and  the  superior 
terjor  and  inferior  posterior  corners  of  the  posterior  spinous  process,  as  seen  on 
id.  (See  Fiaures  3  and  4  of  Theoretical  Kinematics  Char 


The  x-ray  anthropometry  procedure  was  as  follows: 

a.  Fix  the  anatomical  and  transducer  mounts  to  the  head  and  the  spine 

as  done  for  each  acceleration  run.  Repeatability  of  mount  positioning  was  established 
by  the  use  of  unique  individual  anatomical  mounts  for  each  subject  as  previously 
explained.  -  - 

b.  Small  spherical  lead  anatomical  markers  were  placed  at  the  superior 
edge  of  each  auditory  meatus  and  at  the  infraorbital  notches,  as  shown  in  Figure  8, 

c.  A  precision  lead  linear  calibration  rule  was  placed  in  the  mid-sagittal 
plane  as  shown  in  Figure  8  and  Figure  9. 

d.  A  lateral  x-ray  of  the  skull  including  head  mount  and  mouth  mount 
screws  and  the  lead  anatomical  markers  was  taken.  Unfortunately,  reproduction  of 
these  x-rays  does  not  permit  demonstration  of  the  head  and  mouth  screws  which  are 
very  small. 

e.  A  laterql  x-ray  of  the  cervical  spine  including  the  spine  mount,  screw 
holes  and  the  first  thoracic  vertebra  was  taken  next,  as  shown  in  Figure  9.  Unfor¬ 
tunately,  reproduction  of  these,  x-rays  does  not  permit  demonstration  of  the  T] 
vertebral  body  • 

The  x-rays  were  acceptable  only  if:  the  calibration  device,  the  anatomical 
markers  or  features  and  the  mount  screw  holes  or  screws  were  clearly  visible  with 
good  resolution.  Occasionally  this  required  repetition  of  an  x-ray  but  every  effort 
was  made  to  minimize  exposure  for  eacn  volunteer.  Acceptable  x-rays  were  then 
submitted  for  data  reduction.  (See  X-ray  Anthropometry  Data  Reduction  in  Data 
Processing  Chapter) 

The  anatomical  coordinate  systems  which  were  then  established  using  these 
data  points  are  described  under  mat  heading  in  the  Theoretical  Kinematics  section  . 
The  head  coordinate  system  is  illustrated  in  Figure  2  and  the  spine  coordinate  system 
is  Figures  3  and  4  of  tnat  section. 

Once  the  coordinate  system  was  established  on  the  x-rays,  the  precise  location 
of  each  of  the  transducers  within  it  could  be  determined.  Furthermore,  the  heed 
C.G.  could  also  be  established  within  the  head  anatomical  coordinate  system,  thus 
permitting  transformations  of  the  transducer  data  to  the  C.G.  using  rigid  body 
mechanics. 

5.  Location  of  the  head  center  of  gravity  for  each  volunteer. 

Head  center  of  gravity  measurements  on  cadavers  were  undertaken  in  1968  by 
NAMRL  because  of  the  paucity  of  useable  data  available  at  that  time.  A  study 
performed  by  Tulane  University  on  sixteen  cadavers  determined  the  head  C.G. 


locations  within  the  head  anatomical  coordinate  system,  (8)  The  measurements  which 
are  previously  unpublished  are  presented  in  Figure  10.  This  study  will  be  reported  in 
detail  elsewhere. 

Another  study  performed  by  Becker,  on  five  cadavers,  also  previously  unpublished, 
gave  comparable  results  which  are  also  presented  in  Figure  10.  (9) 

One  of  the  findings  both  by  the  Tulane  group  and  by  Becker  was  that  no  useable 
data  exist  concerning  head-neck  C.G.  location.  Therefore,  all  transformations 
in  this  report  used  the  averaged  head  C.G.  location,  as  determined  by  the  Tulane 
measurements  because  of  the  larger  population  size  and  presented  in  Figure  10. 

While  It  would  appear  much  more  realistic  to  use  head-neck  C.G.,  quantitatively 
useable  data,  as  noted,  are  not  known  to  exist. 

It  is  of  some  interest  to  note  the  small  variance  of  head  C.G.  locations  when 
expressed  in  the  head  anatomical  coordinate  system,  which  also  is  previously  unpub¬ 
lished.  It  would  appear  that  such  a  coordinate  system  may  be  of  considerable  utility 
in  biodynamics  research  in  view  of  the  fortuitous  location  of  its  origin  near  the 
measured  head  C.G. 

The  averaged  data  from  the  Tulane  study  for  sixteen  head  C.G.  locations  in  the 
head  anatomical  coordinate  system  are: 


X  =  +0.883  cm,  Z  =  +2.145  cm  (8) 

The  averaged  data  from  Hie  N.AMRL  study  for  five  head  C.G.  locations  in  the 
head  anatomical  coordinate  system  are: 

X  -  +1.33  cm,  Z  =  +2.09  cm  (9) 
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FIGURE  10  Center  of  Gravity  Expressed  in  the  Head 
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DATA  PROCESSING  CHAPTER 
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SENSOR  DATA  PROCESSING 

Sensor  data  processing  began  with  the  creation  of  analog  tapes  containing,  for 
each  transducer  channel, nine  calibration  steps  of  known  voltage  levels  provided  by 
a  calibrator-amplifier, and  the  voltage  output  from  the  given  transducer. 

These  tapes  were  recorded  at  the  experimental  site  at  60  ips  on  an  Ampex  FR  1800L 
fourteen  channel  FM  instrumentation  recorder,  which  had  an  FM  flat  band-pass 
of  DC  to  20  kHz.  After  recording,  the  tapes  were  transported  in  cases  shielded 
against  magnetism  to  NAMRL,  Pensacola,  where  an  identical  Ampex  FR  1800L  was 
utilized  for  playback. 

Servo  control  of  tape  motion  was  employed  in  all  runs  to  insure  that  the  tape  on 
the  playback  FR  1800L  would  be  completely  synchronized  with  the  tape  speed 
variations,  if  any,  of  the  recording  FR  1800L.  Since  the  FR  1800L  uses  IRIG 
Standard  tape  speed  control  frequency  of  100  kHz  for  60  ips  tape  speed,  it  was 
decided  to  use  the  control  track  frequency  as  the  master  timing  clock  in  A/D  con¬ 
version  also.  Playback  of  the  analog  tape  was,  however,  accomplished  at  a  tape 
speed  of  15  ips;  therefore,  the  control  track  frequency  (master  clock)  became  25  kHz. 
This  four-fold  speed  reduction  was  employed  to  facilitate  analog  strip  chart  monitoring 
of  playback  and  to  insure  remaining  within  the  sampling  rates  available  to  the  A/D 
converter. 

At  the  recording  (experimental)  site  an  analog  signal  calibrating  unit  was  employed 
to  provide  the  steps  of  known  precision  voltage  levels  from  -1 .4  to  +  1 .4  volts  which 
were  applied  to  each  sensor  channel.  Upon  completion  of  the  calibration  steps, 
recorder  input  channels  were  returned  to  the  sensor  channels  and  their  respective 
signal  conditioning  amplifiers. 

For  the  duration  of  the  calibration  steps,  a  1.4  volt  control  gate,  DC  level  shift, 
(Figure  DP-2)  was  applied  to  one  track  of  the  recording  FR  1800L.  Similarly,  when 
the  output  from  the  transducers  was  recorded  and  for  a  short  time  before  each  run,  a 
second  control  gate  was  applied  to  this  same  channel .  The  second  gate  was  concluded 
at  the  end  of  the  run.  These  two  control  gates  provided  the  triggering  signal  for  the 
A/D  sampling  periods  and  rates. 

The  playback  FR  1800L  was  physically  located  remote  from  the  A/D  converter  and 
computer  facility  at  NAMRL,  Pensacola.  Fifty  ohm  coaxial  transmission  cables 
approximately  150  feet  long  were  required  for  playback  of  the  analog  data  to  the 
A/D  converter.  Amplifiers  for  each  playback  channel  were  employed  to  interface 
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the  playback  FR  1800L  with  these  transmission  lines.  These  amplifiers  had  provision 
for  adjusting  the  analog  signal  level  in  gain  increments  of  0,5,  1 .0  and  2.0.  This 
interfacing  unit  also  provided  the  means  of  using  the  25  kHz  servo  control  signal 
(master  clock)  to  direct  sampling  of  the  A/D  data.  Sampling  rates  were  triggered 
by  the  gates  contained  on  the  control  track.  (See  Figure  DP-0) 

The  precision  of  the  servo  control  track  (master  clock)  and  of  the  control  gates  was 
such  that  repeated  digitizing,  over  a  period  of  two  years,  of  the  same  run  resulted 
in  sampling  of  the  same  portions  of  the  analog  signal.  This  was  demonstrated  by 
comparison  plots  of  the  sorted  data. 

The  band-pass  of  these  interfacing  amplifiers  was  DC  to  at  least  10  kHz,and  the 
gain  increments  were  maintained  to  within  a  few  millivolts. 

The  A/D  converter  offered  the  option  of  digitizing  7  to  14  channels.  Since  the 
experimental  analog  tapes  contained  up  to  10  channels  of  sensor  data,  the  14  channel 
option  was  elected. 

The  A/D  converter  digitized  each  channel  sequentially  and  thus  introduced  a 
small  time  skew  of  one  channel  relative  to  another.  The  time  required  to  sample  an 
individual  channel  was  37.5  microseconds.  This  introduced  a  maximum  time  skew 
across  14  channels  of  approximately  500  microseconds.  Since  the  playback  speed 
was  1/4  the  data  acquisition  speed,  the  maximum  possible  time  skew  introduced  was 
125  microseconds.  Compared  to  the  basic  band-pass  characteristics  of  the  transducers, 
this  is  considered  to  be  an  insignificant  source  of  error  in  the  data. 

The  first  step  in  the  A/D  procedure  consisted  of  a  visual  evaluation  of  the  data 
on  each  channel.  The  analog  tape  was  played  backend  strip  charts  were  made. 

These  strip  charts  were  carefully  examined  to  determine  how  many  samples  of  cali¬ 
bration  and  data  would  be  taken  to  insure  adequate  data  acquisition.  In  addition, 
a  check  was  made  to  determine  if  the  signal  on  any  channel  exceeded  the  upper  and 
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at  the  time  of  A/D  conversion. 

Sampling  rates  were  set  at  least  10  times  the  highest  anticipated  data  frequency: 
200  samples/second  for  calibration  and  2000  samples/second  for  data  at  the  time  of 
the  experiment.  With  playback  speed  1/4  of  record  speed,  these  rates  became 
50  samples/second  and  500  samples/second  respectively. 

An  assembly  language  program  on  the  UNIVAC  418  digital  computer  was  used  to 
collect  data  from  the  A/D. 


FIGURE  DP-0 


BETHBH  -  The  A/D  Conversion  Program 

BETHBH  is  an  assembly  language  program  which  collects  digital  data  from  an 
A/D  converter,  stores  it  on  a  FASTRAND  drum  and  then  records  on  digital  tape  a 
permanent  file  of  the  digitized  data  for  the  run  being  processed.  Because  a  large 
amount  of  data  is  arriving  in  real  time,  c  double  buffered  procedure  is  used  to  allow 
continuity  of  input  and  to  facilitate  immediate  data  transfer  onto  the  drum.  BETHBH 
does  not  alter  the  data  in  any  way;  it  serves  only  as  a  data  transfer  operation. 

Input  to  BETHBH  is  the  digital  output  from  the  A/D  converter  and  a  card  containing 
the  following  parameters. 

1 .  The  run  number 

2.  Output  tape  name 

3.  Number  of  channels  to  be  converted 

4.  Number  of  records  of  calibration  data 

5.  Number  of  samples/ecord 

6.  Number  of  records  of  data 

7.  Number  of  samples/record 

8.  Number  of  seconds  to  wait  between  calibration  A/D  conversion  and 
data  A/D  conversion. 

This  input  is  utilized  in  the  following  A/D  procedure: 

1.  The  BETHBH  program  is  loaded  on  the  UNIVAC  418. 

2.  The  analog  tape  is  loaded  on  the  FR  1800L  tape  recorder. 

3.  A  switch  is  set  on  the  UNIVAC  418  which  will  enable  the  computer  to 
accept  digitized  data  as  soon  as  it  is  available. 

4.  The  analog  tape  is  started.and  the  first  control  gate  appears.  The  A/D 
converter  receives  the  clock  signal  which  triggers  sampling  and  digitizing  of  the 
calibration  steps  at  50  samples/second .  The  computer  accepts  from  the  A/D  the 
number  of  samples  of  calibration  data  for  each  channel  which  were  specified  in  the 
input  card. 

5.  The  computer,  having  satisfied  its  requirements  for  channel  calibration 
data,  stays  away  a  prescribed  number  of  seconds.  The  first  control  gate  ends  and 
the  A/D  suspends  sampling. 

6.  The  second  control  gate  triggers  sampling  and  digitizing  at  the  rate  of 
500  samples/second.  The  computer  accepts  from  the  A/D  the  number  of  samples  for 
each  channel  which  was  specified  in  the  input  card. 

7.  As  a  final  step  BETHBH  writes  an  output  tape  containing  all  of  the  data 
collected  on  the  drum. 

BETHBH  is  a  multi-run  program  enabling  th’s  procedure  to  be  followed  on  up  to 
20  runs  without  reloading  the  program. 
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BETH8H  Validation 

As  a  validation  of  the  A/D  conversion  system,  two  completely  arbitraty  test  tapes 
were  digitized  at  regular  intervals  (approximately  every  lu  runs)  throughout  a  given 
group  of  runs.  The  first  of  these  test  tapes  contained  triangular  waves  of  known 
amplitudes.  This  tape  was  digitized,  sorted  and  plotted.  The  plotted  data  were  then 
compared  to  analog  strip  charts  of  the  test  tape.  A  second  tape  contained  alternating 
zero  and  full  scale  DC  levels.  It  was  digitized,  sorted  and  processed  by  a  statistical 
program  (STATAD)  which  calculated  the  mean,  variance  ana  standard  deviation  of 
the  levels.  Evaluation  of  the  output  from  these  two  test  tapes  gave  a  very  precise 
measure  of  the  performance  of  the  system.  If  an  indication  of  impaired  performance 
was  detected,  all  digitized  data  subsequent  to  the  last  successful  run  of  the  test  tapes 
was  discarded.  The  system  was  then  checked /md  the  hardware  or  software  defects 
causing  the  problem  were  determined  and  corrected.  The  test  tapes  were  then  run 
again.  If  they  were  found  to  measure  satisfactory  performance  of  the  system,  the  A/D 
process  was  then  resumed.  The  test  tape  data  after  A/D  conversion,  sorting  and 
plotting  were  retained  as  reference  data. 

The  data  were,  at  the  conclusion  of  the  operation  of  the  program  BETHBH,  resident 
in  digital  form  on  digital  magnetic  tape,  unsorted  as  to  channel.  A  separate  tape 
was  used  for  each  run  and  these  tapes  were  and  are  maintained  as  a  primary  reference. 
Furthermore,  the  original  analog  tape  for  every  run  is  presently  preserved  in  cabinets 
which  are  protected  against  magnetism,  dust,  humidity  and  temperature  extremes. 

It  was  then  necessary  to  sort  the  digital  data  by  channel,  which  was  accomplished 
by  operation  of  the  BETH02  program  on  the  data. 
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BETH02  -  The  Data  Sort  Program 

BE'iH02  is  a  FORTRAN  program  that  sorts  and  reorders  14  channel  A/D  output 
data  as  specified  by  the  user. 

Input  to  BETH02  consists  of  the  output  tape  from  A/D  conversion  and  parameter 
cards  specifying  which  channels  contain  data,  the  desired  order  for  the  data,  and 
the  variable  name  to  be  assigned  to  each  channel .  The  A/D  output  tape  contains 
data  ordered  exactly  as  converted;  that  is,  first  data  value  for  each  of  fourteen 
channels,  second  data  value  for  each  of  fourteen  channels  and  so  on.  This  input  is 
utilized  in  the  following  sort  procedures 

1.  Each  data  value  on  the  input  tape  contains  channel  identification 
information  in  the  high  order  six  bits  and  digitized  data  in  the  low  order  twelve  kits. 
BETH02  drops  off  the  channel  identification  information. 

2.  Using  the  channel  identification  information,  BETH02  sorts  the  data. 

All  fourteen  channels  are  sorted  regardless  of  whether  they  contain  data. 

3.  The  sorted  data  are  stored  on  the  drum. 

4.  The  sorted  data  for  the  channels  which  contain  data  are  pulled  from  the 
drum  in  the  specified  order,  and  variable  names  are  added. 

5.  An  output  tape  is  written  containing  the  sorted  data  in  the  specified  order. 
BETHBH  and  BETHQ2  Validation 

As  a  validation  of  the  sort  procedure  and  further  verification  of  the  A/b  procedure, 
all  of  the  BETH02  output  channels  were  plotted.  A  scaling  factor  of  0.00069053396 
was  used  on  the  data  which  transformed  each  data  value  from  analog-to-digital  units 
to  voltage,  for  plotting  purposes. 

This  conversion  permitted  the  comparison  of  B0THO2  plots  of  the  data  (Figure  DP-1) 
with  the  analog  strip  chart  plots  of  the  same  data  (Figure  DP-2)  generated  from  the 
original  analog  tape  at  the  time  of  A/D  conversion.  These  comparisons  of  each 
channel  of  each  run  are  on  file  and  show  no  degradation  of  the  analog  signal  by  the 
A/D  conversion  process  or  by  the  sorting  process.  It  is  considered  that  this  is  tne  best 
possible  verification  of  the  absence  of  artifact  due  to  hardware  and/or  software  to 
this  point.  If  a  channel  was  noted  to  have  defective  data  for  any  reason,  it  was 
excluded  from  further  processing  until  the  cause  of  the  defect  could  be  determined. 

Several  runs  suffered  from  sudden  transducer  failure  midway  through  the  run. 

While  no  such  runs  were  used  in  this  report,  options  3  or  4  of  program  OUTVAR 
could  be  used  to  yield  perfectly  valid  results  due  to  the  considerable  redundancy  of 
the  instrumonration  system  design. 

A  single  digital  tape  was  used  to  record  the  BETH02  output  from  each  run  and 
these  tapes  were  and  are  maintained  as  a  primary  reference. 

The  datu  were  then  subjected  to  a  scaling  program,  DSCALE. 
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DSCALE  -  The  Transducer  Scaling  Program 

DSCALE  is  a  FORTRAN  program  which  performs  three  major  functions.  The  first 
function  of  DSCALE  is  to  locate  and  identify  points  of  experimental  interest  within 
the  input  data  stream.  These  points  of  interest  are: 

1.  Tq  ,  the  time  of  a  specified  pulse  of  known  amplitude  and  duration  on  the 
sled  accelerometer  channel.  Tq  provides  a  means  for  determining  the  location  of  any 
frame  of  film  relative  to  the  time  of  occurrence  of  the  sensor  data  .(See  Photographic 
Data  Processing) 

2.  The  peak  sled  acceleration. 

3.  DPTZ,  a  computer  defined  time  zero  which  represents  the  beginning  of 
the  data  of  experimental  interest. 

DSCALE's  second  function  is  to  convert  the  data  from  the  system  of  units  which  is 
the  analog-to-digital  converter's  representation  of  voltage  to  physical  units  of 
acceleration  and  angular  velocity  (cm/m secs  and  radians/msec).  This  conversion 
involves  linear  least  squares  analyses  of  calibration  data  for  the  tape  recorder  - 
A/D  system,  the  range  amplifiers,  and  the  accelerometers  and  rate  gyros. 

After  the  data  a.o  converted,  the  last  function  of  DSCALE  is  to  edit  the  data  and 
to  save  only  that  portion  of  experimental  significance. 

The  input  to  DSCALE  is  a  digital  tape  containing  nine  calibration  steps  of  known 
voltage  level,  and  data  from  each  of  tne  accelerometer  and  rate  gyro  channels. 

This  tape  is  the  sorted  output  data  from  the  analog-to-digital  conversion.  Transducer 
output  is  expressed  in  a  system  of  units  which  is  the  A/D  converter's  representation 
of  voltage.  Additional  input  includes  identification  specifications  which  define  TQ 
and  DPTZ;  the  high,  low  and  null  input/output  voltages  for  each  range  amplifier;  ° 
and  the  manufacturer's  calibration  for  each  transducer.  DSCALE  uses  this  input  as  it 
does  the  following: 

1.  All  of  the  transducer  output  data  are  multiplied  by  0.00069053396.  This 
transforms  the  data  from  analog-to-digital  units  to  voltage. 

2.  The  calibration  steps  for  each  analog  recorder  channel  are  examined,  and 
step  averages  are  computed  for  each  channel.  An  examination  is  made  to  determine 
if  each  step  average  lies  within  predetermined  limits*  Steps  which  do  not  fall  within 
these  limits  are  deleted.  The  deletion  of  more  than  six  steps  for  any  channel  aborts 
the  run.  All  valid  step  averages  are  maintained  for  future  usage. 

3.  The  sled  accelerometer  data  is  examined  and  T  is  located.  T  is  defined 

o  o 

to  be  the  first  data  point  of  a  pulse  that  has  a  minimum  amplitude  of  -0.5  volts  and 
a  duration  of  3.5  msec. 
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4.  The  sled  peak  acceleration  is  located  by  examining  a  10  point  running 
average  of  all  of  the  sled  acceleration  data  after  Tq.  The  point  of  peak  acceleration 
is  where  this  running  average  reaches  its  mcoo'mum. 

5.  A  backward  running  average  is  begun  from  the  sled  peak  acceleration. 
When  this  average  falls  below  40  percent  of  the  peak  acceleration,  the  averaging 
stops  and  30  data  points  are  stepped  back.  The  data  value  reached  is,  by  input 
definition.  Data  Processing  Time  Zero  (DPTZ).  DPTZ  is  thus  defined  to  be  the  data 
point  which  occurs  15  msec  prior  to  the  time  at  which  40  percent  of  the  averaged 
sled  peak  acceleration  is  reached. 

6.  For  each  channel,  the  32  data  points  immediately  before  DPTZ  are 
averaged  to  get  a  number,  BASE,  which  is  utilized  in  rate  gyro  scaling. 

7.  A  least  squares  linear  fit  is  executed  for  each  channel  using  the  calibration 
step  averages  previously  calculated  as  the  dependent  values,  and  input  voltages  to 
the  calibrator  as  independent  values.  The  calibrator  input  voltages  were  supplied  on 
log  sheets  which  accompanied  each  analog  tape  and  were  incorporated  into  punch 
cards  as  input  to  DSCALE.  The  slope  (B2)  and  intercept  (B1)  from  the  least  squares 

fit  for  a  given  channel  are  used  to  transform  each  data  point  (D)  of  that  channel . 

A  =  (B2)(D)+B1  (1) 

where  A  is  the  corrected  data  point. 

8.  A  second  least  squares  fit  for  each  channel  is  used  to  incorporate  the 
range  amplifier  settings.  The  high,  low  and  null  input/output  voltages  for  each  range 
amplifier  were  supplied  on  log  sheets  which  accompanied  each  analog  tape  and  were 
incorporated  into  punch  cards  as  input  to  DSCALE.  The  resulting  slope  (B6)  and 
intercept  (B5)  are  used  to  transform  each  corrected  point  (A)  of  a  given  channel . 

S  =  (B6)(A)+B5  (2) 

where  S  is  the  corrected  data  point. 

9.  A  final  least  squares  fit  is  done  to  convert  voltages  to  accelerations  and 
angular  velocities.  This  fit  uses  the  appropriate  manufacturer's  calibration  for  a 
given  transducer.  The  manufacturer's  calibration  was  incorporated  into  punch  cards 
as  input  to  DSCALE.  The  slope  (B8)  and  intercept  (B7)  of  this  final  fit  are  applied 
to  each  corrected  data  point  (S)  of  a  given  transducer. 

F  =  (B8)(S)  +  B7  (3) 

where  F  is  a  completely  scaled  data  point. 

10.  The  above  described  transformation  from  D  to  F  is  executed  in  one  step. 

It  may  readily  be  seen  that  given: 

A  *  (B2)(D)+B1 
S  =  (B6)(A)  +  B5 
F  =  (B8)(S)+B7 
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(4) 

(5) 

(6) 


A  final  slope  (B12)  and  intercept  (Bll)  may  be  calculated* 

B12  =  (B8)(B6)(B2)  (7) 

Bll  =  (B8)(B6)(B1)  +  (B8)(B5)+B7  (8) 

The  total  data  transformation  for  any  data  point  D  of  a  given  transducer  is 

F  =  Bll -BASEL +  (B12)(D)  (9) 

where  BASEL  is  an  offset  constant. 

BASEL  =  O  for  accelerometers 

BASEL  =  Bll  +  (B12)  (BASE)  for  rate  gyros  (10) 


where  BASE  was  previously  defined  as  being  the  average  of  32  data  points 
immediately  preceding  DPTZ  and  BASEL  is  the  scaled  BASE. 

11.  Each  scaled  accelerometer  data  point  is  multiplied  by  0.000980616  to 
convert  it  from  G's  to  cm/msecs  and  the  scaled  rate  gyro  data  is  multiplied  by 
0.00001745329  to  convert  it  from  degrees/sec  to  radians/msec. 

12.  The  final  action  of  DSCALE  is  to  create  a  digital  output  tape  containing 
1,000  corrected  data  values  for  each  transducer.  The  output  for  each  transducer 
begins  with  data  point  DPTZ  and  consists  of  the  data  for  500  msec  after  DPTZ. 

DSCALE  Validation  ' 

As  a  validation  of  the  scaling  procedure  and  the  computer  program,'  all  of  the 
DSCALE  output  was  plotted  by  a  digital  computer -driven  Stromberg-Carlson  4020 
plotter  (Figures  DP-3  and  DP-4)  ana  carefully  compared,  channel  by  channel,  with 
its  corresponding  BETH02  plot  (Figure  DP-l)  and  analog  strip  chart  (Figure  DP-2). 
Plots  showing  discrepancies  were  investigated  and  the  cause  determined  and  corrected. 
These  plots  were  ana  are  maintained  on  file. 

Some  of  the  plots  presented  herein  were  originally  designed  to  verify  transducer 
and  photographic  scaling.  An  examination  of  the  PDH  (derived  from  tne  head  angular 
position  as  computed  from  photographic  data)  vs  RDH  (derived  from  the  scaled  head 
rate  gyro)  comparison  plots  (Plots  633  through  673)  illustrate  the  precision  and 
accuracy  of  sensor  and  photographic  scaling. 

A  single  digital  tape  was  used  to  record  the  DSCALE  output  from  each  run  and 
these  tapes  were  and  are  maintained  as  a  primary  reference. 

Figure  DP-5  illustrates  the  flow  of  sensor  data  to  a  DSCALE  output  tape  compatible 
with  the  UNIVAC  1108  computer  system. 
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SCALED  HEAD-MOUNTED  SENSOR  DATA  DSCALE 
SUBJECT  03  9.5G  358G/SEC 


While  the  sensor  data  were  being  digitized,  sorted  and  scaled,  a  simultaneous 
effort  was  underway  to  acquire  scaled  photographic  data  (PHD  AT  A)  and  to  establish 
the  x-ray  anthropometry  and  initial  conditions  necessary  for  subsequent  processing  4 
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PHOTOGRAPHIC  DATA  PROCESSING 


The  filming  of  each  human  exposure  to  impact  acceleration  began  before  the 
placement  of  the  subject  on  the  sled  and  continued  until  the  end  of  the  experiment. 
The  sequence  of  events  filmed  may  be  classified  in  four  areas: 

1.  An  identification  board  was  placed  in  the  sled-mounted  camera  field 
of  view.  This  identification  board  contained  the  sled  run  number/  the  human  run 
number,  the  subject  identification  number,  the  nominal  G  level  and  rate  of  onset, 
the  camera  framing  rate,  and  the  Julian  date,  along  with  accumulator  pressure 
and  volume. 

2.  An  inverse  L-shaped  calibration  device  which  had  a  tripod  mounting 
which  fit  into  precisely  located  and  machined  pin  holes  on  the  slea  was  mounted. 

The  design  was  such  that,  when  placed,  the  calibration  device  was  located  in  a 
plane  corresponding  to  the  mid-sagittal  plane  of  the  volunteer  subject  and  had 
straight  lines  perpendicular  and  parallel  to  gravity.  The  centers  of  the  targets  which 
defined  the  lines  were  12  inches  apart.  The  cameras  were  actuated  momentarily  to 
record  the  board  and  calibration  device  on  the  beginning  of  each  data  film. 

Figure  DP-6  illustrates  such  a  frame  of  photo  data  calibration.  The  plumb  line  used 
to  verify  the  calibration  device  is  visible  at  right. 


FIGURE  DP-6  A  Frame  of  Film  Demonstrating  the  I  .D.  Board 
and  Calibration  Device 
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3.  The  subject  was  then  placed  on  the  sled  and  positioned  prior  to  the  onset 
of  acceleration  (Figure  DP-7) . 


FIGURE  DP-7  A  Frame  of  Film  Demonstrating  the  Volunteer's 
Position  Prior  to  the  Onset  of  Acceleration 


hooH  nri ’n0A£CelerfMi^ on  was  initiated,  and  the  resulting  response  of  the  volunteer's 
during  a^run^  WQS  me^*  ^'9ure  DP”8  shows  the  point  of  maximum  neck  extension 
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FIGURE  DP-0  A  Frame  of  Film  Demonstrating  the  Resulting 
Response  of  the  Volunteer's  Head  and  Neck 

The  arbitrary  Tq  pulse  on  the  sensor  data  AXC  channel  gave  a  means  for 
determining  the  location  of  any  frame  of  film  relative  to  the  time  of  occurrence 
of  the  sensor  data  (since  the  pulse  could  be  precisely  located  in  time  in  each 
sensor  channel)  and  simultaneously  reset  the  camera  timers  to  0.  Therefore,  the 
first  non-zero  camera  timer  reading  occurred  within  2  msec  of  the  Tq  pulse  for 
frame  rates  of  500  F.P.S.  which  were  standard  for  all  runs.  First  motion  was  marked 
by  the  appearance  of  a  light  which  was  a  flashbulb  triggered  by  the  closing  of  a 
microswitch  as  the  ram  left  the  resting  position  (See  Figure  DP-8  just  in  front  of  the 
forehead) . 

After  the  completion  of  each  run,  the  exposed  film  was  mailed  to  the  Photographic 
Laboratory  at  Eglin  AFB  for  precision  development  and  printing.  Film  from  both  the 
posterior  and  lateral  sled-mounted  cameras  was  available  on  most  runs.  On  some 
runs,  a  supplementary  lateral  camera  provided  color  coverage  of  the  experiment  at 
varied  frame  rates. 

The  developed  film  negatives  were  then  forwarded  to  the  Mathematical  Services 
Laboratory  at  Eglin  for  data  reduction. 

Photographic  data  reduction  was  accomplished  by  use  of  Benson  Lehner  Telereadex 
Type  29E  equipment  interfaced  with  an  IBM  523  summary  punch.  When  the  horizontal 
and  vertical  cross-wires  of  the  device  were  aligned  so  that  their  intersection  was  on 
a  desired  point,  analog  voltages  were  generated  which  represented  the  coordinates 
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of  that  point  relative  to  any  other  point  selected  as  origin*  All  targets  were 
measured  relative  to  the  sled  target  origin  on  data  runs.  These  voltages  were  then 
digitized  and  the  digitized  coordinates  were  punched  on  standard  80  column  cards. 
Table  1  presents  card  format  and  Table  2  and  Figure  DP-9  present  data  point 
identification.  One  card  was  punched  per  point.  This  card  contained,  in  addition 
to  the  coordinates,  point  identification  and  information  obtained  from  the  identifi¬ 
cation  board  seen  at  the  beginning  of  each  film.  This  latter  information  was  input 
to  the  summary  punch  from  the  Telereadex  key  board  and  a  set  of  constant  input 
switches. 

The  specific  procedure  utilized  by  the  Mathematical  Services  Laboratory  was 
as  follows: 

1 .  The  film  negative  was  mounted  and  projected  on  the  operator's  screen. 

2.  The  identification  board  was  read  and  required  information  input  into 
the  machine  via  the  Telereadex  key  board  and  constant  switches. 

3.  Ten  calibration  frames  containing  the  calibration  device  were  reduced. 

A  time  code  of  600000  was  entered  for  all  calibration  frames.  The  sled  target  origin 
was  the  X  and  Z  zero  point. 

4.  The  film  was  then  advanced  to  the  first  frame  with  a  human  subject  and 
a  non-zero  time  reading.  This  frame  corresponds  to  the  first  frame  of  film  after  the 
Tq  pulse,  some  time  prior  to  acceleration. 

5.  The  film  was  then  advanced  to  the  first  frame  in  which  the  first  motion 
light  appeared. 

6.  The  film  was  then  retreated  25  frames,  at  which  point  data  reduction 
began.  The  sled  target  origin  was  the  X  and  Z  zero  point.  The  six  digit  camera 
timer  time  was  read  on  each  data  frame. 

7.  Data  reduction  of  every  frame  was  continued  for  300  frames. 

The  above  procedure  described  was  designed  to  insure  that  sufficient  precise 
data  were  reduced.  The  film  negatives  were  carefully  stored  in  individual  cans  for 
each  run  and  were  and  are  maintained  in  a  controlled  environment  as  a  primary 
reference. 

The  computer  analysis  of  photographic  data  began  with  the  receipt  of  the  reduced 
data  cards  in  Pensacola.  Computer  programs  were  then  utilized  to  select  the  pre¬ 
acceleration  data  necessary  for  determination  of  initial  conditions  and  to  determine 
which  portion  of  the  reduced  data  was  of  experimental  interest.  The  ten  frames  of 
cal  (oration  data  and  ten  frames  of  pre-acceleration  data  were  used  to  develop 
horizontal  and  vertical  scaling  constants  (see  program  ICONDS).  These  scaling 
constants  and  the  remaining  reduced  photographic  data  were  input  to  a  photographic 
scaling  program,  PHDATA,  and  were  used  to  accomplish  the  following: 

1 .  The  X  and  Z  coordinates  of  each  point  from  each  frame  and  the  corres¬ 
ponding  frame  time  were  read  into  digital  computer  storage  arrays. 
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2.  A  check  was  made  to  eliminate  any  data  having  a  camera  timer  reading 
prior  to  DPTZ  (see  program  DSCALE).  This  insures  the  same  beginning  point  for 
sensor  and  photographic  data. 

3.  A  check  was  made  for  missing  or  out  of  order  data  cards.  Missing  or  out 
of  order  data  cards  aborted  the  run . 

4.  The  horizontal  and  vertical  scaling  constants  were  applied  respectively 

to  each  of  the  X  and  Z  coordinates.  # 

5.  An  output  digital  tape  was  written  containing  the  X  and  Z  scaled 
coordinates  and  the  corresponding  times.  These  tapes  as  well  as  all  the  photo  data 
cards  were  and  are  maintained  as  a  primary  reference. 

No  interpolation  was  done  for  photographic  data.  When  a  photographic  target 
was  temporarily  obscured  by  subject's  anatomy  or  for  any  other  reason,  the  readings 
for  all  frames  affected  were  deleted.  This  can  be  observed  in  the  plots  of  PDH  vs 
ROH.  Plots  639,  641,  644,  670,  671  and  672  are  examples. 
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Subject  Identification 

ID  Board 
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Telereadex 

TABLE  2 


PHOTOGRAPHIC  DATA  POINT  IDENTIFICATION 


CALIBRATION 


Origin  (0,0) 
PtO 
Pt  5 
Pt  6 


Center  of  target  on  chair  back 
Center  of  upper  left  target  on  calibration  device 
Center  of  upper  right  target  on  calibration  device 
Center  of  lower  right  target  on  calibration  device 


DATA 


Origin  (0,0) 
Pt  1 
Pt  2 
Pt  3 
Pt  4 


Same  as  calibration 
Center  of  left  Tj  mount  target 
Center  of  right  T  j  mount  target 
Center  of  head  mount  target 
Center  of  mouth  mount  target 
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X-RAY  ANTHROPOMETRY  DATA  PROCESSING  AND  INITIAL  CONDITIONS 


Coordinate  measurements  were  performed  at  Eglin  AFB  on  the  same  Benson- 
Lehner  Telereadex  machines  which  were  utilized  for  the  photographic  data 
reduction.  (See  Photographic  Data  Processing)  The  X  and  Z  grid  of  the  Telereadex 
device  established  the  coordinate  reference  about  the  arbitrary  zero  point  or  origin. 
The  card  format  for  digitized  x-ray  data  was  identical  to  the  digitized  photographic 
data  with  the  addition  of  the  x-ray  number  in  card  columns  69  through  74,  as  shown 
in  Table  1 .  Table  3  provides  x-ray  data  point  identification  for  the  skull  and  Table  4 
similarly  for  the  spine. 

There  is  one  difference  in  calibration  between  x-ray  and  photographic  data 
that  arises  from  the  fact  that  a  lead  calibration  rule  was  used  for  the  x-ray  instead 
of  an  L-shaped  calibration  device.  Slight  differences  between  the  X  and  Z  of  rhe 
reader  machine  could  be  a  source  of  error  if  the  rule  were  read  only  once.  Therefore, 
the  two  points  of  the  rule  were  read  twice;  once  with  the  rule  image  aligned  sub¬ 
stantially  along  the  Z  axis  of  the  Telereadex  and  once  with  the  rule  substantially 
along  the  X  axis. 

In  some  cases  the  x-ray  was  not  a  "true"  lateral  or  there  may  have  been  subject 
asymmetry.  A  best  estimate  of  the  Frankfort  plane  line  was  derived  from  reading 
the  two  external  auditory  meatus  point  marker  locations  and  the  two  infraorbital 
notch  point  marker  locations.  The  average  of  these  readings  was  the  best  estimate 
of  a  line  in  the  Frankfort  plane.  In  those  cases  where  the  points  overlapped,  the 
readings  were  made  twice  to  maintain  consistency. 

All  x-rays  and  x-ray  cards  were  and  are  maintained  in  files  as  a  primary 
reference. 
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TABLE  3 

X-RAY  DATA  POINT  IDENTIFICATION 


SKULL  X-RAY 

Point  Identifier  Number  Description 


01 

02 

03 

04 


Z  reading  of  one  calibration  point  (skull  data) 

Z  reading  of  second  calibration  point  (skull  data) 
X  reading  of  one  calibration  point  (skull  data) 

X  reading  of  secsnd  calibration  point  (skull  data) 


05 

06 

07 

08 


External  auditory  meatus  superior  edge  center 
marker 

External  auditory  meatus  superior  edge  center 
marker 

Infraorbital  notch  marker 
Infraorbital  notch  marker 


09 

10 

11 

12 

13 

14 


Head  mount  superior  anterior  screw  hole  point 
Head  t  Hint  superior  middle  screw  hole  point 
Head  t.iount  superior  posterior  screw  hole  point 
Head  mount  interior  anterior  screw  hole  point 
Head  mount  inferior  middle  screw  hole  point 
Head  mount  inferior  posterior  screw  hole  point 


15 

16 

17 

18 

19 

20 


Mouth  mount  superior  anterior  screw  hole  point 
Mouth  mount  superior  middle  screw  hole  point 
Mouth  mount  superior  posterior  screw  hole  point 
Mouth  mount  inferior  anterior  screw  hole  point 
Mouth  mount  inferior  middle  screw  hole  point 
Mouth  mount  inferior  posterior  screw  hole  point 
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TABLE  4 


CERVICAL  SPINE  X-RAY 


Point  Identifier  Number 

51 

52 

53 

54 

55 

56 

57 


61 

62 

63 

64 

65 

66 


Description 

Z  reading  of  one  calibration  point  (spine  data) 

Z  reading  of  second  calibration  point  (spine  data) 
X  reading  of  one  calibration  point  (spine  data) 

X  reading  *>f  second  calibration  point  (spine  data) 

Superior  anterior  corner  of  the  body  T] 

Superior  posterior  corner  of  the  posterior  spinous 
process  of  T  j 

Inferior  posterior  corner  of  the  posterior  spinous 
process  of  T  j 

T]  mount,  superior  anterior  screw  hole  point 
T]  mount,  superior  middle  screw  hole  point 
T]  mount,  superior  posterior  screw  hole  point 
fj  mount,  inferior  anterior  screw  hole  point 
Tj  mount,  inferior  middle  screw  hole  point 
Tj  mount,  inferior  posterior  screw  hole  point 
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ANXRAY  -  The  X-ray  Anthropometry  Program 

ANXRAY  is  a  FORTRAN  program  which  executes  the  computations  necessary 
to  locate  and  orient  the  six  accelerometers  used  in  the  two  dimensional  experiments 
within  their  respective  anatomical  coordinate  systems. 

Input  to  ANXRAY  consists  of  computer  cards  containing  digitized  x-ray  data. 
ANXRAY  used  this  input  in  the  following  procedure: 


1 .  The  digitized  data  cards  are  read  and  X,  Z  coordinates  are  obtained  for 
each  point. 

2.  Points  01  through  04  are  used  to  obtain  scaling  information  for  points 
05  through  20.  (See  XRAY  Point  Identification  Table) 

|"(Z3-Z4)s -(Zl -Z2)s 


cx  = 


cz  = 


(Z3  -  Z4)2  (a'1  -  X2)s  -  (X3  -  X4)a  (Z1  -  Z2)a 
» 

(X3  -  X4)a  -  (XI  -  X2)a 


.  (15.24)  (1) 


1% 


(X3  -  X4)s  (Z1  -  Z2)s  -  (XI  -  X2)s  (Z3  -  Z4)s 


.  05.24)  (2) 


where  15.24  cm  -  6  in  and  XK,  ZK  are  the  coordinates  of  point  K. 

3.  The  values  of  the  remaining  16  points  are  scaled  by  these  conversion 

factors: 


XK  =  (XK).(CX)  (3) 

ZK  =  (ZK).(CZ)  (4) 

4.  The  position  and  orientation  of  the  head  anatomical  coordinate  system 
relative  to  an  arbitrary  reference  frame  established  by  the  Telereadex  reduction 
system  are  determined. 

AH  =  (X5  +  X6)/2  (5) 

BH  =  (Z5+Z6)/2  (6) 

0H  =  ARCTAN  fZ5lZ6^Z7-Ze| 

1X5  +  X6  -  X7  -  X8j  w 


Careful  consideration  was  given  to  appropriate  quadrant  assignment  for  the 
related  angle  which  was  obtained  by  the  FORTRAN  ARCTAN  routine  in  all  equations 
of  all  programs  employing  it. 

5.  Points  09  through  14  are  used  to  compute  the  position  and  orientation  of 
the  head  mount  relative  to  the  arbitrary  reference  frame. 

AD  =  (X9+X10  +  X11  f  X12  +  X13  +  X14)/6 
BD  =  (Z9+Z10  +  Z11 +Z12+Z13  +  Z14)/6 
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(8) 

(9) 


For  the  orientation  of  the  head  mount,  the  coordinates  far  the  superior  screw  holes 
are  first  averaged  with  the  coordinates  of  the  corresponding  inferior  screw  holes: 


XXI  =  (X9  +  X12)/2 

00) 

ZZ1  =  (Z9  +  Z 1 2)/2 

(11) 

XX2  =  (X10  +  X13)/2 

(12) 

ZZ2  =  (Z10  +  Z13)/2 

(13) 

XX3  =  (XI 1  +XU)/2 

(14) 

ZZ3  =  (Zll  +Z14)/2 

(15) 

This  yields  the  coordinates  for  three  points  on  the  midline  of  the  mount.  A  least 
squares  technique  is  then  used  to  calculate  the  slope  of  the  best  fit  line  using  these 
points: 

0  D  =  ARCTAN 


6.  Points  J5  through  20  are  used  in  a  similar  manner  to  calculate  the  position 
and  orientation  of  the  mouth  mount  relative  to  the  arbitrary  reference  frame:  AM, 

BM  and  0M. 

7.  The  position  and  orientation  of  the  head  mount  in  the  head  anatomical 


coordinate  system  are  then  calculated: 

XHD  =  (AD  -  AH)  COS  (0H)  -  (BD  -  BH)  SIN  (0H)  (17) 

ZHD  =  (AD  -  AH)  SI  N  (0H)  +  (BD  -  BH)  COS  (0H)  (18) 

0HD  =  0D  -  0H  (19) 

8.  Similarly,  the  position  and  orientation  of  the  mouth  mount  in  the  head 
anatomical  coordinate  system  become: 

XHM  =  (AM  -  AH)  COS  (0H)  -  (BM  -  BH)  SIN  (0H)  (20) 

ZHM  =  (AM  -  AH)  SIN  (0H)  +  (BM  -  BH)  COS  (0H)  (21) 

0HM  =  0M-0H  (22) 


9.  It  is  now  possible  to  utilize  a  constant  taken  from  the  transducer  mount 
drawings  and  the  previously  calculated  location  and  orientation  of  the  head  and 
mouth  mounts  to  calculate  the  position  and  orientation  of  the  individual  accelero¬ 
meters  in  the  head  anatomical  coordinate  system. 

LI  =  1.676654  cm 
XAFM  =  XHM  +%  LI  COS  (0HM) 

ZAFM  =  ZHM-%  LI  SIN  (0HM) 

OAFM  =  0HM 


;XX)  fcZZ)  -3&XXZZ)' 
£XX)s  -  3CCXXd)  * 
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(23) 

(24) 

(25) 


XASM  =  XHM  -  %  LI  COS  (0HM) 

(26) 

ZASM  =  ZHM  LI  SIN  (0HM) 

(27) 

OASM  =  0HM 

(28) 

XAFD  =  XHD  +  h  LI  COS  (0HD) 

(29) 

ZAFD  =  ZHD  -  h  LI  SIN  (0HD) 

(30) 

OAFD  =  0HD 

(31) 

XASD  =  XHD  -  h  LI  COS  (0HD) 

(32) 

ZASD  =  ZHD  +  %  LI  SIN  (0HD) 

(33) 

OASD  =  0HD 

(34) 

10.  The  head  span  is  calculated: 

HSP  =  f  (AM  -  AD)2  +  (BM-BD)V 

(35) 

11.  The  skull  x-ray  processing  is  completed  with  the  calculation  of  the 
head  span  line  orientation  relative  to  the  +Z  axis  of  the  head  anatomical  coordinate 
system.  From  0H,  the  orientation  of  the  head  anatomical  coordinate  system  +X  axis 
relative  to  the  x-ray  +X  axis,  the  orientation  of  the  head  anatomical  +Z  axis  is 
calculated: 

0HZ  =  0H  +  1T/2  (36) 

where  the  sign  is  dependent  on  the  sign  of  0H .  Then,  the  orientation  of  the  head 
span  line  relative  to  the  x-ray  +X  axis  is: 

AHS  =  ARCTAN  [(BM  -  BD)/(AM  -  AD)]  (37) 

Finally,  the  desired  orientation  is: 

AHS  =  AHS  -  0HZ  (38) 

12.  For  spine  x-ray  anthropometry,  points  51  through  54  are  used  to  obtain 
scaling  information  for  points  55  through  57  and  61  through  66.  The  procedure  is 
exactly  the  same  as  that  used  to  obtain  scaling  information  for  the  skull  x-ray  . 

13.  The  coordinates  of  points  55  through  57  and  points  61  through  66  are 
scaled  by  these  conversion  factors.  The  procedure  is  the  same  as  for  the  skull  x-ray  . 

14.  The  position  and  orientation  of  the  T^  anatomical  coordinate  system 
relative  to  the  arbitrary  reference  frame  are  calculated: 


AN  =  X55 
BN  =Z55 
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(39) 

(40) 


(41) 


r  BN  -  [  (Z56  +  Z57)/2] ' 
0N  =  7T  +  ARCTAN  |_AN  _  [  (X56  +  X57)/2] . 


15.  Points  61  through  66  are  used  to  compute  the  position  and  orientation 
(AT,  BT  and  0T)  of  the  T|  mount  relative  to  the  arbitrary  reference  frame  in  the 
same  manner  described  for  the  head  mount. 


16.  The  position  and  orientation  of  the  Tj  mount  relative  to  the  Ti 
anatomical  coordinate  system  are  computed: 

XNT  =  (AT  -  AN)  COS  (0N)  -  (BT  -  BN)  SIN  (0N)  (42) 

ZNT  =  (AT  -  AN)  SIN (0N)  +  (BT  -  BN)  COS(0N)  (43) 

0NT  =  0T-0N  (44) 


17.  Finally,  the  position  and  orientation  of  the  Tj  accelerometers  are: 


XAFN  =  XNT  +  %  LI  COS(0NT)  (45) 

ZAFN  -ZNT  -%  LI  SIN(0NT)  (46) 

OAFN  =  0NT  (47) 

XASN  «  XNT  -  h  LI  COS (0 NT)  (48) 

ZASN  =  ZNT  +%  LI  SIN  (0NT)  (49) 

OASN  =  0NT  (50) 


The  output  from  ANXRAY  consists  of  punch  cards  containing  the  values  computed. 
These  cards  were  input  to  the  ICONDS  program;  ICONDS  added  the  ANXRAY 
computed  values  to  the  ICONDS  output  file. 

ANXRAY  Validation 

ANXRAY  validation  was  accomplished  as  follows: 

1.  An  examination  was  made  of  the  horizontal  and  vertical  scaling  constants. 
Any  substantial  difference  would  indicate  improper  calibration  of  the  photographic 
data  system. 

2.  An  examination  was  made  of  the  head  span  (XHSP)  computed  from  the 
skull  x-ray.  This  variable  was  one  of  the  first  and  best  indications  of  incorrect 
data  reduction  and  could  readily  be  compared  withHSP  from  Program  TROUTP. 

3.  The  scaled  coordinates  of  each  point  were  plotted  to  provide  a  quick 
indication  of  data  reduction. 


4.  The  scaled  coordinates  were  then  backscaled  to  the  scale  of  the  x-ray 
and  plotted.  These  plots  were  overlaid  on  their  corresponding  original  x-rays.  In 
all  cases,  the  comparison  was  very  good. 
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5.  Finally,  a  straight  edge  was  used  to  approximate  the  anatomical 
coordinate  systems  on  the  scaled  coordinates  plots.  A  protractor  and  ruler  were  then 
used  to  measure  directly  the  positions  and  orientations  which  were  computed  in 
ANXRAY.  There  was  no  essential  difference  noted  between  the  measured  values 
and  the  computed  values.  It  Is  not  considered  either  feasible  or  necessary  to  further 
validate  these  data* 


ICONDS  -  The  Initial  Conditions  Program 


I  ICONDS  is  a  FORTRAN  program  which  calculates  the  horizontal  and  vertical 

scaling  constants  for  reduced  photographic  data  and  establishes  the  initial  positions 
of  the  subject  with  respect  to  the  Telereadex  reference  coordinate  system;  and  the 
initial  orientations  of  the  subject  head  and  neck  with  respect  to  the  laboratory 
reference  coordinate  system. 

ICONDS  utilizes  three  sources  of  input  data: 

1 .  Average  of  pre-acceleration  transducer  data  for  each  of  the  head,  mouth 
and  T]  mounted  accelerometers. 

j 

|  2,  Punch  cards  containing  the  reduced  photographic  target  coordinates  for 

i<  ten  frames  of  calibration  and  ten  frames  of  pre-acceleration  photographic  data.  For 
the  calibration,  targets  were  located  on  an  inverted  L-shaped  device  of  known  size 
;  with  arms  perpendicular  and  parallel  to  gravity  as  shown  in  Figures  DP- 10  and  DP-6. 

For  the  data,  two  targets  were  located  on  the  T  |  mount  and  one  each  on  the  head 
i  mount,  the  mouth  mount  and  the  sled  as  shown  in  Figure  DP-9. 

3.  X-ray  anthropometric  data  consisting  of  the  coordinates  and  orientation 
of  the  head  and  mouth  accelerometers  relative  to  the  head  anatomical  coordinate 
system;  and  the  coordinates  and  orientation  .of  the  Tj  accelerometers  relctive  to  the 
j  T]  anatomical  coordinate  system.  (See  Program  ANXRAY). 

I 

s  The  first  step  in  the  execution  of  ICONDS  is  the  calculation  of  average  coordinates 

fj  from  the  10  frames  of  calibration  for  the  upper  left  target,  point  0;  for  the  upper 

■  right  target,  point  5;  and  for  the  lower  right  target,  point  6,  as  shown  in  Figure 

S  DP-10,  Similarly,  average  coordinates  from  the  10  frames  of  data  are  calculated  for 

j  the  left  Tl  target,  point  1;  the  right  T  i  taiget,  point  2;  the  head  target,  point  3;  and 

|  the  mouth  target,  point  4,  as  shown  in  Figure  DP-9. 

The  average  coordinates  from  the  calibration  frames  are  used  to  convert  the 
remaining  reduced  data  from  Telereadex  units  to  centimeters  and  to  correct  for  possible 
angular  difference  between  the  laboratory  reference  coordinate  system  and  Telereadex 
reference  coordinate  system.  As  a  first  step  in  the  conversion  process,  ICONDS 
calculates  two  angles.  ALPHmX  and  ALPHAZ  as  shown  in  Figure  DP-  10! 

|  ALPHAX  =  ARCTAN  [  (Z5  -  Z0)/(X5  -  X0)]  (1) 

j  ALPHAZ  =  ARCTAN  [  (X6  -  X5)/(Z5  -  Z6)j  (2) 

j  The  small  magnitude  of  these  two  angles  permit  averaging  to  obtain 

’  ALPHA  =  (ALPHAX  +  ALPHAZ)/2  (3) 

This  angle  ALPHA  is  defined  as  the  angular  difference  between  the  Telereadex 
reference  coordinate  system  and  the  laboratory  reference  coordinate  system.  An 
angle  Beta  in  the  Telereadex  reference  coordinate  system  would  be  (BETA  -  ALPHA) 
in  the  laboratory  reference  coordinate  system.  Due  to  possible  variation  in  the 
horizontal  and  vertical  measurement  of  the  photographic  reduction  system,  two 
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FIGURE  DP-  10  Illustration  of  the  Angular  Quantities  used  to  Relate  T< 
Coordinate  System  Z  Axis  to  Laboratory  Reference  Coordinate 
(Gravity).  ALPHAX  and  ALPHAZ  are  exaggerated  for  illust 


conversion  factors  are  calculated.  Since  the  distances  from  point  0  to  point  5 
and  from  point  5  to  point  6  are  both  one  foot, 

1  foot  -  (X5  -  XO)  COS  (ALPHA)  horizontal  Telereadex  units  (4) 

1  foot  =  (Z5  -  76)  COS  (ALPHA)  vertical  Telereadex  units  (5) 

The  average  coordinates  of  the  targets  shown  in  Figure  DP-9  are  divided  by  these 
constants  to  convert  them  to  feet  and  then  multiplied  by  30.48  to  convert  them  to 
centimeters. 

Specific  discussion  of  parameter  derivation  in  ICONDS  will  be  limited  to  those 
parameters  which  are  utilized  in  the  calculation  of  the  output  variables  presented. 

ADH  is  the  angle  formed  by  the  line  connecting  the  center  of  the  head  target  and 
the  center  of  the  mouth  target  with  the  4X  ax  is  of  the  laboratory  reference  coordinate 
system,  as  shown  in  Figure  DP-11. 


ILLUSTRATION  OF  VARIABLES  ADH,  RN, 
AND  POINTS  X34,  Z34,XN,ZN 


FIGURE  DP- 11 

ADH  =  ARCTAN  [(Z3  -  Z4)/(X3  -  X4)]  -  ALPHA  (6) 

X34,  Z34  are  the  coordinates  of  the  midpoint  of  the  line  joining  the  head  and 
mouth  targets  with  respect  to  the  Telereadex  reference  coordinate  system. 
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X34  =  (X3  +  X4)/2  (7) 

Z34  =  (Z3  +  Z4)/2  (8) 

XN,  ZN  are  the  coordinates  of  rtte  midpoint  of  the  line  joining  the  two  T. 
targets  with  respect  to  the  Tetereadex  reference  coordinate  system* 

XN  =  (XI  +  X2)/2  (9) 

ZN  =  (Zl+Z2)/2  (10) 


It  should  be  noted  that  X34,  Z34  and  XN/  ZN  are  corrected  to  the  laboratory 
reference  coordinate  system  prior  to  their  boing  utilized  as  initial  conditions.  (See 
Program  TROUTP) 

RN  is  the  angle  formed  by  the  line  connecting  the  two  T  j  targets  with  the  +X  axis 
of  the  laboratory  reference  coordinate  system. 

RN  =  ARCTAN  [(Z2  -  Z1)/(X2  -  XI)]  -  ALPHA  (11) 

ICONDS  is  a  multi-run  program  processing  any  number  of  runs  without  reloading 
the  program.  In  order  to  expedite  tne  processing  of  ICONDS,  the  input  cards  are 
read  from  a  card  image  digital  tape  which  contains  the  ICONDS  input  data  for  all 
of  the  runs.  The  output  from  ICwNDS  was  a  digital  tape  containing  60  initial 
condition  parameters  which  were  either  derived  by  ICONDS  or  provided  by  ANXRAY. 

ICONDS  Validation 

ICONDS  validation  involved  three  different  objectives: 

1.  Validation  of  the  horizontal  and  vertical  scaling  constants  developed  by 
ICONDS  for  scaling  photographic  data. 

2.  Validation  of  the  initial  conditions  which  were  computed  from  pre¬ 
acceleration  photographic  data. 

3.  Validation  of  the  x-ray  anthropometric  data. 

The  best  test  for  evaluating  the  horizontal  and  vertical  scaling  constants  is  the 
examination  of  scaled  photographic  data  plots,  particularly  when  those  plots  compare 
the  scaled  photographic  data  with  independently  scaled  sensor  data.  An  examination 
of  the  PDH  (derived  from  the  head  angular  position  as  computed  from  photographic 
data)  vs  RDH  (derived  from  the  scaled  rate  gyro)  comparison  plots  633  through  673 
illustrates  the  precision  and  accuracy  of  photographic  and  sensor  scaling. 

The  initial  conditions  which  were  computed  from  pre-acceleration  photographic 
data  include  the  scaled  initial  coordinates  of  the  photographic  targets  and  initial 
condition  values  which  were  verified  in  the  validation  of  scaling  constants.  The 
values  computed  from  these  coordinates  involve  very  basic  equations  (midpoint  of  a 
line  joining  two  targets,  slope  of  a  line  connecting  two  targets).  Those  computations 
were  carefully  checked  for  accuracy. 
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ufed  in  the  program  ANXRAY 
hese  values  were  input  to 

The  ICONDS  output  tape  was  and  is  maintained  as  a  primary  reference* 

, ,  Jl'9Mre  DP- 12  illustrates  the  flow  of  photographic  and  x-ray  data  to  a  UNIVAC 
j  108  output  tape* 


Finally,  the  x-ray  anthropometric  values  were  cor 
(See  Program  ANXRAY)  and  were  validated  therein. 
ICONDS,and  ICONDS  only  added  them  to  firs  output 
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FIGURE 


The  scaled  data,  photographic  and  sensor,  were  then  subfected  to  output 
variables  and  transformation  programs.  Initial  conditions  and  x-ray  anthropometry 
were  required  by  these  programs. 
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SENSOR  OUTPUT  VARIABLES  AND  SUBSEQUENT  TRANSFORMATIONS 
OUTVAR  -  The  Output  Variables  Program 


OUTVAR  Is  a  FORTRAN  program  which  takes  scaled  accelerometer  and  rate 
gyro  data  from  an  experimental  run  and  combines  this  with  the  appropriate  initial 
conditions  data  to  compute  various  kinematic  descriptors  of  subject  motion. 

Input  to  OUTVAR  is  the  digital  output  tape  from  the  DSCALE  program,  the  digital 
output  tape  from  the  ICONDS  program  and  a  set  of  data  cards  which  include  run 
numbers,  run  processing  options,  and  output  medium  choices. 


Typically,  a  set  of  scaled  input  data  (produced  by  program  DSCALE)  consists 
of  the  variables  AFM,  ASM,  AFD,  ASD,  AFN,  ASN,  RVH,  RVN  and  AXC  where 
each  of  these  three-character  "labels"  refers  to  a  time  series  of  approximately  1000 
sensor  values  spaced  0.5  msec  apart. 

Since  some  experimental  runs  may  not  possess  this  set  in  total,  OUTVAR  will  not 
attempt  to  compute  any  quantity  which  is  a  function  of  a  missing  variable.  This 
selectivity  of  processing  is  implemented  by  subjugating  the  program's  control  to  a 
run  option  number  which  is  one  of  several  parameters  used  to  govern  the  program 
execution.  The  four  options  available  to  OUTVAR  have  the  following  meanings: 

Option  1  ....  process  the  run  as  having  all  variables. 

Option  2  ....AFN,  ASN,  RVN  and  AXC  are  the  only  variables 
available,  that  is,  there  are  no  transducers  at  the  head  or  mouth. 


meters. 

accelerometers. 


Option  3  . .  ..all  variables  are  present  except  the  head  accelero- 
Option  4,,  ..all  variables  are  present  except  the  mouth 


For  each  experimental  run  the  set  of  scaled  accelerometer  and  rate  gyro  data  has 
an  associated  set  of  "initial  conditions"  information  (see  program  ICONDS).  In 
performing  its  functions,  OUTVAR  requires  selected  portions  of  this  data. 

When  all  the  control  parameters,  scaled  run  data,  and  Initial  conditions  infor¬ 
mation  are  read  in,  OUTVAR  begins  processing  a  run.  The  following  discussion  will 
be  limited  to  those  variables  computed  by  OUTVAR  which  are  presented  in  the  plots 
and/or  variables  required  in  the  computation  of  presented  variables.  The  OUTVAR 
list  in  the  Variable  Definitions  section  gives  precise  variable  definitions  and  cross- 
reference  to  appropriate  figures. 


OUTVAR  calculates  the  angular  acceleration  of  the  head,  RAH,  and  the  angular 
acceleration  of  the  spine,  RANI,  by  smoothing  and  differentiating  the  head  and  the 
spine  angular  velocities  (RVH  and  RVN  respectively)  measured  by  the  rate  gyros. 
The  smoothing  consists  of  an  eleven-point  central  moving  average  on  RVH  and  RVN 
and  is  followed  by  a  three-point  central  differentiation  formula. 
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The  eleven-point  central  moving  average  on  the  set  )((1)  through  X(1002) 
applied  as  follows: 


/  (2i  -  1)  for  i  =  1  to  5  (1) 


/II  for  i  =6  to  997  (2) 


/  (2005  -  2i )  for  i  =  998  to  1002  (3) 


The  three-point  central  differentiation  formula  for  the  ith  point  is: 

Xi+1  X,.!  X;-Xf-1  Xj»  ]  «  Xj 

XDOTj  =  - = - + - (4) 

2h  2h  2h 

where  h  is  the  interval  step  size  (0.5  msec). 

Inserting  the  values  of  RVH  and  RVN  in  place  of  the  X  values  in  the  algorithms 
above  results  in  RAH  and  RAN  respectively,  as  represented  by  equation  (4). 

RAH  as  used  in  the  co-location  corrections  was  obtained  by  this  method.  The 
RAH  which  is  presented  in  the  plots  was  obtained  by  filtering  RVH  with  a  second- 
order  numerical  filter  with  a  damping  ration  (£)  of  0,6  and  a  natural  frequency 
(^n/2ir)  of  50  Hertz.  The  filtered  RVH  was  then  differentiated  in  accordance  with 
equation  (4). 

RPN  and  RPH,  angular  displacements  of  the  spine  and  head,  are  computed  by 
integrating  RVN  and  RVH.  In  both  instances  the  trapezoidal  rule  is  employed  and 
the  ICONDS  angles  ADH  and  RN  are  used  as  initial  conditions. 
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The  trapezoidal  rule  algorithm  on  the  set  X(l)  through  X(1002)  is  applied  as  follows 
Yj  =  PQ  (where  PQ  is  the  initial  position) 

Y2  =  Y,  +  (X,  +X2)  •  (h/2)  (5) 


Yi  -  Yi  - 1 


<xi- 


+  x. ) 


(h/2) 


(6) 


Y1002  =  Y1001  +  ^X1001  +  X1002^  *  ^ 

where  h  is  the  integration  step  size  (0.5  msec). 

OUTVAR  corrects  the  head,  mouth  and  spine  acceleration  for  non-colocation  of 
the  anterior  and  superior  accelerometers  within  each  mount  system.  Theoretically, 
it  is  desired  that  the  anterior  and  superior  accelerations  be  a  measure  of  orthogonal 
acceleration  components  at  one  specific  point  on  each  mount;  however,  due  to  the 
physical  impossibility  of  mounting  two  accelerometers  at  the  same  point  they  are 
mounted  so  that  they  measure  orthogonal  accelerations  at  two  distinct  points  which 
are  equidistant  from  and  colinear  with  the  throretically  desired  point  (Figure  DP-13). 


D  =  0  9019  9 


FIGURE  DP-13  Geometry  of  Accelerometer  Pair  Within  the  Mouth  Mount 
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By  utilizing  the  angular  acceleration  and  angular  velocity  of  both  the  head  and 
spine,  a  pair  of  corrected  accelerations  are  computed  for  each  mount  which 
represent  the  desired  orthogonal  accelerations  at  the  midpoint  of  the  line  |oining 
the  centers  of  the  accelerometers. 

Given  the  configuration  of  accelerometers  as  shown  in  Figure  DP-13,  where  the 
anterior  acceleration  is  recorded  for  point  Po>  and  the  superior  for  point  Pj,  the 
anterior  and  superior  accelerations  at  point  r  are: 


AFMp 

=  afmP9 

+  D.co3 

(8) 

ASMp 

=  asmp1 

+  D.a 

(9) 

where  oU“  the  angular  velocity  of  the  head,  RVH,  0=  the  angular  acceleration 
of  the  head,  RAH,  and  D  is  the  distance  from  the  midpoint  of  the  accelerometer 
to  the  midpoint  of  the  line  connecting  the  centers  of  the  two  accelerometers. 

Similar  corrections  are  made  for  the  head  and  spine  mounts;  except  that  in  the  case 
of  the  spine  corrections,  to  is  RVN  and  a  is  RAN.  These  "corrected"  accelerations 
are  assigned  the  same  three-character  names  as  were  the  input  accelerations  and 
subsequent  program  computational  requirements  referencing  the  individual  mount 
accelerations  are  furnished  with  the  corrected  accelerations. 

By  two  successive  applications  of  the  trapezoidal  algorithm  to  the  sled  accelera¬ 
tion  data  AXC,  with  initial  velocity  and  position  zero,  we  derive  the  sled  velocity 
XVC  and  position  XDC  relative  to  the  laboratory  reference  coordinate  system. 

As  output,  OUTVAR  is  able  to  deliver  some  seventy  variables  representing 
angular  and  linear  acceleration,  velocity  and  displacement  of  head  and  spine 
instrumentation  coordinate  systems, 

OUTVAR  is  a  multi-run  program  enabling  these  procedures  to  be  executed  for 
thirty  runs  without  reloading  the  program. 

OUTVAR  Validation 

To  validate  the  computing  software  of  OUTVAR  for  algorithm,  programming  and 
computational  errors,  simulated  data  were  generated. 

The  input  data  required  by  OUTVAR  consists  of  DSCALE  and  ICONDS  data,  both 
on  digital  tape. 

DSCALE  data  designed  to  simulate  a  fully  instrumented  run  should  consist  of  six 
subject  accelerdtion  data  streams  (AFM,  ASM,  AFD,  ASD,  AFN,  ASN),  two  subject 
angular  velocity  data  streams  (RVH  and  RVN)  and  one  sled  acceleration  data  stream 
(AXC). 

ICONDS  simulation  data  must  consist  of  60  quantities  describing  the  simulated 
initial  conditions  and  including  simulated  x-ray  anthropometric  data. 

A  model  will  be  described  which  was  the  basis  for  generating  the  required 
simulated  test  data. 
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Consider  an  apparatus  of  five  rods  (Lx  ,  L  ,  L  ,  L4  ,  L  )  connected  as  shown 
Jn  Figure  DP-14.  In  the  list  of  definitions  which  follow,  all  fixed  numerical 
quantities  are  chosen  to  bear  a  close  resemblance  to  typical  experimental  values 
and  to  keep  calculations  as  simple  as  possible. 


Lg  ^ 


FIGURE  DP-14  Illustration  of  Simulation  Model 

1 .  Lx  will  act  as  the  line  joining  the  center  of  the  head  mount  to  the 
center  of  the  mouth  mount.  Let  its  length  be  30  centimeters,  a  number  resembling 
typical  headspans. 

2.  Ls  will  represent  the  line  which  joins  the  two  corner-mounted  photo¬ 
targets  of  the  spine  mount, 

3.  Lg  will  be  the  line  joining  the  center  of  Lg  to  the  center  of  Lt  .  Its 
value  shall  be  20  centimeters,  a  value  close  to  those  computed  from  initial  conditions 
data. 


4.  L^  and  will  represent  the  lines  joining  the  centers  of  the  anterior 
and  superior  accelerometers  at  the  mouth  and  head  mounts  respectively.  They  will 
cross  i1  at  their  centers  and  in  both  cases  it  will  be  assumed  that  the  sensitive  axes 
of  the  anterior  and  superior  accelerometers  are  respectively  colinear  with  and  perpen¬ 
dicular  to  these  lines.  The  length  of  L^  is  equal  to  the  length  of  LB  which  is  equal 
to  .90199  cm,  the  value  used  in  program  OUTVAR. 

It  would  also  be  possible  to  add  a  sixth  rod  L*  whose  center  would  coincide  with 
the  juction  of  Lg  and  Lg  and  which  would  have  the  same  length  as  and  Lg  and 
would  represent  the  line  joining  the  centers  of  the  spinal  accelerometers.  Since 
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this  rod  would  not  move  relative  to  Lfl  it  is  only  necessary  to  keep  its  relative 
position  in  mind  as  seen  in  Figure  DP-15, 


FIGURE  DP-15  Illustration  of  Spine  Mount  Simulation  Data 

Programs  OUTVAR  and  TROUT  do  not  see  the  complex  linkages  uniting  experi¬ 
mental  components  so  that  separate  simplified  driving  functions  may  be  applied 
to  produce  linear  acceleration  and  angular  acceleration. 
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Let  the  input  sled  acceleration  be  given  by. 


a(t) 

=  (-250g/sec) 

*(t) 

for  0  st  s4Q  msec 

(10) 

a(t) 

~  (250g/$ec)  • 

(t)  -20g 

for  40  msec  st  s80  msec 

(ID 

c(t) 

=  0 

for  80  msec  sts  500  msec 

(12) 

That  is,  the  sled  acceleration  waveform  will  be  an  isosceles  triangle  with  a  base 
of  80  milliseconds  and  height  of  -  iOg's  ,  representing  a  nominal  iOg  run  at  onset 
of  250 g's  per  second. 

It  will  be  assumed  that  this  sled  acceleration  is  passed  unchanged  to  the  generator 
apparatus  and  is  applied  in  such  a  way  that  only  translation  of  the  apparatus  is 
effected  and  that  there  is  no  rotation  due  to  its  action. 

A  controlled  rotational  effect  will  be  superimposed  on  the  apparatus  as  follows: 

At  the  end  of  the  experiment,  500  msec  after  its  start,  the  orientation  of  the  apparatus 
is  to  be  -70  degrees  relative  to  the  laboratory  reference  coordinate  system  X  axis, 
having  sustained  a  rotation  of  -50  degrees  about  the  spine  point  PQ  (Figure  DP-15!) 
while  not  having  sustained  any  internal  relative  rotations.  The  average  angular 
velocity  of  both  L  and  L  is  therefore  -  100  degrees/sec.  In  order  that  a  non-zero 
angular  acceleration  be  involved,  let  the  spine  and  head  rods  La  and  L^^  sustain  an 
angular  velocity  given  by 

a)(t)  =  (800  degrees/sec2)  •  (t)  Osi  s  250  msec  (13) 

co(t)  =  (-800  degnees/sec3)  •  (t)  +  400  degrees/sec 

250  msec  sit  s;500  msec  (14) 

This  results  in  the  required  average  angular  velocity  of  -  100  degrees/sec. 

Figures  DP-16  and  DP-17  depict  both  the  linear  and  angular  driving  functions. 

A  program  was  written  to  derive  the  appropriate  data  streams  for  both  DSC  ALE 
and  ICONDS  simulated  data  when  the  data  generator  apparatus  is  driven  by  these 
functions  with  the  initial  conditions. 

Once  the  simulated  ICONDS  and  DSCALE  data  had  been  generated  and  recorded 
on  digital  tape  in  the  proper  formats,  program  OUTVAR  was  run  using  these  tapes 
as  input.  Correction  tor  non-colocation  was  specified.  The  critical  test  for 
OUTVAR's  accuracy  is  in  the  output  variables  XDH  and  ZDH,  the  X  and  Z  coordinates  . 
of  the  displacement  of  the  center  of  tha  head  relative  to  the  laboratory  reference 
coordinate  system.  If  these  values  are  correct  then  OUTVAR  is  performing  properly. 

The  predicted  values  for  these  quantities  after  500  msec  into  the  "model" 
experiment  are 

X  -  -161.638  cm 

Z  =  6.840  cm 
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The  values  computed  by  OUTVAR  are 
X  -  -161.616  cm 
Z  =  6.813  cm 

which  gives  rise  to  an  absolute  difference  of 

|ax|  =  0.022  cm 

|az|  =  0.027  cm 

The  source  of  these  small  errors  is  most  likely: 

1.  The  colocation  correction  of  simulated  accelerations  is  computed  using 
angular  accelerations  derived  from  the  simulated  accelerations  prior  to  colocation 
correction.  This  will  lead  to  small  computational  errors. 

2.  Small  errors  are  contributed  by  round  off. 

Evaluation  of  the  simulated  OUTVAR  output  data  demonstrated  that  the  difference 
in  computed  values  and  predicted  values  could  be  considered  insignificant. 

The  output  from  OUTVAR  was  a  digital  tape  containing  approximately  70  variables 
for  up  to  30  runs.  These  tapes  were  and  are  maintained  as  a  primary  reference. 

When  the  OUTVAR  processing  was  completed,  x-ray  anthropometric  data  and 
center  of  gravity  data  were  used  to  transform  data  to  anatomical  coordinate  systems 

and  to  the  head  C.G. 


TROUT  -  The  Sensor  Data  Transformation  Program 

TROUT  is  a  FORTRAN  program  which  accepts  input  data  from  OUTVAR,  DSCALE 
and  ICONDS,  and  converts  accelerations  to  head  anatomical  coordinates  utilizing 
ANXRAY  data  to  locate  anatomical  reference  points.  It  will  also,  on  request, 
transform  the  data  to  the  head  principal  coordinate  system. 

Input  to  TROUT  consists  of: 

1,  Co-located  acceleration  variables  AFM,  ASM,  AFD,  ASD,  AFN,  ASN; 
measured  angular  velocities  of  the  head  and  neck  RVH  and  RVN;  and  derived 
angular  accelerations  RAN  and  RAH  from  OUTVAR  and  DSCALE. 

2„  Anthropometric  data  such  as  OAFM,  OAFD,  XAFM,  ZAFM,  XASM, 
ZASM,  XAFD,  ZAFD,  XASD,  ZASD,  AHS.  (See  program  ANXRAY) 

5.  Other  data  such  as  ANGCG,  the  angle  between  the  head  principal 
coordinate  system  +X  axis  and  the  head  anatomical  coordinate  system  +X  axis; 

XCGR,  ZCGR,  which  are  the  X  and  Z  displacements  of  the  head  principal  coordinate 
system  origin  relative  to  the  head  anatomical  coordinate  system. 

The  major  calculating  algorithm  in  TROUT  is  a  subroutine  which  transforms 
accelerations  measured  at,  and  relative  to,  a  coordinate  system  at  one  point  in  a 
rigid  body  to  accelerations  which  would  be  measured  at  and  relative  to  a  different 
coordinate  system  at  a  second  point  in  the  rigid  body.  This  transformation  utilizes 
a  standardized  form  of  the  formula  for  rigid-body  kinematic  acceleration  plus  a 
rotation  of  the  new  acceleration  vector  into  the  new  coordinate  system. 

Specific  variable  deriviation  will  be  limited  to  those  variables  presented  in  the 
plots  and/or  those  variables  required  in  the  computation  of  variables  presented  in  the 
plots:  IFM,  ISM,  4FM,  4SM,  1FD,  1SD,  4FD,  4SD,  4XM,  4ZM.  See  TROUT  list 
in  Variable  Definitions  section  for  precise  definitions  and  cross-reference  to  appro¬ 
priate  figures. 

Program  TROUT  first  calculates  the  magnitude  P  and  direction  0  of  the  vector 
from  the  first  coordinate  system  origin,  where  the  acceleration  vector  is  known  to  the 
point  in  the  rigid  body  at  which  it  is  to  be  calculated. 

Then  using  the  known  acceleration  vector  along  with  the  angular  acceleration 
(RAH)  and  velocity  (RVH),  program  TROUT  calculates  the  acceleration  vector  at  the 
second  point.  This  is  the  new  acceleration  vector  but  it  is  expressed  in  terms  of  the 
old  coordinate  system.  The  final  step  is  to  perform  a  rotation  of  coordinates  by  the 
appropriate  angle  so  that  the  new  acceleration  vector  is  expressed  in  terms  of  the  new 
coordinate  system. 


1 
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Using  the  input  quantities  as  required  the  following  quantities  are  calculated: 


\  =  (XAFM  +  XASM)/2 
Zl  =  (ZAFM+ZASMJ/2 

P1  =  <  +  z“'.) 

e  =  OAFM  +  ARCTAN  (Za  /  X% ) 


(1) 

(2) 

(3) 

(4) 


where  P1  and  01  represent  the  magnitude  and  direction  of  the  vector  from  the 
center  of  the  mouth  mount  ft-  the  head  anatomical  coordinate  system. 

Then,  the  acceleration  vector  at  the  head  anatomical  coordinate  system  origin 
(1FM,  ISM)  is  calculated  from  the  mouth  acceleration  data  (AFM,  ASM)  as  follows: 

1FM  \  /COS0  -  SIN ^  I/- AFM\  /COS6l  -  SIN0a  V  /-<Rvl-^t 

J  1  l  pi  k  i  a  /"/-vr  n  |  |  RAH  J 


ISM 


SIN0  COS0I  L\  ASM  /  1  \SIN6.  COS9 


where  0  =  -OAFM. 


To  proceed  from  the  head  anatomical  coordinate  system  to  the  head  principal 
coordinate  system,  the  vector  (XCGR,  ZCGR)  is  rotated  to  the  head  principal 
coordinate  system: 


3  = 


/COS0  -  SIN0\  XCGR 
SIN0  COS#)’  ZCGR 


(6) 


where  0  -  ANGCG  and  (X2  ,  Z2)  is  the  vector  after  rotation. 

P2  =  (X‘  +  z‘  )4 

q3  =  -ANGCG  +  ARCTAN  (Z2/X2) 


(7) 

(8) 


where  P2  and  02  represent  the  magnitude  and  direction  of  the  vector  from  the  head 
anatomical  coordinate  system  origin  to  the  head  principal  coordinate  system  origin. 
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Then  the  acceleration  vector  at  the  head  center  of  gravity  is 

|4FM\  f  COS0  -  SIN0\  flFM  _/COS63-  SIN6a\  f-(RVH)5 


i4SM/  \SIN0  COS 


JlsA 


SINS  COS0 J\  RAH 


(9) 


where  0  -  ANGCG. 

The  above  procedure  is  repeated  using  the  head  acceleration  data  (AFD,  ASD)  to 
generate  (1FD,1SD)  and(4FD/  4SD).  The  new  quantities  relating  to  the  head  mount 
are: 


X3  =  (XAFD  +  XASD)/2  (10) 

Z,  =  (ZAFD  +  ZASD)2  (11) 

P3  =(x3a  +Z32)'/2  (12) 

e.  =  OAFD  +  ARCTAN  (Z./X.)  (13) 

3  3  3 

where  P3  and  Og  represent  the  magnitude  and  direction  of  the  vector  from  the  center 
of  the  head  mount  to  the  head  anatomical  coordinate  system. 


where  0  =  -OAFD 
and 


where  0  =  ANGCG 


Finally,  a  simple  rotation  from  the  head  principal  coordinate  system  to  the 
laboratory  reference  coordinate  systems  gives 

4XM\  /COS0  -  SIN0 

4Zm]  \SIN0  COS0 

where  0  =  RPH  •<  AHS  -  ANGCG  -  ir/2 


MFM' 

4SM 


(16) 
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Trout  Validation 


The  input  data  required  by  TROUT  consists  of  DSCALE,  ICONDS,  AND  OUTVAR 
data,  all  on  digital  tape. 

To  test  the  computing  software  of  TROUT  for  algorithm, programming,  and  computa¬ 
tional  errors,  the  simulated  DSCALE  and  ICONDS  data  used  to  validate  OUTVAR  and 
the  resulting  simulated  OUTVAR  output  data  (See  description  of  program  OUTVAR) 
were  used  as  input  to  TROUT*  In  addition,  the  simulated  head  principal  coordinate 
system  was  defined  to  have  its  origin  at  X  =  2cm,  Z  =  2cm  relative  to  the  head 
anatomical  coordinate  system  origin  and  +X  axis  oriented  20  degrees  clockwise  from 
the  +X  axis  of  the  head  anatomical  coordinate  system. 

Once  TROUT  had  been  run  on  the  model  data,  verification  followed  with  checkout 
of  two  specific  areas: 

1 .  The  accelerations  computed  for  the  head  anatomical  and  head  principal 
coordinate  systems  had  to  be  Hie  same  whether  computed  from  mouth  mount  or  head 
mount  data.  Both  sources  did  indeed  produce  the  same  results  for  each  of  the 
coordinate  systems  in  question. 

2.  The  fact  that  results  match  for  both  sources  does  not  mean  the  accelera¬ 
tions  are  correct.  Furthermore,  it  says  nothing  of  the  T|  anatomical  coordinate 
system.  By  checking  the  final  positions  of  each  coordinate  system  against  predicted 
values  the  entire  program  becomes  verified.  That  is,  if  all  coordinate  systems  have 
accurate  positions  then  the  accelerations  and  velocities  leading  to  these  positions 
are  assumed  accurate. 

Table  5  gives  the  x  and  z  coordinates  of  position  of  each  coordinate  system  origin 
relative  to  the  laboratory  reference  coordinate  system. 
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TABLE  5 


Name 

TROUT  Positions  at  1/2  sec  (centimeters) 

Variable  Predicted  Value  Computed  Value 

A 

T]  Anatomical 

3XN 

-160.93 

-160.95 

.02 

Coordinate  System 

3ZN 

4.93 

4.96 

.03 

Head  Anatomical 

3XM 

-160.93 

-160.91 

.02 

Coordinate  System 

3ZM 

4.93 

4.97 

.04 

Head  Principal 

6XM 

-158.11 

-158.08 

.03 

Coordinate  System 

6ZM 

4.69 

4.71 

.02 

Any  significant  errors  in  4XM  and  4ZM  which  are  acceleration  variables  would 
manifest  themselves  as  part  of  the  errors  in  6XM  and  6ZM  which  are  final  position 
variables.  It  can  be  seen  by  examination  of  the  final  column  that  there  is  essentially 
no  error  introduced  into  the  data  by  the  TROUT  program. 

TROUT  is  a  multi-run  program  processing  up  to  30  runs  at  a  time.  The  output 
from  TROUT  is  a  digital  tape  containing  up  to  36  variables  for  each  of  the  runs 
processed.  These  tapes  were  and  are  preserved  as  a  primary  reference. 

Figure  DP-  18  gives  a  summary  illustration  of  the  data  flow  from  DSCALE  to 
TROUT. 


118 


•  -  -v. ' v ;.gr?  yj ' y^v- sv~ «>wirciwa i  ■■. < ’yyH .  t\ •  |j  -.ii)fiuw,i4u-i 


PHOTOGRAPHIC  OUTPUT  VARIABLES  AND  TRANSFORMATIONS 

TROUTP  -  The  Photographic  Data  Transformation  Program 

TROUTP  is  a  FORTRAN  program  which  computes  the  motion  of  the  head 
anatomical  coordinate  system,  the  Tj  anatomical  coordinate  system  and  the  head 
principal  coordinate  system  from  photographic  data.  Each  system's  motion  is 
described  in  terms  of  the  X  and  Z  position  of  its  origin  relative  to  the  sled  reference 
coordinate  system. 

Furthermore,  the  position  vectors  to  the  head  anatomical  coordinate  system  and 
the  head  principal  coordinate  system  from  the  origin  of  the  Tj  anatomical  coordinate 
system  are  computed.  These  position  vectors  are  in  terms  of  magnitude  and  orientation 
relative  to  the  +X  axis  of  the  laboratory  reference  coordinate  system. 

Input  to  TROUTP  consists  of  the  following: 

1 .  The  scaled  X  and  Z  coordinates  of  photographic  targets  1  through  4  and 
their  corresponding  time  in  msecs  (TMS)  from  approximately  250  frames  of  film. 

These  coordinates  are  relative  to  the  Telereadex  reference  coordinate  system,  (See 
Photographic  Data  Processing) 

2.  The  output  from  the  ICONDS  program  which  includes  the  ANXRAY 
anthropometric  data. 

3.  The  components  of  displacement  of  the  origin  of  the  head  principal 
coordinate  system  relative  to  the  head  anatomical  coordinate  system. 


Specific  discussion  of  the  variables  computed  in  TROUTP  will  be  limited  to  those 
variables  presented  in  the  plots  and/or  those  variables  required  in  the  computation 
of  variables  presented.  See  the  TROUTP  list  of  the  Variable.  Definitions  Section 
for  precise  definitions  and  cross-reference  to  appropriate  figures. 

The  first  function  of  TROUTP  is  to  express  the  input  photographic  data  relative  to 
the  sled  reference  coordinate  system.  TROUTP  rotates  the  input  coordinates  through 
the  angle  ALPHA  (ICONDS  input)  which  is  defined  to  be  the  magnitude  of  the 
angular  displacement  of  the  Telereads  *•  reference  coordinate  sy'em  relative  to  the 
laboratory  reference  coordinate  system.  The  rotated  values  ar  ius  relative  to  the 
sled  reference  coordinate  system  since  their  origin  is  the  sled  pU  graphic  target 
and  their  orientation  is  that  of  the  laboratory  reference  coordinate  system. 

Once  the  photographic  data  has  been  expressed  in  the  sled  reference  coordinate 
system,  the  displacement  of  the  head  target  center  relative  to  the  mouth  target 
center  is  computed: 

HSP  =  [  (X3  -  X4)S  +  (Z3-Z4)3]*  (1) 

Angles  RN  and  ADH  are  then  computed: 

RN  =  ARCTAN  [(Z2  -  Z1)/(X2  -  XI)]  (2) 

ADH  =  ARCTAN  [(Z3  -  Z4)/(X3  -  X4)]  (3) 
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RN  is  the  angle  formed  by  the  line  connecting  the  two  T j  targets  with  the  X  axis 
of  the  laboratory  reference  coordinate  system.  (Figure  DP-  11),  ADH  is  the  angle 
formed  by  the  line  connecting  the  head  target  and  the  mouth  target  with  the  X  axis 
of  the  iaooratory  reference  coordinate  system  .  (Figure  DP-11) 

RN  is  used  to  develop  the  angle  NRN  which  is  the  angular  displacement  of  the 
+X  axis  of  the  T]  anatomical  coordinate  system  relative  to  the  +X  axis  of  the 
laboratory  reference  coordinate  system. 

NRN  =  RN  -(-0.2  -  OAFN)  (4) 

where  OAFN  (ANXRAY  input)  is  the  angular  displacement  of  the  +X  axis  of  the 
transducer  T]  instrumentation  coordinate  system  relative  to  the  +X  axis  of  the  T ] 
anatomical  coordinate  system;  and,  from  the  geometry  of  the  spine  mount,  the  line 
connecting  the  two  T]  targets  has  an  angular  displacement  of  0.2  radians  relative 
to  the  +X  axis  of  the  transducer  Tj  instrumentation  coordinate  system. 

ADH  is  used  to  develop  the  angle  PDH  which  is  the  angular  displacement  of  the 
+X  axis  of  the  head  anatomical  coordinate  system  relative  to  the  +X  axis  of  the 
laboratory  reference  coordinate  system. 

PDH  =  ADH  -  (  tt/2  -  AHS)  (5) 

where  AHS  (ANXRAY  input)  is  the  angular  displacement  of  the  -X  axis  of  the 
photographic  head  instrumentation  coordinate  system  relative  to  the  +Z  axis  of 
the  head  anatomical  coordinate  system. 

Having  computed  the  orientations  of  the  anatomical  coordinate  systems  relative 
to  the  laboratory  reference  coordinate  system,  rROUTP  computes  the  components  of 
displacement  (XNA,  ZNA)  of  the  origin  of  the  Tj  anatomical  coordinate  system 
relative  to  the  origin  of  the  sled  reference  coordinate  system.  From  the  rotated  X 
and  Z  coordinates  of  the  two  Tj  targets,  the  coordinates  of  the  center  of  the  spine 
mount  relative  to  the  sled  reference  coordinate  system  are  computed: 

XN  =  (XI  +  X2)/2  (6) 

ZN  =  (Z1  +  Z2)/2  (7) 

Similarly,  from  XAFN,  ZAFN,  XASN,  ZASN  (ANXRAY  input),  the  coordinates 
of  the  center  of  the  spine  mount  (XNC,  ZNC)  relative  to  the  T ]  anatomical 
coordinate  system  are  computed: 

XNC  =  (XAFN +XASN)/2 
ZNC  =  (ZAFN +ZASN)/2 


(8) 

(9) 


Therefore,  the  coordinates  XNA,  ZNA  are 

|XNA\  Jx N 
iZNA/  ~(ZN 

where  0  -  NRN. 


COS 0  -  SI N/3f  (XNCl 
SIN0  COS0)m  \ZNC/ 


00) 


The  next  function  of  TROUTP  is  to  compute  the  components  of  displacement  of  the 
origin  of  the  head  principal  coordinate  system  (XCL,  ZCL)  relative  to  the  origin  of 
the  sled  reference  coordinate  system.  From  the  rotated  X  and  Z  coordinates  of  the 
head  target  and  of  the  mouth  target,  the  coordinates  of  the  mid-point  of  the  line 
connecting  the  head  and  mouth  targets  relative  to  the  sled  reference  coordinate 
system  are  computed: 


X34  =  (X3+X4)/2 
Z34  =  (Z3+Z4)/2 


(11) 

(12) 


(See  Figure  DP-13)  Similarly,  from  XAFM,  XASM,  XAFD,  XASD,  ZAFM,  ZASM, 
ZAFD  and  ZASD  (ANXRAY  input),  the  coordinates  of  the  midpoint  of  the  head- 
mouth  line  (XC,  ZC)  relative  to  the  head  anatomical  coordinate  system  are  computed: 


XC  =  (XAFM  +  XASM  +  XAFD  +  XASD)/4 
ZC  =  (ZAFM  +  ZASM  +  ZAFD  +  ZASD)/4 


(13) 

(14) 


The  coordinates  of  the  origin  of  the  head  anatomical  coordinate  system  (XHA,  ZHA) 
relative  to  the  sled  reference  coordinate  system  may  be  derived  as  follows: 


|XHAL/X34)  !COS0  -  SIN0]  | 
[ZHAp[Z34 J"  \SIN0  COS  01 

where  0  -  PDH. 


XC 

ZC 


(15) 


Finally,  if  XCG,  ZCG  represent  the  coordinates  of  the  origin  of  the  head 
principal  coordinate  system  relative  to  the  head  anatomical  coordinate  system  (Input 
to  TROUTP), 


/XCL  __ 

/XHAi  +  1 

(ZCL  ' 

IZHA)  +  l 

(16) 


where  0  =  PDH 


From  XNA,  ZNA  and  XCL,  ZCL,  TROUTP  computes  the  magnitude  of  displacement 
of  the  head  principal  coordinate  system  origin  relative  to  the  T]  anatomical  coordinate 
system  origin. 


MCT  =  [(XCL  -  XNA)2  +  (ZCL  -  ZNA)2)  % 


(17) 
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The  polar  angle  (ACT)  of  MCT  is  also  computed: 

ACT  =  ARCTAN  [(ZCL  -  ZNA)/(XCL  -  XNA)] 

It  should  be  noted  that  a  special  version  of  TROUTP,  TROUTH,  calculates  the 
same  variables  using  the  same  algorithms,  but  is  designed  to  accept  as  input, 
scaled  photographic  data  obtained  from  the  filming  of  the  subject's  normal  head  nod 
rather  than  the  usual  photographic  data  from  the  acceleration  experiments. 

TROUTP  Validation 

Validation  of  TROUTP  was  based  on  simulated  data  generated  by  a  model  essentially 
similar  to  the  model  used  to  generate  data  for  OUTVAR  validation.  (See  program 
OUTVAR  for  detailed  description  of  the  model)  There  were,  however,  certain 
differences  necessary  to  accommodate  photographic  input  data  and  the  structure  of 
the  program  TROUTP, 

The  first  difference  was  that  no  linear  acceleration  was  used.  The  model  testing 
OUTVAR  and  TROUT  used  a  linear  acceleration  a(t)  as  well  as  a  rotation  a>(t)  as 
driving  functions.  a(t)  was  necessary  since  sensor  data  consisted  of  acceleration 
signals,  one  component  of  which  was  a  linear  acceleration  which  ultimately  would 
contribute  to  position  relative  to  a  fixed  lab  frame.  Photographic  data  is  referenced 
to  a  moving  frame  (the  origin  of  the  sled  reference  coordinate  system)  and  hence 
does  not  involve  linear  translation.  Therefore,  only  co (t)  was  chosen  as  a  driving 
function,  the  same  a)(t)  that  was  used  before. 

6o(t)  =  (-800  degrees/sec3)^  for  0  ss  t  250  msec  (19) 

co(t)  =  (800  degrees/sec  )t  -  400  degrees/sec 

for  250  msec  s  500  msec  (20) 

The  second  difference  involved  the  internal  geometry  of  the  model.  The  T| 
anatomical  system  was  chosen  to  be  at  coordinates  (2cm,  2cm)  relative  to  a  coordinate 
system  formed  by  the  Juction  of  L3  and  Ls  and  having  axes  along  these  lines.  The 
orientation  of  the  Tj  anatomical  coordinate  system  would  coincide  with  the  laboratory 
reference  system  as  an  initial  condition.  The  head  anatomical  coordinate  system  as 
used  in  testing  TROUT  stayed  the  same  (  i.e.  2cm  to  the  right  of  the  junction  of  L 
and  L  on  L  )  with  the  orientation  of  the  laboratory  reference  coordinate  system  as 
an  initial  condition.  The  head  principal  coordinate  system  was  chosen  to  be  at 
(2cm,  2cm)  relative  to  the  head  anatomical  coordinate  system  .  Figure  DP-19 
represents  all  of  these  positions  at  t  =  0. 
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Origin  of  Head 
Principal  Coordinate  System 
/ 


The  positions  of  the  various  coordinate  systems  at  t  =  500  msec  were  to  be 
compared  with  predicted  values  from  the  model.  Due  to  a  fault  in  a  printing  routine 
it  was  necessary  to  perform  the  check  at  the  250th  data  point  (correspond  to  498 
milliseconds  into  the  run)  instead  of  the  251st  point.  The  following  is  a  table  of 
computed  values  against  predicted  values,  for  TROUTP  variables  referenced  in  this 
presentation.  It  can  easily  be  seen  that  all  errors  are  on  the  order  of  10” 6  centi¬ 
meters  and  probably  represent  errors  due  to  the  truncation  in  the  predicted  values. 
These  are  considered  insignificant. 
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TABLE  5 


COMPARISON  OF  COMPUTED  vs  PREDICTED  VALUES 
(at  t  =  498  msec) 


Variable 

Name 

Computed 

Value 

Predicted 

Value 

|A| 

PDH 

-*87263  rad 

-.87263  rad 

• 

XHA 

IV.  47774  cm  * 

19.47775  cm 

lxlO"13 

ZHA 

4.96160  cm 

4.96156 

4x1 0“6 

NRN 

-.87264  rad 

-.87264  rad 

- 

XNA 

2.56346  cm 

2.56346  cm 

m 

e 

ZNA 

-1.19527  cm 

-1.19526  cm 

1x10 

XCL 

22.29541  cm 

22.29540  cm 

lxl0"B 

ZCL 

4.71518  cm 

4.71518  cm 

■a 

MCI 

20.59814  cm 

20.59813  cm 

lxlQ“6 

ACT 

.29103  rad 

.29103  rad 

•m 
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FIGURE  D 


POST-EXPERIMENT  CONDITIONING  OF  CERTAIN  DATA 
PRESENTED  IN  THIS  REPORT 


Gilbert  Willems  M.S.E.E.  and  William  Anderson  M.A.PHys. 
INTRODUCTION 

'  '  The  majority  of  the  data  collected  during  the  human  impact  experiments 
conducted  at  Wayne  State  University  was  of  sufficient  quality  to  be  useable  in 
post-experiment  analysis  without  further  conditioning.  However,  the  quality  of 
a  small  subset  of  this  data  was  found  to  be  enhanced  considerably  by  judicious 
application  of  certain  conditioning  techniques,  such  as  filtering  and  time  correlation; 
these  techniques  and  the  justification  for  their  use  are  detailed  in  the  remainder  of 
this  section. 

APPLICATION  OF  FILTERING  TECHNIQUES 

Filtering  was  introduced  into  some  of  the  data  to  perform  one  of  two  functions: 
reduction  of  excessive  noise  or  as  a  pre-smoother  prior  to  performing  differentiation. 
Of  the  recorded  sensor  data,  only  the  sled  acceleration  traces  were  sufficiently 
noisy  to  require  filtering.  This  excessive  noise  has  been  attributed  to  the  fact  that 
an  undamped,  high  bandwidth  accelerometer  was  hard  mounted  to  the  sled;  conse¬ 
quently,  not  only  was  its  own  natural  frequency  excited,  but  sled  resonances  were 
also  sensed,  yielding  traces  sufficiently  noisy  to  hamper  data  analysis. 

Several  variables,  particularly  photometrically  derived  displacements  were 
differentiated.  The  original  variables  while  not  sufficiently  noisy  as  ro  require 
filtering,  did  contain  some  high  frequency  jitter,  which  would  tend  to  be  accentuated 
by  differentiation.  This  latter  phenomenon  can  readily  be  understood  heuristically 
if  note  is  made  of  the  fact  that  the  derivative  of  a  function  at  any  given  point  is  a 
measure  of  the  slope  of  the  function  at  that  point.  Therefore,  since  high  frequencies 
are  associated  with  steep  slopes,  the  noise  exaggeration  effect  of  differentiation 
occurs.  To  minimize  this  effect,  ail  variables  subjected  to  differentiation  were  first 
appropriately  filtered. 

Selection  of  a  filtering  criterion  is  dependent  on  two  important  requirements: 

1 .  The  filter  bandwidth  must  be  sufficient  to  "pass"  all  significant  signal 
components. 

2.  Distortion  in  time  due  to  phase  shift  must  be  avoided. 

Appropriate  filter  cutoff  frequencies  were  selected  by  first  performing  a  spectral 
analysis  via  standard  Fast  Fourier  Transform  (FFT)  techniques,  in  order  to  establish 
the  spectral  location  of  the  significant  signal  components.  Time  axis  distortion  was 
avoided  by  properly  accounting  for  time  delays  due  to  filter  phase  shift,  and 
correcting  the  time  axis. 
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The  filters  chosen  for  this  effort  are  representable  by  one  of  the  following  two 
second  order  transfer  functions: 


output 

input 


52  +  ki 


Wi  toi 


s  +  1 


S3  +2£a 


(l)  or  -SH*Ei!L  = 
input 


1 


(2) 


s  +  1 


Ub  Ui 


$3 


2C. 

ox, 


s  +  1 


In  these  transfer  functions  Wl  and  w  g  are  the  numerator  and  denominator 

break  frequencies  respectively  (In  radians/sec),  and  £  and  £a  are  the  damping 

ratios;  s  is  the  Laplace  operator.  The  second  characterization  allows  a  40  db/decade 
rolloff  beyond  the  passband.  The  Inclusion  of  the  numerator  term  on  the  first  transfer 
function  allows  the  insertion  of  an  additional  "notch"  anywhere  along  the  frequency 
axis,  ParticaSIy,  however,  the  insertion  of  the  notch  is  feasible  only  if  inserted 
significantly  beyond  the  denominator  breakpoint,  since  the  rolloff  is  changed  to  a 
constant  attenuation  beyond  the  notch. 

The  derivation  of  a  digital  algorithm  to  represent  the  indicated  transfer  functions 
is  readily  obtained  from  the  state  space  representation  of  the  cited  filters,  both  of 
which  are  represented  in  block  diagram  form  in  Figure  1 .  In  this  figure  the 
integrator  outputs  are  defined  as  state  variables  x1  and  xa  .  If  a  low  pass  mechani¬ 
zation  is  desired,  is  used  as  the  output;  if  in  addition,  a  notch  is  desired,  E0ut  is 
the  output.  From  the  figure,  the  vector-matrix  filter  equations  are  readily  obtained: 


L  X8  j 


1 


*1 

X,* 

'  0  ' 

1 

+ 

J 

~*s. 

-Was. 

LcOs3  “2C 3  ^8 * 

E0ut  is  obtainable  as  a  linear  combination  of  xx  and  xg  : 

Eout  =  ^  Ei«  +^2CiC01-2£8ie8) 


(3) 


(4) 


<0- 


The  above  vector-matirx  differential  equation  is  of  the  general  form 

><  =  Ax  +  Bu 

whose  solution  is  given  by: 

x(t)  =  e  x(o)  +  J  Q  e  "  ^  Bu  C7)  d^ 


(5) 

(6) 
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Since  the  data  of  interest  is  discrete  (due  to  the  digitization  process)  the  solution 
for  the  discrete  form  of  subject  set  of  differential  equation  is  needed.  A  recursive 
relation  exists  that  allows  me  computation  of  the  system  state  at  time  (k  +  1)T  based 
on  knowledge  of  the  state  at  time  icT;  k  is  a  running  index  and  T  is  the  sampling 
interval.  The  recursive  relation  is  given  by  a  formula.  (Reference  1) 

x(k  +  l)T  =  G(T)  x  (kT)  =  H(T)  u  (kT)  (7) 

where  G(T)  =  eAt;  H(T)  =  J  T  e  ^  dr  B  (8) 

The  solution  of  the  subject  difference  equation  reduces  essentially  to  computing 

e  and  its  integral  over  T ,  a  tedious,  but  straightforward  task  which  is  rendered 
simpler  by  the  following  substitutions: 


Ca  40a  =  a  >  60  =  cu3 
The  resulting  solution  then  becomes: 

x  i  (k+l)  'Au  Aja*  xx(k) 

_x  (k+1)  A^  A^  _x5(k)_ 


where: 


rI  - r? 


<4?  F» 
Wa2  Faa 


An  *  e  -QfT (cos  WT  +-J-sin  wr  ) 
A,  =  —  (  e  sin  coT  ) 


A  =  -  W  +  Wa> 

. . .  u 

C 0 


“  sin  cut] 


= 


e  “  [  cos  ^T  -  a/ w  sin  cuT  ] 


F*  = 


CO  (a?  +  co2) 


(asin  co T  +  co  cos  coT)  +  1 

<?  +  CO3 


F*  =  Al. 
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(17) 


From  the  original  equation,  E  .(ic)  is  inferred  to  be 


Eout<k> 


E,.<k>  +f' 


x  (k)  +  - 


x3(k) 


Therefore,  the  digital  algorithm  consists  of  mechanizing  the  three  equations 
for  x^k  +  1),  x3(k  +1)  and  £0(k).  Required  inputs  are  the  sequence  to  be  filtered 

EJn(k)#  and  the  initial  states  x  ^0)  and  x  a (0) .  x^O)  represents  the  value  of 

E 4  (k)  at  t  =  0,  whereas,  x_(0)  is  the  initial  slope  of  E,  (k).  The  latter  is  often 

1  S  2  1  ltt 

not  known,  in  which  case  it  can  be  estimated  or  set  to  zero;  this  may  sometimes 
cause  a  transient  for  a  few  sampling  intervals,  a  common  filter  initialization  problem. 
In  actual  filtered  data  x8(0)  was  zero  and  transients  have  not  presented  a  problem. 

As  previously  noted,  amplitude  distortion  is  avoided  by  selecting  a  filter 
breakpoint  that  lies  beyond  the  significant  signal  frequencies  as  determined  via 
spectral  analysis.  Phase  shift,  however,  poses  a  different  problem  since  there  is 
always  some  phase  shift  for  all  non-zero  frequencies.  Filter  phase  shift  causes  a 
time  delay  in  the  data,  defined  by: 


Time  delay 


T4  (secs) 

F 

given  frequency 


08) 


where  /  =  phase  shift  in  degrees 

F  =  frequency  in  degrees/sec 

It  is  obvious  from  the  above  equations  that  if  we  can  somehow  obtain  a  linear  phase- 
frequency  relationship,  then  T  is  a  constant  since  it  would  be  defined  by  the  slope 
of  the  straight  line  represented  by  the  Tfl  equation  given  above.  Filters  with  such 
a  linear  phase  characteristic  in  the  signal  region  are  called  "constant  group  delay 
filters"  since  each  frequency  component  is  delayed  the  same  amount  T4  .  The  only 
distortion  caused  by  such  filters  is  a  pure  time  shift  (delay)  by  the  amount  Td;  there¬ 
fore,  the  data  are  easily  corrected  by  applying  an  equivalent  leading  shift  to  the 
time  axis. 

The  phase  shift  for  any  second  order  transfer  characteristic  is  given  by: 

/ (degrees)  =  tan"1/ 2^^  \ 

\  CO^-OJ8/ 

Where  g  and  ^  are  the  damping  ratio  and  natural  frequency  of  the  ^th  second 
order  transfer  characteristic  and  w  is  the  particular  frequency  for  which  we  apply 
fi  .  For  a  numerator  term,  /  is  positive,  and  negative  for  a  denominator  term. 
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Phase  linearity  is  a  function  of  the  specific  filter  parameters  selected. 

Appropriate  parameters  can  be  selected  heuristically  without  undue  difficulty  with 
the  aid  of  normalized  frequency  response  plots  for  the  general  filter  forms.  (2)  All 
filters  used  for  data  presented  in  this  report  were  selected  with  above  considerations 
in  mind.  For  example,  the  sled  accelerometer  data  were  filtered  twice;  first  by  a 
iOO  Hz  filter  with  a  notch  at  500  Hz,  then  by  a  40  Hz  low  pass  section  notched  at 
100  Hz.  The  composite  frequency  response  is  shown  in  Figure  2,  where  it  is  seen 
that  flat  response  nas  been  preserved  well  beyond  20  Hz,  the  established  upper  limit 
of  signal  frequencies.  The  phase  shift  in  this  range  is  shown  in  Figure  3,  which 
establishes  that  good  phase  linearity  has  been  obtained. 

The  actual  filter  parameters  for  each  filtered  variable  included  in  this  report 
are  shown  in  Table  1.  In  each  case  the  time  delay  was  automatically  computed 
during  the  filtering  process  and  the  appropriate  corrective  time  shift  introduced 
in  the  respective  time  axis.  Sample  "before  and  after"  traces  for  sled  acceleration 
are  shown  in  Figures  4  and  5  . 

RUN  CLASSIFICATION  CRITERIA 

Due  to  the  variability  of  accelerator  performance,  a  continuum  of  sled  accelera¬ 
tion  and  rates  of  onset  were  obtained  at  Wayne  State  University,  rather  than  the  few 
preset  discrete  levels.  For  this  reason,  some  grouping  logic  was  required  so  that 
manageable  subsets  of  data  could  be  obtained.  The  approach  used  was  to  establish 
several  "nominal"  levels  of  each  variable  as  follows: 

"a"  level  -  2,  2.5,  3,  4,  5,  6,  7,  8,  9,  10 

rate  of  onset  (g's/sec)  -  100,  150,  200,  250,  300,  350,  400,  500, 

600,  700,  800,  900 

This  procedure  required  assignment  of  a  specific  sled  acceleration  and  rate  of 
onset  level  to  each  run,  which  presented  an  additional  problem  due  to  the  nature 
of  the  sled  acceleration  data.  The  filtered  accelerometer  traces  such  as  shown  in 
Figure  5  are  a  considerable  improvement  over  the  original  data.  With  few  exceptions, 
the  onset  side  is  of  sufficient  quality  for  accurate  determination  of  rate  of  onset. 
However,  there  are  enough  variations  (low  frequency  oscillations)  around  the  peak 
to  make  an  accurate  assessment  of  the  time  of  occurrence  of  the  peak  acceleration 
difficult.  Therefore,  a  new  algorithm  was  created  to  yield  the  "best"  (in  a  least 
squares  sense)  ideal  triangular  waveform  approximation  to  the  filtered  sled  accelera¬ 
tion  traces.  This  algorithm  consisted  of  the  following: 

1.  Determine  the  peak  acceleration  based  on  the  filtered  data. 

2.  Compute  the  20%  and  80%  points. 

3.  Using  all  data  between  these  points,  compute  the  best  straight  line  fit 
in  a  least  squares  sense.  This  was  done  independently  for  the  positive  and  negative 
slopes  of  the  traces. 

4.  The  intersection  of  these  lines  was  used  to  define  the  peck  acceleration 
and  the  time  of  occurrence  of  the  peak.  The  positive  slope  of  the  fitted  line  defined 
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the  rate  of  onset.  This  technique  was  evaluated  by  comparing  it  with  previously 
used  ones  and  found  to  be  significantly  superior  to  previous  data  in  that  excellent 
time  correlation  among  similar  runs  was  obtained. 

The  data  obtained  were  fitted  into  one  of  the  previously  mentioned  categories 
if  it  was  within  +  10%  of  same.  If  not,  it  was  left  unclassified.  in  those  instances 
where  the  +  10%"criteria  overlapped  into  two  possible  nominal  categories,  the  data 
were  fittecTto  the  group  to  which  it  was  nearest  auantitatively.  For  purposes  of  this 
report,  the  above  criteria  were  broadened  somewtiat,  in  order  to  include  a  repre¬ 
sentative  number  of  experiments  in  each  desired  category.  For  example,  runs 
nominally  assigned  to  the  250  g  and  350  g  levels  were  grouped  with  the  300  g  runs, 
and  assigned  to  the  latter  level. 

Selection  of  sled  acceleration  parameters  for  the  data  presented  in  this  report 
was  based  on  the  above  criteria,  since  these  have  yielded  groupings  with  very  small 
dispersions  in  terms  of  peak  acceleration,  rate  of  onset,  arid  time  to  peak  for  similar 
run  groupings. 

TIME  CORRELATION  BETWEEN  PHOTO  AND  SENSOR  DATA 


Two  different  systems  were  used  to  collect  dr,ra  during  the  WSU  impact  accelera¬ 
tion  experiments:  (1)  the  photo  system;  and  (2)  the  sensor  system.  Within  either 
system  all  the  recorded  data  were  found  to  be  well  aligned  with  respect  to  time. 
Between  systems,  however,  many  runs  have  exhibited  a  considerable  misalignment, 
a  phenomenon  uttributable  to  the  fact  that  a  central,  synchronized  timing  system 
such  as  IRIG  was  not  employed. 


It  was  necessary,  therefore,  to  find  an  accurate  estimate  for  the  timing  misalign¬ 
ment  between  systems,  which  was  performed  as  follows:  The  values  for  angular  dis¬ 
placement  of  the  head  derived  from  the  photo  data  (ADH)  are  generally  in  good 
agreement  with  the  values  derived  from  the  rafro  gyro  (RPH),  except  for  a  time  shift. 
To  determine  the  shift,  the  photo  values  were  compared  with  the  rate  gyro  values 
for  all  time  shifts  +  50  msec  from  nominal  in  0.5  msec  increments.  The  criterion 
used  for  evaluating  the  comparison  was  the  cross  correlation  coefficient,  calculated 
as 


P  xy  (t)  = 


Where 


n 


2x. 

Ex  s-  1 


n 


Eyt  s 
n 


(20) 


P  xy  (t)  =  cross  correlation  coefficient  as  a  function  of  time  shift, 
xl  =  ith  photo  data  point. 
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yi  =  sensor  data  point  correspondent  to  the  tth  photo  data  point  after 
photo  data  has  been  shifted. 

n  =  total  number  of  photo  data  points. 

f  =  time  shift  of  photo  data  ranging  from  -  50  msec  to  +  50  msec,  In 
0.5  msec  increments. 

As  P  xy  (t)  increases  to  1  the  relationship  between  x  and  y  values  becomes  more 
closely  approximated  by  an  equation  of  the  form 

xi  -  °y t  +  b  (21) 

As  P  xy  (t)  decreases  from  1  the  linear  relationship  becomes  less  able  to  account 
for  differences  between  the  x  and  y  values. 

Thus,  the  best  agreement  between  the  two  data  series  occurs  at  the  time  shift  r 
at  which  P  xy  (r)  is  maximum. 

A  computer  program  was  written  to  calculate  P  xy  (t)  for  the  different  T  values 
for  each  run,  and  to  list  the  run  number,  the  best  estimate  of  the  time  shift  for  the 
photo  data,  and  the  corresponding  cross  correlation  value. 

Results  for  the  data  presented  in  this  report  are  summarized  in  Table  2. 
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TABLE  2 


RESULTS  OF  THE  TIME  CORRELATION  ANALYSIS  BETWEEN  PHOTO 
AND  SENSOR  DATA  FOR  HEAD  ANGULAR  DISPLACEMENT 


SUBJECT  NO. 

PEAK  G 

ROO 

PHOTO  TAU 
(Correction  to 
Time  Axis ) 

CROSS  CORRELATION 

007 

6.8 

360 

6.000 

1.000 

004 

6.9 

390 

5.500 

1.000 

008 

5.5 

300 

1.500 

.999 

011 

5.5 

280 

6.500 

1.000 

016 

5.5 

250 

4.000 

1.000 

010 

7.2 

320 

4.000 

1.000 

012 

5.7 

340 

5.500 

1.000 

016 

7.4 

450 

5.000 

1.000 

004 

8.9 

630 

4.000 

.999 

007 

8.5 

350 

3.500 

1.000 

004 

5.9 

330 

8.000 

.999 

017 

5.6 

340 

2.500 

1.000 

011 

8.3 

300 

-3.500 

1.000 

010 

8.9 

330 

-3.500 

.999 

012 

8.1 

370 

3.000 

.999 

003 

5.9 

560 

11.000 

1.000 

017 

8.0 

410 

1.500 

.999 

011 

5.8 

510 

9.000 

1.000 

015 

7.0 

380 

2.500 

1.000 

013 

7.2 

320 

3.500 

1.000 

017 

8.9 

560 

2,000 

1.000 

016 

5.8 

630 

9.000 

1.000 

011 

8.9 

380 

-2.500 

1.000 

003 

7.2 

400 

1.000 

1.000 

016 

8.4 

440 

.000 

.999 

010 

6.0 

560 

8.500 

1.000 

004 

6.1 

620 

12.000 

1.000 

009 

8.4 

560 

3.000 

1.000 
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TABLE  2  (Continued) 


SUBJECT  NO.  PEAK  G  ROQ  PHOTO  TAU  CROSS  CORRELATION 


RESULTS  CHAPTER 


RUN  GROUPINGS 


All  runs  In  this  report  were  conducted  at  Wayne  State  University  between 
30  January  and  23  April!  969. 

For  this  report,  forty-one  runs  were  selected  on  the  basis  that:  (a)  no  transducer 
or  camera  failures  had  occurred  during  the  run;  (b)  all  data  were  available  from  both 
photographic  and  transducer  instrumentation  systems;  (c)  the  runs  were  fully  instru¬ 
mented  with  a  sled  peak  acceleration  category  of  6  G  or  higher;  (d)  all  data 
processing  had  been  completed  in  time  to  include  the  runs  in  this  report.  There  are 
numerous  additional  runs  which  were  excluded  by  not  appearing  in  all  4  categories. 

The  runs  were  grouped  by  volunteer  sitting  height  percentile  and  by  subject  and 
sled  peak  acceleration  category  as  noted  in  Table  1,  and  by  volunteer  as  shown  in 
Table  2. 

To  show  repeatability,  two  nearly  identical  runs  each  on  three  subjects  are 
presented.  One  subject  had  two  pairs  of  nearly  identical  runs.  The  subjects  and 
their  runs  are  presented  in  Table  3. 

The  overall  results  are  presented  in  755  black  and  white  computer-drawn  plots, 
numbered  consecutively  from  1  through  755  and  individually  labeled.  They  are 
listed  in  the  plot  table  of  contents  and  presented  in  a  separate  plot  section.  The 
variable  symbols  are  defined  In  detail  in  the  Theoretical  Kinematics  section  of  the 
report.  Careful  attention  was  given  to  clarity  of  presentation  to  facilitate  com¬ 
parison  of  the  curves  on  each  plot.  In  several  cases  this  required  plotting  fewer 
data  words  than  available  on  the  data  tapes  to  prevent  overlap  of  the  curve  symbols. 
The  appearance  of  the  data  Is  unchanged  as  seen  by  comparing  identical  variable 
curves  in  which  both  every  word  was  plotted;  and  every  sixth  word  was  plotted,  The 
plotting  density  is  as  follows: 

Plot  Plot  Density 


1  -  36 

Every  sixth  data  word  (N=6) 

78-86 

128  -  136 

367  -  429 

178  -  366 
430  w  516 
521,  568,  592 
612  -  632 


1  Plot  Density 

Every  other  data  word  (N5^) 


3 7  -77  .  IV; 
87  -  127  :r 
137-177 
493  -  520 
522  -  567 
569  -  591 
593  -  611 
633-  755 


Every  data  word  (N=l) 


N  x  0.5  msec,  where  N  is  indicated  above.  For  variables  derived  from  photographic 
data  the  relation  is  N  x  2.0  msec.  ...* 

Variable  Discussion 

Plots  1  through  27  present  the  performance  of  the  accelerator,  grouped  by  sled 
peak  acceleration  and  rate  of  onset  of  acceleration.  The  sled  acceleration  as 
presented  has  been  filtered  as  explained  in  the  section  oh  signal  conditioning. 
Because  of  the  demonstrated  repeatability,  the  grouping  by  sled  peak  acceleration 
and  rate  of  onset  of  acceleration  is  an  excellent  way  of  establishing  categories  for 
presentation  of  the  variables  measured  from  the  experimental  subjects. 

Plots  28  through  177  present  the  Cartesian  components  and  magnitudes  of  the 
linear  acceleration  computed  at  the  head  C.  G.  Each  variable  is  presented  grouped 
first  by  sled  peak  acceleration  category  and  then  with  the  same  curves  as  individual 
plots.  The  Cartesian  components  4XM,  4ZM  are  quite  detailed  and  complex.  The 
peaks  and  time  of  peak  repeatability  is  illustrated  on  the  plots  grouped  by  sled 
acceleration  category.  The  resultant  magnitude,  4XZ,  of  the  component  4XM  and 
4ZM  has  two  characteristic  peaks.  The  first  peak  is  sharp  occurring  at  about  75 
milliseconds  after  OATZ  for  all  runs.  The  repeatability  within  sled  acceleration 
grouping  is  illustrated.  The  maximum  resultant  of  about  27  G's  occurred  in  plot  176, 
subject  11,  8.9G,  380  g/sec. 


The  linear  velocity  data  measured  on  the  subject  is  presented  on  plots  178 
through  303.  Each  variable  is  plotted  grouped  by  subject  and  then  grouped  by  sled 
peak  acceleration.  The  Cartesian  components  and  mugnitude  of  the  head  C.  G. 
linear  velocity  relative  to  the  sled  are  represented  by  XCV,  ZCV,  and  XZV.  The 
magnitude  XZV  demonstrates  two  successive  peaks  which  are  repeatable  in  peak  and 
time  for  DATZ  at  a  given  run  level.  The  peaks  systematically  increase  in  value  with 
Increasing  sled  peak  acceleration  level.  The  maximum  velocity  of  about  480 
centimeters  per  second  occurred  on  plot  240  on  subjects  07,  9.6G,  630  G/sec  arid 
12,  9.8G,  530  G/sec. 

The  analogous  linear  velocities  of  the  head  relative  to  T]  are  XNV,  ZNV,  and 
NNV,  These  variables  illustrate  the  velocity  of  the  head  C.  G.  with  the  velocity 
of  the  first  thoracic  vertebra  subtracted.  The  only  change  is  d  relatively  small 
reduction  of  the  velocity  peaks. 

The  linear  displacement  of  the  head  C.  G.  rotative  to  the  anterior  superior 
corner  of  the  first  thoracic  vertebra  is  shown  in  plots  304  through  366  as  Cartesian 
components  XCN,  ZCN  and  magnitude  MCT.  The  X  component  XCN  shows 
apparent  increase  in  amplitude  with  sled  acceleration.  The  Z  component  ZCN 
decreases,  and  then  becomes  negative  or,  stated  differently,  goes  below  a  line 
perpendicular  to  gravity  and  passing  through  the  T]  origin,  with  increasing  sled 
acceleration.  This  result  was  anticipated  due  to  the  increasing  anterior  bending  of 
the  head  and  neck  with  higher  acceleration.  The  magnitude  MCT  shows  little 
change.  The  small  change  present  is  a  decrease  in  value  up  to  3  to  3.5  centimeters. 
For  these  three  variables  the  response  is  apparently  very  repeatable  and  character¬ 
istic. 

The  plots  of  MCT  as  a  function  of  time  should  be  considered  in  conjunction  with 
the  headnod  and  dynamic  event  plots  of  MCX,  MCZ  arranged  by  subject,  plots 
493  through  516.  Both  the  headnod  (HEDNOD  DATA)  and  the  dynamic  event 
(PHOTOGRAPHIC  DATA)  plots  illustrate  segments  approximating  the  arc  of  a  circle. 
There  are  no  apparent  major  differences  between  the  arcs  from  headnod  and  the 
dynamic  event.  Time  increases  with  each  data  point  from  the  upper  left  to  the  lower 
right  on  the  piofs.  The  plots  are  presented  to  illustrate  the  geometry  of  the  impact 
induced  motion  (PHOTOGRAPHIC  DATA)  as  compared  with  the  geometry  of  a 
voluntary  headnod  (HEDNOD  DATA).  Also,  these  can  be  used  with  a  protractor  to 
determine  MCT  for  any  angular  increment  desired,  with  the  zero  point  of  the  pro- 
tractor  placed  on  the  origin  of  the  plot  coordinate  system. 

The  angular  acceleration  of  the  head,  RAH  is  shown  in  plots  367  through  387. 
There  is  a  characteristic  triphasic  response  lasting  from  50  to  225  milliseconds  after 
DATZ.  The  response  is  repeatable  from  subject  to  subject  within  sled  acceleration 
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level.  The  first  peak  increases  with  increasing  sled  peak  acceleration.  The  max¬ 
imum  angular  acceleration  of  approximately  91,000  degrees  per  second  squared  was 
reached  on  subject  12,  9.8G,  530  G/sec,  plot  387. 

The  angular  velocity  of  the  head  (RVH)  is  presented  in  plots  388  through  408. 

The  response  illustrates  the  same  degree  of  repeatability  as  RAH.  Also  an  increasing 
response  with  increasing  sled  acceleration  is  evident.  The  shape  of  the  curve  and 
time, of  peaks  are  consistent  with  its  differentiated  value  RAH.  The  maximum  value 
of  about ;  1,700  degrees  per  second  is  slightly  less  than  30  radians  per  second  on  plot 
408,  subject  07,  9.6G,  630G/sec.  A  higher  peak  from  a  run  not  included  in  this 
.gteUp:(vy^;:;prqv.i6.c!S'|y  reported  (1).  ..  v  'V..-'  ..  Si 

The  angular  velocity  of  the  head  relative  to  T]  (RHN)  is  illustrated  f.% plots  409 
through  420.  RHN  is  angular  Velocity  of  head  (RVH)  minus  the  angular  velocity  of 
T|  (RVN).  The  effect  of  this  subtraction  is  minimal.  Therefore,,  RHN  has  a  shape 
and  amplitude  similar  to  RVH.  The  most  consistent  result  is  that  RHN  is  slightly 
smaller  in  peak  amplitude  than  RVH.  This  was  anticipated  since  regardless  of  the 
tightness  of  torso  restraint,  Tj  will  have  some  small  anterior  displacement  and., 
rotation  yielding  an  angular  velocity. 

Angular  displacement  of  the  head,  PDH,  is  demonstrated  in  plots  430  through  450. 
All  curves  of  PDH  show  a  characteristic  and  repeatable  pattern.  The  maximum 
change  in  PDH  from  initial  value  to  peak  value  was  91  degrees  occurring  on  plot  450, 
subject  07,  9.6G,  63QG/sec.  Also,  the  maximum  angular  displacement  is  increased 
with  increasing  peak  acceleration  of  the  sled. 

The  angular,  displacement  of  T j  (NRN)  is  shown  in  plots  451  through  471 1  These 
plots  show  a  much  smaller  angular  displacement  of  Tj  as  compared  with  the  angular 
displacement  of  the  head.  Presumably  this  is  a  function  of  the  restraint  system, 
among  other  factors. 

The  angular  displacement  of  the  head  relative  to  T]  (PHN)  is  calculated  by  sub¬ 
tracting  NRN  from  PDH,  plots  472  through  492.  PHN  has  only  minor  changes  from 
PDH.  Essentially,  the  peak  of  PHN  is  less  than  PDH  due  to  the  small  angular  dis¬ 
placement  of  T] . 

The  plots  from  1  through  516  contain  substantially  ail  the  results  of  the  report. 
Repeatability  Discussion 

Repeatability  is  demonstrated  by  plots  517  through  611.  The  runs  thus  selected 
were  paired  according  to  the  following  criteria: 


a)  identical  subject 

b)  essentially  identical  sled  peak  acceleration 

c)  similar  rate  of  onset  of  sled  acceleration 

Four  pairs  of  runs  were  thus  identified.  It  should  be  noted  that  the  initial  con¬ 
ditions  within  a  selected  pair  of  runs  were  not  precisely  controlled.  However,  they 
were  carefully  measured  and  are  represented  for  all  variables  by  the  first  value  on 
the  plots.  This  is  an  important  consideration  since  differences  in  initial  conditions 
from  run  to  run  on  the  same  subject  may  be  a  source  of  variation  in  the  ensuing 
response. 

"r- .  .  •  •  •  u 

The  repeatability  plots  include  those  of  the  acceleration  transformed  to  the  head 
C.  G.  from  the  mouth  and  head  accelerometer  pairs.  This  is  the  4FM  vs.  4FD  com¬ 
parison,  and  the  4SM  vs.  4SD  comparison.  It  is  therefore  possible,  with  the  two 
plots  from  the  paired  runs  to  compare  the  repeatability  of  several  functions. 

In  addition,  all  variables  which  are  presented  in  plots  I  through  492  except  for 
XVC  and  XEJC  are  presented  for  the  4  pairs  selected.  A  completely  repeatable 
response  would  be  illustrated  by  identical  traces  of  the  paired  observations  un  each 
plot.  Review  of  each  plot  shows  only  minor  differences  between  paired  observations. 
The  largest  of  these  differences  can  be  seen  to  be  associated  with  variations  of  the 
initial  conditions  evident  from  the  first  values  on  the  plots.  In  most  cases  each 
variable  is  repeatable  to  approximately  the  same  degree  as  the  sled  acceleration 
trace,  in  these  plots  every  digital  data  word  is  plotted. 

Quality  Control  Discussion 

Qualify  control  is  illustrated  by  the  group  of  plots  from  612  to  755.  HSP  is  the 
headspan  distance  from  the  mouth  target  to  the  head  target  and  is  illustrated  in  plots 
612  through  632.  It  should  be  a  constant  throughout  the  experimental  run.  Also  the 
value  should  be  identical  for  a  given  subject  for  all  runs.  Plots  612  to  623  grouped 
by  subject  show  that  the  maximum  difference  in  the  initial  value  of  HSP  is  only  5 
millimeters.  For  most  subjects  the  difference  is  Oven  less.  This  smalt  distance 
approaches  the  effective  resolution  of  the  photographic  system.  Review  of  HSP 
grouped  by  sled  peak  acceleration  category,  plots  624  to  632,  show  the  amount  of 
variation  developing  during  the  dynamic  event.  This  variation  remains  quite  small 
for  most  subjects  for  most  runs.  In  the  worst  case,  maximum  change  was  5%  of  the 
HSP  value.  As  such  it  confirms  the  following  presumptions: 

a)  the  response  of  the  head  was  two  dimensional  in  the  mid-sagittal  plane. 

b)  the  head  and  attached  mounts  acted  as  a  rigid  body. 
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c)  the  effects  of  error  introduced  by  the  photographic  data  system,  including 
scaling,  were  satisfactorily  controlled. 

d)  the  mounts  were  placed  on  the  head  in  the  same,  place  from  run  to  run  on 
the  same  subject. 

Review  of  the  headspan  plots  also  discloses  those  precise  data  points  where  photo- 
,  graphic  data  acquisition  was  incomplete.  This  is  indicated  by  gaps  in  the  individual 
^  traces.  It  should  be  noted  that  these  gaps  are  carried  throughout  all  variables  derived 
from  photographic  data.  No  interpolation  has  been  used  except  as  required  fc 
.  differentiation  routines  as  defined  in  the  signal  conditioning  section.  Table  T.  Even 
In  those  cases  the  gaps  are  maintained  in  the  finally  presented  differentiated  variable. 
The  reason  for  the  gaps  is  that  the  targets  were  not  acquired  by  the  camera,  usually 
because  some  portion  of  the  mount  straps  or  the  anatomy  occluded  the  target. 

A  comparison  of  angular  displacement  of  the  head),  PDH,  derived  from  photo- 
graphic  data;  and  angular  displacement  of  the  head,  RDH,  derived  from  the  rate 
;  gyroseope  Output  by  integration,  is  presented  In  plots  633  to  673.  By  definition 
they  should  be  identical.  Some  minor  differences  can  be  seen  by  reviewing  the 
plots.  This  comparison  vatdafes: 

/  i'  . .  i*’  ’  ,  '  '* 

a)  the  ability  of  the  sensor  data  system  to  supply  accurate  angular  displace¬ 
ment  information. 

b)  the  scaling  and  accuracy  of  both  sensor  and  photographic  data. 

c)  the  validity  of  the  measured  response  due  to  the  fact  that  the  measure¬ 
ment  systems  are  independent. 

d)  the  two  dimensional  nature  of  the  response. 

e)  the  minimal  error  due  to  parallax  in  the  photographic  data  system. 

Plots  674  through  755  are  presented  as  evidence  of  validation  of  the  coordinate 
system  transformation  procedure  applied  to  the  measured  data.  The  variables  in 
these  plots  are  components  of  the  acceleration  at  the  center  of  gravity  of  the  head 
derived  from  the  mouth  mount  acceleration  and  the  head  mount  acceleration  which 
are  independently  measured.  By  definition  the  measurements  should  be  identical. 

The  plots  illustrate  the  correspondence  and  limited  variation  between  the  two 
independent  measurements.  The  lack  of  variance  tends  to  verify: 


a)  coplanar  placement  of  the  accelerometers,  in  the  plane  of  induced  motion. 

b)  accurate  analog  recording  of  the  data. 

c)  accurate  analog  to  digital  conversion  of  the  data. 

d)  accurate  scaling  of  the  data.  . 

e)  accurate  execution  of  the  programs  which  transform  the  data. 


f)  stability  of  the  mounts  on  the  head  and  at  the  mouth. 

g)  motion  of  the  head,  and  attached  mounts  and  transducers  as  a  rigid  body. 

h)  accuracy  of  the  x-ray  anthropometry  methodology. 

In  summary  the  results  are  a  selected  sample  of  41  runs  illustrating  major  variables 
and  features  of  the  data.  The  29  variables  presented  represent  only  a  sample  of 
those  available.  The  runs  presented  represent  about  25%  of  those  eventually  avail¬ 
able. 

The  best  representation  of  the  results  is  the  plotted  data.  Much  more  information 
is  inherent  in  these  plots  and  the  digital  data  bank  from  which  they  were  derived 
than  could  be  described  in  a  single  monograph. 

Spectral  Characteristics  of  the  Sampled  Data  by  William  Anderson 

The  Fourier  coefficients  of  the  measured  photo  time  series  were  calculated.  In 
all  cases  the  amplitudes  of  the  Fourier  components  for  frequencies  greater  than 
30  HZ  were  at  noise  levels.  This  substantiates  the  assumption  that  a  photo  sampling 
rate  of  500  frames/sec  is  adequate,  and  also  indicates  that  a  30  HZ  low-pass  filter 
would  be  suitable  for  removing  noise  frcm  the  signal. 

The  Fourier  coefficients  of  the  measured  rate  gyro  time  series  were  also  calcu¬ 
lated.  The  rate  gyro  signals  normally  had  Fourier  components  at  noise  levels  for 
frequencies  larger  than  4Q  HZ.  Consequently,  a  50  HZ  low-pass  filter  was  used  to 
remove  noise  from  the  signals.  The  spectral  analysis  of  the  data  series  also  demon¬ 
strated  that  the  frequency  response  of  the  rate  gyres  (64  Hz,  +4  Hz)  was  adequate. 


The  Fourier  coefficients  of  the  measured  accelerometer  time  series  were  calcu¬ 
lated.  The  frequencies  at  which  the  amplitudes  dropped  into  the  noise  regions  were 
dependent  on  the  peak  sled  acceleration  and  the  amount  of  mount  slippage,  if  any, 
that  was  present. 

For  some  data  it  was  suspected  that  the  frequency  response  of  the  accelerometer 
was  attenuating  the  signal.  The  accelerometer  calibrations  yielded  a  range  of 
natural  frequencies  from  M2  Hz  to  165  Hz,  and  a  range  of  damping  ratios  from  0  5 
to  0.8.  However,  when  the  dafa  was  compensated  for  the  accelerometer  attenuation 
no  significant  change  was  observed.  Thus,  it  may  be  concluded  that  the  frequency 
characteristics  of  the  measured  accelerometer  data  have  been  preserved. 

It  should  be  mentioned  here  that,  except  for  sled  acceleration,  the  only  filtering 
done  was  to  smooth  data  enough  to  calculate  time  derivatives. 

Volunteer  Results  Discussion 

The  two  hundred  and  thirty-six  human  exposures  among  the  total  volunteer  group 
of  17  subjects  ware  conducted  without  injury.  Complete  records  of  the  examinations 
at  time  of  qualification,  at  time  of  departure  from  the  program  and  at  time  of  follow¬ 
up  one  or  more  years  after  cessation  of  the  experiments  disclosed  no.  pathological 
changes.  Examination  before  and  after  each  run  are  on  record  and  revealed  only 
minor  surface  abrasion  of  a  few  areas  of  skin  covered  by  the  restraint  system.  In 
every  case  the  subject  continued  his  usual  duty  and  off-duty  activities  after  the  post¬ 
run  examination. 
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RUN  DISTRIBUTION  BY  VOLUNTEER  SITTING  HEIGHT 
AND  SLED  ACCELERATION  CATEGORY 


Percent!  »e  Group 
12th  -  30th 


47th  -60th 


87th  -  96th 


Subject 

Category 

Actual 

No. 

Acc 

.  Onset 

Acc. 

Onset 

017 

6 

300 

5.6 

340 
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7 
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7 
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7 
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8.4 

440 

003 
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9.5 
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9 
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330 

004 

9 

500 

8.9 

630 

004 

9 
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7 
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7.1 
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on 

8 

300 

8.3 

300 
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8 
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8 

300 

8.5 

350 

Oil 

9 

300 

8.9 
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9 

300 
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RUN  DISTRIBUTION  BY  SUBJECT  AND  SLED  PEAK  ACCELERATION 


12th  to  30th  Percentile 


Subject  No. 

. 

Actual 

Percentile 

Acc. 

Onset 

013 
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47th  to  60th  Percentile 
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5.5 

250 

7.4 

450 
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5.9 

330 

6.1 

620 
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5.8 
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8.9 
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5.7 
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7.2 
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RUN  DISTRIBUTION  BY  SUBJECT  AND  SLED  PEAK  ACCELERATION 


Subject  No. 


87th  to  96th  Percentile 
Percentile 


Ac 

Acc. 


008 

007 

010 


95%  ... 

5.5 

96% 

6.8 

8.5 

9.6 

96% 

7.2 

8.9 

6.0 

7.1 

TABLE  R-2 


REPEATABILITY  PLOTS 


Subject 


Actual  Acceleration 


004 

004 

010 

016 


8.9  &  9.0 

5.9  &  6.1 
7.2  &  7.1 
8.4  &  8.4 


TABLE  R-3 


Actual  Rate  of  Onset 

630  &  550 
330  &  620 
320  &  510 
440  &  340 
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CONCLUSIONS 


The  considerable  theoretical/  experimental/  data  processing/  and  analytical  effort 
described  has  resulted  in  development  of  a  detailed  and  validated  set  of  data  describ¬ 
ing  the  complete  two  dimensional  kinematic  response  of  the  human  head  relative  to 
the  first  thoracic  vertebra  and  selected  reference  coordinate  systems.  The  described 
responses  are  for  41  -Gy  impact  acceleration  experiments  on  12  human  subject? 
with  torso  and  pelvic  restraint. 

From  careful  examination  of  the  results/  it  would  appear  that: 

a)  the  response  of  the  unrestrained  human  head  and  neck  to  -Gx  impact 
acceleration  is  two  dimensional  in  the  mid-sagittal  plane. 

b)  the  response  is  characteristic/  and  repeatable. 

c)  simultaneous  photographic  and  sensor  measurement  of  the  response  is 
possible,  thus  permitting  cross  validation. 

d)  adequate  presentation  of  the  response  requires  detailed  plotting  of  the 
variables  as  a  function  of  time.  Tabular  presentation  of  variable  peak  magnitudes 
is  not  an  adequate  way  of  presenting  the  complexities  of  human  dynamic  response. 

e)  the  utilization  of  the  theoretical  mechanics  of  rigid  body  motion  is  a 
valid  experimental  design  basis  for  measuring  the  response  of  human  head  and  neck 
to  impact  acceleration. 

From  the  methods  of  measurement  it  is  further  concluded  that: 

a)  the  use  of  x-ray  anthropometry  permits  X,  Y  and  Z  coordinate  description 
of  the  position  of  the  center  of  gravity  of  the  head  within  a  head  anatomical 
coordinate  system. 

b)  the  use  of  x-ray  anthropometry  and  the  center  of  gravity  determination 
makes  it  possible  to  express  the  data  at  the  center  of  gravity  of  the  head  rather  than 
at  the  anatomical  mounts.  This  permits  description  of  the  response  in  a  way  which 
is  independent  of: .  anatomical  variations  from  subject  to  subject;  possibly  from 
.'pecies  to  species;  location  of  instrumentation  mounting. 

c)  the  sampling  rates  of  the  photographic  data  system  and  the  frequency 
response  of  the  sensor  data  system  were  sufficiently  broad  band  to  measure  completely 
the  response,  in  the  reported  experiments. 


d)  the  automatic  data  acquisition,  reduction,  and  analysis  of  the  measure¬ 
ments  makes  it  possible  to  conduct  the  large  number  of  experiments  of  this  type 
necessary  to  achieve  statistical  significance.  The  computer-drawn  plots  of  the  data 
illustrate  the  success  of  the  effort , '  However,  the  variables  presented  in  the  report 
are  only  a  partial  sampling  of  what  is  presently  available.  Exhaustive  and  effective 
use  of  the  data  requires  large  scale  digital  cc  ;rer  interaction  with  the  data  base 
on  a  continuing  basis.  The  specific  interaction  should  be  determined  by  the  specific 
problem  posed. 

Certain  other  conclusions  can  be  drawn  regarding  the  utilization  of  the  data, 
which  can  be  used  for:  (a)  the  evaluation  of  any  proposed  two  dimensional  analog 
of  the  human  head  or  neck  up  t©  9.9G  peak  vehicle  acceleration,  in  the  -Gx 
direction,  and  up  to  800G/sec  rate  of  onset  of  acceleration;  (b)  comparison  of  the 
response  of  non-human  primates  to  the  human  response,  which  establishes  a  means 
for  calibrating  non-human  impact  acceleration  experiments;  (c)  an  exact  measure 
of  the  inertial  input  to  the  various  physiological  system  of  primates,  therefore 
serving  as  a  basis  to  which  to  relate  physiological  responses  to  impact  acceleration. 

Finally,  the  measurement  of  the  kinematics,  numerical  documentation  of  the 
anatomy  by  x-ray,  and  measurement  of  the  center  of  gravity  of  the  head  in  the  head 
anatomy  make  it  feasible  to  develop  a  valid  dynqmic  model  of  the  human  head  and 
neck. 
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